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A Graphical Demonstration of White Light Interference 
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University of Cincinnati, Cincinnati, Ohio 
(Received January 23, 1952) 


A teaching device is described which can be used in lectures on the interference of white light. 
The wave forms of four selected components of white light, in the violet, green, yellow, and red 
regions, are shown one above the other on a cardboard chart, inked in their respective colors. 
A sheet of clear plastic with an identical set of waves drawn on it is placed over the chart so 
that the two sets are in register. The combination represents a beam of white light as it arrives 
at the retina. A retardation of part of the beam is illustrated by displacing the plastic sheet 
laterally. Since the sine waves differ among themselves in wavelength, the various phase 
changes which result among the four components are shown quite clearly. 

The principle has been embodied in slide rule form with the wave trains appearing at a large 
plastic window. Five smaller windows indicate, respectively, path difference in millimicrons, 
dominant wavelength and purity of the resulting light when C.I.E. I!luminant C is used, and 
relative energies and resultant hue for the four components shown. 


INTRODUCTION 


ELECTIVE interference among the compo- 

nents of a beam of polychromatic light 
may be achieved in a number of ways, including 
thin film interference, slit and grating phenom- 
ena, birefringence of polarized light, etc. In 
elementary lectures on light and color, it is 
difficult to present a satisfactory picture of 
these effects in terms of wave phenomena 
because of limitations of student background. 
A graphical means to this end has been evolved 
which appears quite economical of background 
requirements. 

The point is made, in the lecture, that the 
various interference phenomena which result in 
colored light are but separate means to a common 
end, that of causing a part of the beam of white 


* Now at Research Laboratory, American Optical Com- 
pany, Stamford, Connecticut. 


light to traverse a greater optical path than the 
rest of the beam before arriving at the screen or, 
ultimately, the retina. Phase changes among 
various components at the destination are then 
illustrated by the device to be described. 


THE DEMONSTRATION 


In its simplest form the display consists of a 
rectangular sheet of white cardboard of conven- 
ient viewing size and a sheet of transparent 
plastic of the same size. The wave forms of four 
selected components of white light are shown 
one above the other on the cardboard sheet, 
extending across the width thereof. The sine 
wave trains are representative of the violet, 
green, yellow, and red spectral regions, and they 
are inked in their respective colors. Their 
relative wavelengths conform with actual rep- 
resentative wavelengths in these regions. An 
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Fic. 1. Interference slide rule shown for optical path difference of 530 mu corresponding to one complete “green”’ 
wave. Window A shows phase relationships among four selected components. Scales B through F are labeled as follows: 
B—Optical Path Difference in Millimicrons 
C—Dominant Wavelength for C.I.E. Illuminant C 
D—Percent Purity for C.I.E. Illuminant C 
E—Relative Energies for Four Components 
F—Resultant Hue for Four Components. 


Scale F is left blank in the figure. The scale indices are indicated by arrowheads. 


identical set is shown on the plastic sheet such 
that it may be superimposed upon the first set. 

When shown in registration, the two charts 
represent a beam of white light arriving coher- 
ently at the retina. A retardation of part of the 
beam is demonstrated by displacing the sheets 
laterally with respect to one another, whereupon 
the various phase changes become apparent. 
Since the four pairs of wave trains differ in 
wavelength from pair to pair, a retardation of 
one complete wave in the violet region, for 
example, will correspond to retardations of 
partial wavelengths in the other regions. This 
would amount to almost complete destructive 
interference of the ‘‘red’”” wave whose wavelength 
is nearly twice that of the violet. 


REFINEMENTS 


A more refined demonstration has evolved 
from this simpler form, offering the advantages 
of compactness and a more thorough treatment 
of the physical problem. The original cardboard 
chart is prepared as the slide of a slide rule with 
the plastic sheet mounted as a window in the 
body, as in Fig. 1. For various settings of the 
slide, the phase relationships are visible at the 
large window labeled A in the figure. 

A sliding scale (B, Fig. 1) indicates the optical 
path difference corresponding to each position 
of the slide. The scale is drawn on the slide and 
is visible through a plastic window set into the 
body. The scale index is marked on the window 
and is denoted by arrowheads at the top and 
bottom. 


















GRAPHICAL DEMONSTRATION OF WHITE LIGHT INTERFERENCE 


In addition to the wave forms and the optical 
path difference scale, certain other information is 
provided which may be of interest in more 
detailed studies of interference. color. This 
information is embodied in four additional 
scales. As with Scale B, these are drawn on the 
slide and indexed at their respective windows, 
the indices being aligned vertically with that of 
Scale B. On the actual slide rule, the five small 
scales indicated by letters in Fig. 1 are labeled 
as follows: 


B—Optical Path Difference in Millimicrons 

C—Dominant Wavelength for C.1.E. Illu- 
minant C 

D—Percent Purity for C.I.E. Illuminant C 

E—Relative Energies for Four Components 

F—Resultant Hue for Four Components. 


Scales E and F refer to two hypothetical 
beams of light whose original spectral composi- 
tions, before interaction, include radiant energy 
only at the four wavelengths indicated in Scale A. 
Scales C and D refer to the interference effects 
in two equal beams of light whose spectral 
composition corresponds to C.I.E. Illuminant 
C. This standard light source has a continuous 
spectral energy distribution as contrasted with 
the line spectrum in the previous case and is 
similar in appearance to daylight. 


THE E AND F SCALES 


At any setting of the slide, inspection of 
Scale A will reveal to some extent the relative 
amplitudes which may be expected among the 
four components as a result of interference. 
For the position indicated in Fig. 1, for example, 
corresponding to an optical path difference of 
530 mu, the resulting amplitudes are seen to be 
at a maximum in the green and to fall off in the 
order of yellow, red, and violet due to the increas- 
ing phase differences. The progression of ampli- 
tudes is based upon visual estimates of the 
algebraic sum of each pair of wave trains. These 
estimates at best cannot yield a precise evalua- 
tion of relative amplitudes. Moreover, since it 
is energy, or amplitude squared, rather than 
amplitude which is perceived by the eye as a 
measure of brightness, the relationships among 
the four resultant amplitudes are of no direct 
psychophysical interest. 
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It is the purpose of Scale E to afford a more 
explicit interpretation of the various phase 
relationships. This scale is labeled ‘Relative 
Energies for Four Components” and shows four 
trigonometric curves representing energy as a 
function of optical path difference. Each curve 
corresponds to one of the component radiations 
of Scale A and, on the slide rule, is inked in 
its corresponding color for identification pur- 
poses. Assuming that the two interacting beams 
are identical and that equal energies are repre- 
sented at the four wavelengths, it is clear that 
for zero path difference the relative energies of 
the four components after interaction will be at 
100 percent of their maximum values. As path 
difference is increased from zero, the four 
resultant energies begin to fall off unequally 
until the first minimum, in the violet, is reached. 
Thereafter, the remaining curves in turn reach 
their minima and begin to rise. Those portions 
of the curves between zero path difference and 
530 my are not visible in the figure. 


The trigonometric curves of Scale E are of the 
form 


I <cos?(A/2), 


which has the following significance. We assume 
two beams of monochromatic light which are 
identical in phase, energy, and wavelength \. 
If there should be imposed upon the beams a 
geometrical path difference of d my in a medium 
whose index of refraction is m, the optical path 
difference between beams will be md and the 
phase difference A will be given by A=2znd/). 
When the beams are recombined the relative 
intensity J as a function of A is indicated by 
the expression above. The intensity will be a 
maximum when the wave trains coincide and 
vice versa. 

For each of the four values of \ in Scale A, 
a curve has been plotted in Scale E. The wave- 
lengths of the curves are adjusted to accord 
with the units of md in Scale B. 

An obvious application of Scale E is the 
preparation of a Scale F to show the resultant 
hue at each position of the slide. The hue at 
each value of optical path difference is derived 
from the relative energies at that position. In 
preparing Scale F it would be assumed, of course, 
that the original light contained equal amounts 
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of the four selected radiations only. Scale F is 
labeled ‘Resultant Hue for Four Components” 
and in the original model was prepared with 
colored pencil. No attempt was made to repro- 
duce saturation or relative brightness. The 
scale is left blank in the reproduction. 


THE C AND D SCALES 


The usefulness of the slide rule is not limited 
to the case where the beam of light is composed 
of only four spectral lines. In fact, the value of 
the four-line spectrum is purely pedagogic and 
it was desired to generalize the treatment by 
including the case of the continuous spectrum. 
However, the interference colors in a beam of 
light of a continuous spectral nature will not 
necessarily be the same as those arising from 
light of a discontinuous nature even though 
both sources may appear white to the eye. Hence, 
neither Scale E nor Scale F will be applicable to 
the continuous spectrum. Furthermore, white 
light of continuous spectral composition does 
not have a unique meaning. There is a wide 
variety of spectral energy curves representing 
sources which will appear essentially white 
to the eye, and there is no @ priori reason for 
these white sources to yield identical colors 
under identical conditions of interference. 

It was therefore necessary to decide upon some 
specific source of white light of general laboratory 
interest, such as C.I.E. Illuminant C, and to 
have information describing the interference 
colors as a function of A for two equal beams of 
light originating in this source. The interference 
color information is available in the technical 
literature and will presently be dealt with in 
some detail. 

It was also necessary to select some means of 
color specification which would be more specific 
than that of Scales E and F. From the point of 
view of the physicist, a numerical system of 
color specification is to be preferred, and for 
that purpose several excellent conventions are 
available. Of these, the designation “dominant 
wavelength and excitation purity” is perhaps 
most closely akin to the subjective sensory 
aspects of a given color stimulus. This mode of 
designation more than compensates in preciseness 
for what it lacks by not depicting the colors 
visually as does Scale F. 


CHROMATICITY DIAGRAMS 


Interference colors are seen ‘“‘by their own 
light” as it were and fall into the category of 
self-luminous colors or colored light sources. 
They are to be distinguished from non-self- 
luminous or surface colors which are apt to be 
more subject to changes in appearance due to 
changes in background and surrounding illu- 
mination. 

A complete physical description of a colored 
light requires a specification in terms of three 
independent variables such as dominant wave- 
length, purity, and luminance. By relating these 
variables to three independent coordinates in a 
Euclidean three-space, it is possible to represent 
lights of all colors in a so-called ‘‘color space” 
wherein each point designates a_ particular 
color. In certain cases of color measurement, 
dominant wavelength and purity are often of 
more particular interest than luminance or 
luminance differences. In such a case the 
luminance aspect of color may often be treated 
separately or even ignored entirely. This is 
generally true in the study of interference colors 
where dominant wavelength and purity are 
treated jointly under the heading “‘chromaticity”’ 
and where luminance considerations are of 
secondary importance. 

The entire gamut of chromaticities is readily 
represented by a plane figure in the form of a 
graph, each point of which may be specified by 
only two variables, such as dominant wavelength 
and purity. Such a graph is known as a chroma- 
ticity diagram and may be thought of asa plane 
in the color space perpendicular to the luminance 
axis. Each point in the chromaticity plane may 
be thought of as representing an entire family 
of color sensations which differ among themselves 
only in luminance. 

There is considerable latitude in the choice of 
coordinates when a chromaticity diagram is to be 
prepared. The respective coordinates need not 
represent dominant wavelength and purity 
explicitly, and there are indeed certain advantages 
in another treatment which has been widely 
adopted. This method makes use of the principle 
that a given sensation of chromaticity may be 
reproduced by the proper mixture of any three 
primary light sources which are mutually 
independent, i.e., no one of which may be 
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matched in color by a mixture of the other two. 
If it is agreed to ignore luminance differences 
between the given color and the mixture of 
primaries, then the absolute amounts of the 
matching primaries may be ignored in favor of 
two ratios among the three quantities, reducing 
the number of variables from three to two. 
These ratios are usually chosen as the respective 
amounts of two of the primaries relative to the 
total amount, expressed decimally. The ratio 
expressing the amount of the third primary per 
total amount is superfluous since the three 
ratios sum to unity and hence only two of them 
are independent. 

It is thus possible to construct a graph in 
which one axis refers, for example, to the 
relative amount of green in the matching primary 
mixture and the other to the relative amount of 
red. Each color sensation would then be repre- 
sented by a point in the plane, and the proximity 
of this point to the origin of axes would indicate 
the relative amount of blue primary in the 
matching mixture. While the location of each 
point is expressed most simply in terms of the 
coordinate axes, it is possible to find a coordinate 
transformation so that each point is uniquely 
identified with a particular dominant wavelength 
and purity. 

The disposition of chromaticities throughout 
the diagram will depend upon the choice of 
primaries. The position of a certain color sensa- 
tion on one chromaticity diagram will generally 
differ from its position on another diagram whose 
axes are based on a different choice of primaries. 
It is therefore possible to devise as many 
different chromaticity diagrams as there are 
choices of primary sources. 


THE C.IL.E. CHROMATICITY DIAGRAM 


The chromaticity diagram adopted by the 
International Commission on Illumination! in 
1931 has become the most prominent means of 
communicating chromaticity data between 
workers in color measurement. This diagram is 
the product of an elaborate system of color 
specification standardized by the C.I.E. The 

1 Commission Internationale de l’Eclairage, Proceedings 


of the Eighth Session, Cambridge, England, 1931. 


2 Abbreviated ‘‘C.I.E.” here and abroad, from the French 
translation. ’ 


VALUES OF y 








VALUES OF x 


Fic. 2. The C.I.E. chromaticity diagram showing inter- 
ference color reflected by a thin film on a glass surface as 
optical path difference increases from zero to 1200 mu. All 
wavelengths are in mu. 


system is amply described in the literature*-5 
and will be discussed here only to the extent 
that it pertains to the derivation of Scales C 
and D. 

The data from which these scales were derived 
are included in the C.I.E. chromaticity diagram 
of Fig. 2. The x and y axes denote, respectively, 
the amounts of red and green C.I.E. primaries, 
relative to the total, which are needed to match 
a particular color. The convoluted curve contains 
the interference color data used in preparing 
Scales C and D. This curve shows the chromatic- 
ities of interference colors as a function of 
optical path difference. 

In addition there are depicted certain colors of 
especial interest. The point labeled C denotes the 
position of C.I.E. Illuminant C whose coordinates 
are x =0.310 and y=0.316. The horseshoe-shaped 
part of the enclosing curve is the locus of the 
spectral colors, the wavelengths of which are 
indicated at ten-my intervals. The straight line 
connecting the ends of the horseshoe is the 
boundary of the family of purple colors; wave- 
lengths have no direct meaning here. This line 


® Donald J. Lovell, Am. J. Phys. 18, 104 (1950). 

4D. B. Judd, J. Opt. Soc. Am. 23, 359 (1933). 

5 A. C. Hardy, Handbook of Colorimetry (The Technology 
Press, Cambridge, Massachusetts, 1936). 
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together with lines J and K denote a triangle 
which encloses the family of purple colors. 

There is one feature of the C.I.E. chromaticity 
diagram which makes it of tremendous value in 
predicting the results of additive color mixtures. 
If we designate any two colors by their positions 
on the chromaticity diagram and connect these 
positions by a straight line, this line will be the 
locus of all possible additive mixtures of the 
two colors. The position of a given mixture along 
the line will depend upon the relative amounts of 
the two components. The locus of the purples is 
nothing more than a mixture line for the violet 
and red limits of the visible spectrum. The 
entire enclosure therefore represents the locus 
of the most highly saturated colors which have 
physical meaning. Within the enclosure lie all 
real colors and outside is the realm of super- 
saturated or imaginary colors. 


DOMINANT WAVELENGTH AND EXCITATION 
PURITY 


As an application of the additivity property of 
the C.I.E. chromaticity diagram, we note that 
it is possible to draw a straight line from the 


Illuminant C point through any given point in 


the enclosure until it intersects either the 
spectrum locus or the locus of purples. In terms 
of actual color mixing, this means that it is 
always possible to find a spectral color or purple 
which when diluted with Illuminant C will 
match a given chromaticity. 

When a given color is matched by such a 
mixture, the wavelength of the spectral compo- 
nent of the mixture is defined as the dominant 
wavelength of the given color. When the chro- 
matic component of the match is a purple 
rather than a spectral color, the dominant 
wavelength is specified in terms of the comple- 
ment of the purple, or that spectral green whose 
mixture line with the purple passes through the 
Illuminant C point. The dominant wavelengths 
of the purple colors are indicated by adding 
the letter c after the numerical designation 
so as to distinguish them from their green 
complements. 

It should at this point be apparent that as 
we move outward from point C along a straight 
line toward the spectrum locus, the colors 


represented become more saturated while re- 
maining constant in dominant wavelength. The 
purity of a given color may be said to vary from 
zero at point C to 100 percent at the spectrum 
locus, assuming Illuminant C to be chosen as 
having zero saturation. More specifically, if we 
consider the distance from point C to the point 
representing a certain color, then the ratio of 
this distance to the total distance from point 
C to the spectrum locus is defined as the excita- 
tion purity of that color. Excitation purity is 
the psychophysical correlate of saturation and 
will be referred to here simply as purity. 

On the C.I.E. diagram it is therefore possible 
to transform the chromaticity coordinates x 
and y of any color to a specification in terms of 
dominant wavelength and purity, the trans- 
formation involving a fairly simple graphical 
construction and numerical calculation. These 
latter steps are reduced to a matter of simple 
interpolation through the use of a series of 
specially constructed charts appearing in Hardy’s 
Handbook of Colorimetry.® 


INTERFERENCE COLORS 


We have now formulated the basis for trans- 
forming the chromaticity data from x, y coor- 
dinates, in which they generally appear, to the 
more readily visualized dominant wavelength 
and purity. The remaining problem in the 
preparation of the slide rule is the one of selecting 
from the literature that set of interference color 
data which applies to the present case. Before 
this can be done, a distinction must be made 
between two general types of thin film inter- 
ference, one being a special case of the other. 

When a ray of light is reflected from the surface 
of a dense optical medium back into a less dense 
or rare one, a 180° phase change is introduced. 
This phenomenon is partly responsible for the 
phase difference in the beams of light reflected 
from the upper and lower surfaces of a thin film 
in air, the remainder of the phase differ- 
ence being accounted for by any optical path 
difference present. As the film thickness becomes 
vanishingly small, the interfering beams become 
exactly out of phase, since one of the reflected 
beams has suffered a phase reversal while the 
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other has not. The following cases may be listed 
as examples of net phase changes due to reflection: 


a. Reflection by a thin film in air. 

b. Reflection by a thin air film between 
dense media. 

c. Transmission by a thin film of inter- 
mediate refractive index bounded by 
one dense and one rare medium. 


On the other hand, in many cases of inter- 
ference either both of the reflected beams or 
neither of them happen to undergo phase 
reversals, in which case an infinitely thin film 
shows zero phase difference. This is illustrated 
by the following cases: 


d. Transmission by a thin film in air. 

e. Transmission by a thin air film between 
dense media. 

f. Reflection by a thin film of intermediate 
refractive index bounded by one dense 
and one rare medium. 


In the first group of examples, the sequence of 
interference colors as a function of optical path 
difference is distinctly different from the sequence 
resulting in the second group, the colors in the 
respective groups being complementary for a 
given optical path difference. The second group 
of phenomena, involving no net phase differences 
due to reflection, is of more general interest. 
It is under these conditions that the interference 
colors of thin films are very similar to those 
arising in birefringent and photoelastic pheno- 
mena and slit and grating phenomena. It is this 
group of phenomena rather than the first group 
which is treated in the present slide rule, for 
Scale A is so designed that for zero optical path 
difference the interfering beams are exactly in 
phase rather than out of phase. 


THE INTERFERENCE COLOR DATA 


Interference in thin films and birefringent 
media has been a topic of interest to many 
investigators since the time of Boyle. As a 
result of both early and contemporary investiga- 
tions, the periodical literature abounds in data 
specifying interference color as a function of 
film thickness or optical path difference. A few 
of these investigations are listed chronolog- 


ically.*"° The investigations are diverse in 
several respects and provide an adequate choice 
of data to fit various situations. Included in the 
data are the cases of phase changes arising from 
optical path difference alone as well as phase 
changes partly due to reflection by denser media. 
Most of the results are empirical, although recent 
papers have presented data derived on theoretical 
grounds. Interference colors are specified for 
various light sources as well. 

From the wide choice of data available, the 
theoretically derived results of Kubota!® appeared 
a sound choice on which to base Scales C and D. 
Kubota’s data satisfy the requirement of phase 
coincidence for zero optical path difference 
since they are computed for the case of reflection 
by a thin film on glass. Illuminant C is the 
assumed light source, which makes the data 
especially useful in the general case. 

As presented in Kubota’s paper, the data are 
in the form of the convoluted curve shown in 
Fig. 2. The original parameter, optical film 
thickness, has been replaced here by optical 
path difference. The former is the product of 
the film thickness and its index of refraction, 
and, since the light was assumed normally 
incident upon the film, these values had only to 
be doubled to yield values of the new parameter. 
The refractive index of the film is assumed 
constant at all wavelengths, thereby neglecting 
dispersion and widening the generality of the 
results for the present purpose. 

Readers familiar with the C.I.E. chromaticity 
diagram will recognize the progression of colors 
represented by the locus shown in Fig. 2. The 
continuity of the curve is of interest, as well as 
the repeated spectral orders. Starting at point 
C, we imagine an increasing thickness of film to 
be deposited on the glass surface as we trace out 
the locus of the interference colors. The light 
reflected from the surface takes on a pale yellow 
appearance at first, becoming more saturated 


as the locus approaches the spectral wavelength 


6 Lord Rayleigh, Scientific Papers me ty Press, 
— England, 1886), Vol. II, p. 4 


. V. Baud and W. D. Wright, J. aa Soc. Am. 20, 
381 (1930). 
— Buchwald, Ann. Physik 38, 5, 325-339 


(194 
, an Prins and J. J. M. Reesinck, Physica 12, 396 
1946 


10 Hiroshi Kubota, J. Opt. Soc. Am. 40, 146 (1950). 
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of 580 mu. As the optical path difference increases 
from 200 my to 240 my the reflected light 
becomes straw colored and then copper colored, 
having reached a maximum of saturation just 
after passing through 220 mu. A further increase 
in film thickness changes the color to reddish 
brown and then to purple at an optical path 
difference of about 260 my. A change through 
violet, pale blue, and a very pale green brings 
the locus close to point C again, nearly closing 
the first spectral order. A second order ensues 
and then part of a third, bringing the locus 
to the limit of Kubota’s investigation. It is well 
known that a further increase in film thickness 
will yield a few more spectral orders of decreasing 
saturation until the locus returns ultimately to 
point C whereupon no further interference color 
is seen. It is for this reason that interest is con- 
fined to thin films when interference color is 
concerned. 

By stereoptican means, the x and y coordinates 
of various points on the locus were found directly 
from Kubota’s figure. These coordinates were 
then transformed to values of dominant wave- 
length and purity through Hardy’s charts 
which are conveniently prepared for Illuminant 
C. The resulting triple column of figures, showing 
dominant wavelength and purity as a function 
of optical path difference, is represented graph- 
ically in Scales C and D, respectively. The slide 
rule designations for these scales are ‘‘Dominant 
Wavelength for C.I.E. Illuminant C” and 
“Percent Purity for C.I.E. Illuminant C.” 

In Scale C of Fig. 1, a discontinuous segment 
of the curve will be noticed. This segment, 
labeled ‘‘c,’’ designates those interference colors 
whose dominant wavelengths are the comple- 
ments of the spectral greens between 492 and 
567 mu. There are two such segments in Scale C, 
the first one not visible in the figure. Each 
segment corresponds to a segment of the Kubota 
locus which lies in the purple region of the 
chromaticity diagram. The discontinuities occur 
at the intersections of the locus and lines J and 
K and are responsible also for discontinuities in 
the slope of the purity curve of Scale D. It is 


recalled that as far as perception is concerned, 
the interference colors change in a perfectly 
continuous manner with changing path differ- 
ence. The discontinuities in the curves are 
merely curious features arising from the manner 


in which dominant wavelength and purity are 
defined here. 


DETAILS OF CONSTRUCTION 


The present slide rule, constructed of sheets 
of stiff cardboard and clear plastic, measures 
20X28 inches and is large enough to be seen 
throughout the average classroom. The optical 
path difference scale is so constructed that 4 inch 
corresponds to 100 mu. The six-inch range of 
the slide thus accounts for path differences up 
to 1200 mu, the limit of information available 
from Kubota’s paper. Using the above scaling, 
the four wave trains were drawn with wave- 
lengths corresponding to 420, 530, 580, and 
650 mu, as indicated in Fig. 1. 

The stationary set of wave trains and the 
indices for the five scales were drawn on a sheet 
of 0.010-inch clear acetate which was then 
mounted on the body of the slide rule directly 
over the slide. A sheet of 0.060-inch clear 
Plexiglas mounted over the acetate serves to 
protect the figures on it and provides rigidity. 
A sheet of heavy white paper with appropriate 
windows cut in it serves as the face of the slide 
rule and provides the surface fo- the labeling 
of the various windows. 

An adaptation which is here suggested is the 
preparation of the slide rule in the form of a 
pair of color slides for projection purposes. In 
this form the demonstration would be suitable 
for large audiences. 
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Two ballistics problems that may arise in future transportation have been formulated and 
solved in detail. The first problem is that of a vehicle moving along a great circle of the earth; 
the second problem is that of a rocket moving radially away from the earth. In both cases the 
magnitude of the resultant acceleration applied to the body has been set equal to a constant. 
The value of this acceleration will be determined ultimately by the maximum acceleration that 
human beings or material objects can withstand. 


PRIMARY function of transportation is to 

convey human beings rapidly from place 
to place with a minimum amount of discomfort. 
In general, for maximum transportation eff- 
ciency, the highest possible velocity would be 
maintained at all times. Nature, as formulated 
in nonrelativistic physics, prescribes no upper 
limit to possible velocity. On the other hand, 
nature imposes a rigid upper limit to change of 
velocity, or acceleration. Consider some common 
examples illustrating this effect of inertia: 


1. The maximum acceleration of an elevator 
is determined by the amount of physical discom- 
fort that the passengers will tolerate. 

2. An airplane requires a certain minimum 
amount of time and space to accelerate to its 
take-off speed or to decelerate to a stop. 

3. A jet plane or a rocket is limited in its 
minimum turning radius and maximum thrust. 

4. Maximum muzzle velocity in ordnance is 
determined by the material strength of the pro- 
jectile and of the weapon. 

Assume that the following function and its first 
derivative are single-valued and continuous: 


s=s(t), _ 
$= 


when t=0, (1) 


where the dot signifies differentiation of distance 
s with respect to time ¢. Write the identities 


_ ds (2) 
s=—, 
dt 


d’s 


s=— 


dt?’ 


Separate variables in Eq. (6) and integrate: 


2= f 2ds. 
0 


Substitute Eq. (2) in Eq. (7) and integrate: 


=f —— (8) 


Lm] 


Equations (7) and (8) are well known in me- 
chanics. Assume that & in Eq. (8) is constrained 
to the range 


OSISKi) Sc. (9) 


From the meaning of the definite integral, ¢ is 
minimized when § takes on its maximum value 
at all times. In the following problems, let the 
tangential acceleration § be maximized subject to 
the condition that the resultant of all accelera- 
tions applied to the body shall not exceed a con- 
stant c in magnitude. 


THEORY OF THE TRANSIT TUBE 


In order to minimize transportation time be- 
tween two points on a planet, subject to a 
maximum tolerable acceleration c, future man 
probably could not improve upon a nearly straight 





&4 WILLIAM 


PLANET 


PLANET 


(b) 


Fic. 1. (a) Coordinate system of transit tube. (b) Ac- 
celerations acting on transit tube vehicle moving with 
velocity $. 


tube in which his vehicle would travel. Future 
man will encounter difficulties if he should at- 
tempt to sink such a tube to great depths, for 
example, through the center of the earth. Cir- 
cumstances probably will force a compromise, 
that is, the construction of a curved tube along 
the circumference of the earth. Consider a body 
that starts from rest and proceeds along a great 
circle of the earth with a tangential acceleration 
§ such that the resultant of all accelerations 
applied to the body has a constant value c. In 
addition to the tangential acceleration, the exist- 
ing accelerations are those due to gravity and 
centrifugal force. The small, but not negligible, 
acceleration of Coriolis will not be considered. 
Consider a rectangular coordinate system, as 
shown in Fig. 1(a), attached to the vehicle, where 
j is normal to the earth. Formulated mathe- 
matically, see Fig. 1(b), the problem reduces to 
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a solution of the following differential equation 
with the associated boundary conditions: 


si+ (g—#/r)j=c, 


when t=0, (10) 


s=0 


where g and r represent, respectively, the average 
acceleration due to gravity, and the average 
radius of the planet. Equation (10) can be written 
in scalar form 

#+- (g—8/r)?=c?. (11) 


The dots in Eqs. (10) and (11) signify differentia- 
tion of distance s along the tube with respect to 
time ¢. For future reference, differentiate Eq. (11) 
with respect to time to obtain the relations 


dg “( ~) 
—_— = ——- a 
dt r = 
ag 28 ae 
5-4-2) 
r 


dj? +r 
With respect to a perpendicular to the surface of 
the earth, the direction of c, Fig. 1(b), can be 


(12) 


(13) 


+50 


VELOCITY ACCELERATION 
(mi/sec) (ft/sec?) 


DISTANCE 
(mi + 1000) 


20 
TIME (min) 


Fic. 2. Solution of the ballistics problems described in 
this paper. Solid lines represent the motion of the transit 
tube vehicle; dashed lines represent the motion of the 
vertical rocket. 
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represented by the relation 


cj g—#/r 
@=cos-'— =cos"! 
—Cc Cc 


(14) 


Differentiate Eq. (14) with respect to time, 
6=28/r. (15) 


Equation (15) states that the angular velocity of 
the resultant ¢ is twice the angular velocity of 
the vehicle with respect to the center of the earth. 


By the use of Eq. (6), a first integral of Eq.: 
(11) is obtained: 


ds 3 
ds [&—(g—8/r)*}! 
— 2sr—ds 


Y 8 
ml ee os 


rf. £ | g-#/r 
§ =>( sin —sin ). 
2 Cc Cc 


Write Eq. (11) in the form 
(ds/dt)? = c? — (g—8/r)?. (19) 


Separate variables and integrate to obtain the 
relation 


(16) 


(17) 


(18) 


(20) 


ss _——EEE— 
0 [(c—gt+8/r)(c+g—#/n) > 


Make the substitution 


$=[r(c+g) ]' cos¢ 
and transform Eq. (20) into the form 


r\t pr dp 
a 2 
eS o: (1—k? sin’p)! 
# 


o.=sin-f 1 


(22) 


where 


; 23 
saa a 


c+ 
ke =—., 


(24) 
26 


Equation (22) is in the form of Legendre’s 
standard elliptic integral of the first kind and can 
be evaluated readily by means of a table of 
elliptic functions, 


t=(r/2c)*[F(m/2, k) — F(x, k) ]. (25) 
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IE ae 


0.0 13.9 


aa 6 ah US 


19.0 22.0 25.0 


eC SE 


30.1 44.0 
(a) 


13.9 


START 


FINISH 


30.1 

(b) 
Fic. 3. Relative orientation of c. The numbers refer to 
time in minutes. 3(a). With respect to horizontal. 3(b). 


With respect to observer seated in a swinging chair free to 
orient with reference to c. 


RESULTS 


Equations (11), (18), and (25) represent the 
solution of the problem of minimizing transporta- 
tion time along a great circle of a planet subject 
to a constant resultant acceleration c. Figure 2 
represents a plot of these equations for a hypo- 
thetical trip determined by the following con- 
stants applicable to earth: c=40, g=32.16, 
r=20 891 200. (The units are English.) The 
solid lines of Fig. 2 represent the motion of a 
vehicle that accelerates to maximum velocity in 
the first 5000 miles and decelerates to a stop in 
the next 5000 miles. The curves are symmetrical 
about their mid-point. Such a hypothetical trip 
is equivalent to that from New York City, USA, 
to Darwin, Australia, via the Aleutian Islands. 
The direction of the constant resultant c, indi- 
cated by arrows in Fig. 3(a), is the direction that 
an observer in the vehicle would consider “down.” 
Assume that an observer is seated in a swinging 
chair free to orient with reference to the constant 
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(ft /sec*) 
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—(—gt)s 







® @=cos!—g/c. 


resultant acceleration c. He interprets the chang- 
ing direction of ¢ as rotation of the vehicle about 
his stationary seat. Figure 3(b) represents the 
sketches of Fig. 3(a) oriented in accordance with 
the viewpoint of the observer. As in Fig. 3(a), 
the arrows indicate the direction of the constant 
resultant c. The observer experiences an increase 
in his weight of about 25 percent; an individual 
who normally weighs 161 lb would seem to weigh 
200 Ib for the duration of the trip. 

Conditions at critical points that exist on the 
solid curves of Fig. 2, are illustrated in Fig. 3. 
General values of physical conditions that exist 
at nine such points are listed in Table 1. Many 
of these quantities, pertinent to the hypothetical 
trip of Fig. 2, are listed in convenient numerical 
units in Table II. The nine critical points are 
characterized by the following conditions: 

1. To avoid passenger distress caused by sud- 
den impact of force, the initial acceleration § is 
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TABLE I. Values for transit tube at critical points in Fig. 2. 
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applied gradually with a rise time of a few sec- 
onds. The initial point =0 represents conditions 
a moment after the entire acceleration § has been 
switched on, and before the vehicle has moved 
appreciably. The resultant acceleration ¢ is at an 
angle of about 36° with the perpendicular to the 
earth, Fig. 3(a). The observer interprets this 
condition as rotation of the vehicle through an 
angle of about 36°, Fig. 3(b). In terms of the 
hypothetical New York—Darwin trip, the vehicle 
begins to move out of the New York terminal 
toward the northwest with an initial tangential 
acceleration of about 24 ft/sec’. 

2. When ¢=9.2 min, the acceleration curve, 
Fig. 2(a), assumes an inflection, see Eq. (13). 
The centrifugal force at this point is } the force 
of gravity. This point occurs as the vehicle 
passes over Lake Superior en route to Canada. 

3. When ¢=13.9 min, the centrifugal force 
equals the force of gravity. Under these condi- 
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tions the vehicle would tend to circle the earth 
with neither a tendency to fall nor a tendency to 
recede to a greater distance from the earth. This 
point occurs as the vehicle passes a few hundred 
miles northeast of Edmonton, Canada. 

4. When ¢=19.0 min, the vehicle has acquired 
escape velocity; that is, should the vehicle be 
released from its constraint, it would recede to 
infinity with respect to the earth. This point 
occurs as the vehicle passes over the Alaska 
Peninsula. The vehicle has passed the northern- 
most point of its trajectory and has begun to 
turn south. 

5. When ¢=22.0 min, the vehicle, moving to 
the southwest a few hundred miles south of the 
Aleutian Islands, has acquired its maximum 
velocity. The tangential acceleration § has been 
reduced to zero and the acceleration curve, Fig. 
2(a), passes through an inflection point, see Eq. 
(13). Assume that subsequent to this point the 
vehicle begins to decelerate. This assumption is 
equivalent to the requirement that Eq. (12) be 
continuous. Equation (11) would still hold, how- 
ever, if the acceleration § were held at zero for 
some time subsequent to this point and prior to 
deceleration. At this point the resultant accelera- 
tion ¢ is directed radially away from the earth, 
Fig. 3(a). The observer interprets this condition 


as rotation of the vehicle through an angle of | 


180°, Fig. 3(b). 

6. When ¢=25.0 min, the distance s, see Table 
I, is independent of c. This point on the trajec- 
tory occurs approximately half-way between 
Midway Island and Tokyo, Japan. 

7. When t=30.1 min, the vehicle is a few 
hundred miles south of Guam. 

8. When t= 34.8 min, the vehicle leaves North 
New Guinea toward Australia. 


TABLE II. Numerical values for transit tube at 
critical points in Fig. 2. 


Point (mi) 
0.00 0 
2.83 730 
4.89 1842 
6.95 3684 
7.33 4940 
6.95 6216 
4.89 8058 


t ‘s s 
(min) (ft /sec?) (mi/sec) 

0.0 +23.7 

9.2 +33.8 
13.9 +40.0 
19.0 +23.7 
22.0 0.0 
25.0 —23.7 
30.1 —40.0 
34.8 — 33.8 2.83 9170 
44.0 —23.7 0.00 9900 
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Fic. 4. Geometry of the radial rocket. 


9. When ¢=44.0 min, the vehicle has ended its 
trip. The observer interprets conditions at this 
point as rotation of the vehicle through an angle 
of about 324°, Fig. 3(b). In practice, the tan- 
gential acceleration § would be switched off with 
a delay time of a few seconds to avoid passenger 
distress caused by sudden release of weight. 


THEORY OF THE VERTICAL ROCKET 


A problem similar to that of the transit tube 
involves a rocket that starts from rest, rises 
radially from a planet, and recedes to infinity 
against the force of gravity which diminishes as 
the inverse square of the distance from the center 
of the planet. In order to minimize transportation 
time, set the resultant acceleration equal to a 
maximum constant value c. Formulated mathe- 
matically (see Fig. 4) the problem reduces to the 
solution of the following differential equation 
with the associated boundary conditions: 


=P 


§+gr’s?=c, (26) 


| when t=0, 


s=0 


where g and r represent, respectively, the average 
acceleration due to gravity and the average 
radius of the planet. Integrate Eq. (26) and intro- 
duce the boundary conditions to obtain the 
relation 


432 =cs+gr*s-!—cr—gr. (27) 
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Equation (27) can be written in the form 


1/ds\? cs*—r(c+g)s?+¢r’s 
( ) = ; (28) 
2\ dt ? 


s2 





Separate variables in Eq. (28) and rearrange to 
obtain the relation 


(s—0)ds 





——. (29) 


P Qrrv3 
t= 
- ) (r —0)*[4(s—r)(s—ger)(s—0) ]! 


Equation (29) is in a form that can be integrated.! 
Substitute into Eq. (29) the relation 

s=rsin—¢ (30) 
to obtain the relation 


2r\? —ddo 
d= (-) —_—————_— ; 
c J sin’¢(1—k? sin’¢)? 


kR?=g/c. 


(31) 


where 
(32) 


Equation (31) can be written in the form 


2r\3 
dt= (=) E —k* sin’p)'— (1—k? sin’g)-! 
c 
1—k? sin‘ 
-—_———_— |w. (33) 
sin’¢(1—k? sing)! 


Integrate Eq. (33) to obtain the relation 


2r 3 ol 
f= (=) (1—k? sin’¢) d¢ 
c x/2 


t do 
j2 (1—R? sin’g)} 


. 7 


+cotds(1—# sin’)! (34) 


where 


¢$,=sin—(r/s)}. (35) 


Equation (34) can be evaluated by means of a 
table of elliptic functions.? Write Eq. (34) in the 


1Edwin P. Adams and R. L. Hippisley, Smithsonian 
Mathematical Formulae and Tables of Elliptic Functions 
(The Smithsonian Institution, Washington, D. C., 1939), 
first edition, p. 256. 

2G. W. Spenceley and R. M. Spenceley, Smithsonian 
Elliptic Functions Tables (The Smithsonian Institution, 
Washington, D. C., 1947), first edition. 


form 


Qr\? 
= (=) [x- Fis, )- E+E, 8) 
c 


[CY ~J-OO]}: & 


Equations (26), (27), and (36) represent the 
solution of the problem of minimizing transpor- 
tation time for a rocket in vertical ascent subject 
to a constant maximum acceleration c. These 
equations are plotted in dashed lines, Fig. 2, for 
a hypothetical trip determined by the following 
constants applicable to earth: c=40, g=32.16, 
r=20 891 200. As in the previous problem, an 
observer with an initial weight of 161 lb would 
appear to weigh 200 lb for all time during which 
the rocket was accelerating. Only two points are 
considered critical, that is, the inflection point 
of the acceleration curve, Fig. 2(a), and the point 
where escape velocity has been attained. 

The point of escape can be determined from 
the relation 


$5—1 = 2gr?s-, 


(37) 


Eliminate § between Eqs. (27) and (37) to obtain 
the relation 


s=(1+gce)r. (38) 


Substitute Eq. (38) into Eq. (37) to obtain the 
escape velocity 
2gerc\} 
in (—) . 
c+g 


Escape velocity is attained at the point t=27.3 
min. Of secondary interest is the distance s at 
which the curve of_Fig. 2(a) has an inflection 
point, 


ZZ] 


This inflection point occurs at the point t=17.9 
min, 


(39) 


CONCLUSIONS 


The two problems solved in this paper repre- 
sent special cases of a single generalized problem, 
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that is, the problem of a vehicle that moves along 
an arbitrary path away from a planet subject to 
a constant resultant acceleration c. In problems 
of this kind, the imposition of constant resultant 
acceleration seems to be a realistic engineering 
requirement. Consequently, the results presented 


herein might well reflect a solution of future 
transportation problems. The final engineering 
problem must contain Coriolis corrections and 
allow for rate of increase of acceleration.* 

3 The importance of time derivatives of acceleration in 


mechanical and transportation problems has been con- 
sidered by P. LeCorbeiller, Am. J. Phys. 13, 156 (1945). 


The Definitions of Magnetic Flux Density and Field Intensity 


J. V. HuGues* 
Queen’s University, Kingston, Ontario, Canada 
(Received August 9, 1951) 


The present theory of magnetism, in which it is treated as a phenomenon of current elec- 
tricity, suffers from objections which are pointed out. The meanings of B and H on an atomic 
scale are considered, and a method of picturing B and H, inside a ferromagnetic body, without 
mentioning magnetic poles or dipoles, is developed. It is shown that this method satisfactorily 
accounts for the demagnetizing field in certain simple cases and that it also makes clear why H 
occurs in the circuital theorem while B is used in Laplace’s equation. Rigorous definitions of B 
and H at any point in a magnetic medium are then developed. In the concluding section the 
advantages of the suggested method are summarized. 


I. INTRODUCTION 


T one time, all the definitions of magnetic 
quantities were based on the concept of an 
isolated pole. These definitions, including those 
of magnetic induction or flux density B and mag- 
netic field intensity H inside a solid ferromagnetic 
material, were clear-cut and precise.! The system, 
however, suffered from certain objections. More 
recent theory, which treats magnetism as a 
branch of current electricity, overcomes these 
objections, but as usually presented it suffers 
from the following defects: (a) The definitions 
of B and H at a point in a solid ferromagnetic 
body are given in a very imprecise form. (b) 
Largely as a result of (a), it is often difficult to 
tell why B occurs in some results and H in others. 
(c) In arrangements involving ferromagnetic 
materials, the distribution of H is often difficult 
to reconcile with a distribution due solely to 
currents. 

The aim of this paper is to give a theory free 
from these objections. The present definitions of 
B and H are considered, to illustrate the objec- 
tions. By considering the atomic model, a method 


* Now at Southwest Research Institute, San Antonio, 
Texas. 


1G. Burniston Brown, Proc. Phys. Soc. (London) 52, 
577 (1940). 


of overcoming these objections is found, and this 
leads to rigorous definitions of B and H at any 
point, without mention of magnetic poles. 


Units and Abbreviations 


Rationalized MKSC units are used. Symbol 
Bo=permeability of vacuum. The abbreviation 
CC is used to denote ‘‘circulating current.” 


Il. THE ATOMIC CIRCULATING CURRENT 
AND THE EQUIVALENT DIPOLE 


Magnetism is frequently considered as due to 
atomic CCs. This idea of an atomic CC needs 
clarification. Each atom is assumed to contain 
many electrons, revolving in orbits and spinning 
on their axes, which must give rise to magnetic 
effects similar to those due to currents. For most 
purposes the total magnetic effect of an atom can 
be considered the same as that of a current flow- 
ing in a closed loop. It is this equivalent current 
in a closed loop which is meant by the atomic CC. 

The magnetic effects of a CC are, at a distant 
point, the same as those of a magnetic dipole. 
This process of replacing a CC by an equivalent 
dipole is a useful mathematical device, but no 
physical significance should be attached to the 
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magnetic poles thus introduced. It should, more- 
over, be possible to define B and H without the 
necessity of introducing poles by this device. 


Ill. PRESENT DEFINITIONS OF FLUX DENSITY 
1. Flux Density B in Air 


There are three possible methods of defining B 
at a point in air, viz., (B1) in terms of the force 
experienced by a current element; (B2) in terms 
of the couple experienced by a small current loop; 
and (B3) in terms of the charge circulated (or 
emf induced) when the flux linkage with a small 
coil is altered by rotating the coil, or by with- 
drawing it to a point of known B, or by altering 
B to some known value. 

Method (Bi) is the commonest definition. 
Textbooks rarely use (B3) as a definition, since 
considerable work is usually done on magnetic 
fields before induced emf’s are considered. But it 
would require only minor rearrangement to re- 
verse this order and use (B3) as the definition 
of B. The outstanding advantage of this is that 
(B3) is the usual method of measuring B in 
practice, not only in air but also in ferromagnetic 
bodies, whereas the other two methods are 
of little use when ferromagnetic bodies are 
considered. 


2. Flux Density B Inside a Solid 
Ferromagnetic Body 


Modern theory is very weak on this definition. 
If the body is in the form of a bar, or an anchor 
ring, a coil wound on the body can be used to 
give an average value of B over the cross section, 
by (B3) above. This is adequate for most prac- 
tical purposes, and is as far as most textbooks go. 
The determination of B at any specified point in 
the body is usually not considered. Smythe? sug- 
gests that a spherical cavity can be excavated at 
the point, filled with a liquid of the same permea- 
bility as the solid, and that a measuring ap- 
paratus (e.g., a rotating coil) placed in the liquid. 
This could be used as a theoretical definition, but 
details of a method of determining when the 
match of permeabilities was achieved would be 
needed. 


2W. R. Smythe, Static and Dynamic Electricity (Mc- 
Graw-Hill Book Company, Inc., New York, 1939). 
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IV. PRESENT DEFINITIONS OF FIELD INTENSITY 
1. Magnetic Field Intensity H in Air 


For H at a point in air, there are two very 
differently worded, but almost identical, defini- 
tions, viz.: (H1) The field intensity is equal to the 
flux density at the point divided by yo: 


H=B/wo. 


(H2) The field intensity is given by the integra- 
tion of a specified function over all the currents 
causing the field 


H=(1/4r) Silds Xr ]/r*. 


The unit “‘amp-turns/m” for H arises from this 
definition. This is the total H if magnetic ma- 
terial is absent, but if it is present this calculation 
gives only that part of H due to the external 
currents. Since the portion of B due to external 
currents is given by yuo times the same integral, 
this definition is only a variation of (H1). 


2. Magnetic Field Intensity H inside a Solid 
Ferromagnetic Body 


In a ferromagnetic body, the portion of H due 
to the external currents is calculated by (H2) 
above. In addition, there is a portion of H due to 
the magnetism of the body, which is usually in 
the opposite direction to B inside the body, and 
is called the demagnetizing field. The currently 
accepted theory gives no method of calculating 
this demagnetizing field for the general case. 
Method (a) below expresses the demagnetizing 
field in terms of another unknown variable, while 
method (b) deals with a particular case. These 
methods are given to point out their deficiencies. 

(a) To every elemental volume of the body is 
assigned a magnetic moment, and the demag- 
netizing field is obtained by evaluating the field 
due to such a magnetic moment and integrating 
over the body. The connection of this process of 
assigning a magnetic moment to an elemental 
volume with the electric current theory of mag- 
netism is based on the possibility of replacing the 
atomic CCs by magnetic dipoles. The magnetic 
moment per unit volume is the intensity of 
magnetization M, a vector function of position. 

(a’) A variation of this method is to assign M 
as above and then obtain H at the point (total H, 
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not the demagnetizing field) from the equation® 


B=y.H+M. (1) 


(b) This method considers the special case of 
a cylindrical bar of ferromagnetic material, of 
finite length, uniformly magnetized in a direction 
parallel to the axis of the bar.‘ All the atoms are 
assumed to have CCs of the same magnitude, 
and their axes parallel. For one section through 
the bar, these atomic CCs are shown in light lines 
in Fig. 1, and the method suggests that these can 
be replaced by a single CC round the outside of 
the bar, shown in heavy line in Fig. 1. The entire 
bar can thus be replaced by an equivalent sole- 
noid. The number of amp-turns/m of this equiva- 
lent solenoid is chosen to give the same H at one 
point outside the bar as is actually given by the 
magnetized bar. Then H at other points outside 
the bar, and also the torque on the bar when 
placed across a weak magnetic field, are cor- 
rectly given by this equivalent solenoid. It is 
readily shown that by attributing equal and oppo- 
site poles to the ends of the bar, along with 
plausible assumptions about the action of poles, 
one can equally well determine both the external 
H and the torque. The demagnetizing field at any 
point inside the bar is now calculated from these 
poles, completely ignoring the equivalent solenoid. 

It will be seen that methods (a) and (a’) cannot 
give the demagnetizing field unless the value of 
M at every point of the body can be found. 
Method (b) only deals with a special case, so that 
none of the methods can find the demagnetizing 
field in the general case. But even accepting this 
fact, there are other objections to these defini- 
tions: (1) Both (a) and (b) find it necessary to 
introduce poles. In (a) the word “pole” may not 
be mentioned, but magnetic dipoles are intro- 
duced either explicitly or by obtaining the 
mathematics by analogy with dielectric theory, 
which is based directly on electric dipoles. In 
(b) the poles are only an alternative for the 
external field, but become an essential part of the 
calculation of the demagnetizing field. Method 
(a’) might be worded to avoid poles, but it 
reduces Eq. (1) to the status of a definition, 


3Or B=(H+M), depending on the definition of M. 

‘It is recognized that an elongated spheroid is more 
amenable to mathematical treatment, but for descriptive 
work the uniformly magnetized cylinder is more convenient. 


Fic. 1. Circulating currents in a magnetized bar. Light 
lines: atomic circulating currents; heavy line: principal 
circulating current; dotted lines: reverse circulating 
currents. 


thereby removing any possibility of deducing 
this equation. (2) If the bar of method (b) is per- 
manently magnetized, then H is in the same 
direction as B outside the bar, but in the opposite 
direction inside the bar, the change of direction 
of H occurring suddenly at the end faces of the 
bar. Such a sudden change of direction does not 
occur for H of the equivalent solenoid, and the 
method fails to explain this discrepancy. This 
objection is equally applicable to (a), since the 
methods have much in common. For this bar, 
method (a) would give M constant, and the 
Poisson distribution equivalent to the set of 
dipoles is then a pole on each end of the bar, 
equal and opposite, agreeing with method (b). 
The difference between the methods is that (a) 
immediately replaces the atomic CCs by poles (in 
the form of dipoles) and then sums the effects, 
while (b) attempts to sum the atomic CCs and 
then converts the result to a distribution of poles. 

It appears, therefore, that both methods pay 
lip service to the idea of CCs, but ultimately have 
to resort to the idea of poles; and both produce 
a result which they cannot explain in terms of 
currents only. 


Vv. A SUGGESTED ATOMIC INTERPRETATION OF 
B AND H: BASIC PRINCIPLES 


The above definitions show that B has physical 
significance (affecting forces on currents, induced 
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emf’s) but they do not explain why H is needed. 
The need for H arises when ferromagnetic ma- 
terials are studied. When an anchor ring of the 
material is wound with a toroidal primary and a 
small secondary coil, a change in primary current 
produces a change in B, the average of which over 
the cross section is measured by the secondary 
coil. Flux density B thus obtained varies in a 
nonlinear manner with the primary current, and 
an extension of the experiment shows that it is 
the amp-turns/m of the primary on which B 
depends. This quantity, the amp-turns/m of the 
primary, is given a name and symbol, H. Since 
for a toroid wound in air B is always proportional 
to H, as just defined, it is reasonable to assume 
that for cases in air where B varies from point 
to point, H also varies and is a field vector. 

It is usual to represent B by lines of flux drawn 
so that their density (lines/m? or webers/m?) 
gives the magnitude of B. It is also usual to draw 
lines of H with a similar convention, but in the 
following a different convention is used. Lines of 
H will be drawn so that their density (lines/m?) is 
Mo times the magnitude of H. This has the ad- 
vantage that, in air, lines of flux can be considered 
as either lines of B or lines of H. In a ferromag- 
netic body the convention will be maintained 
that the lines of H have a density yo times the 
magnitude of H, but lines of B and of H will no 
longer coincide. 

In this and the next section we must imagine 
that we can see the atomic CCs, like wire loops 
carrying currents, and see the lines of flux like 
pieces of string threading these CCs. After a clear 
picture has been obtained, by this process, of 
what represents B and what H, we can go on to 
consider how B and H can be measured when we 
cannot “see” lines of flux but have to rely on 
instrument readings. 

On this picture, the ferromagnetic body is re- 
placed by a set of atomic CCs in vacuum. If B 
and H could be observed on this scale, they would 
both vary rapidly from point to point, with lines 
of flux representing both B and H and with 
H=B/yo. But when we speak of B and H ina 
ferromagnetic body, we refer to averages over an 
area large compared to atomic dimensions, so 
that the local variations around the atoms do not 
appear in the results. If this averaging were 
carried out in a straightforward manner H would 
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still remain equal to B/yo. The fact that this 
relationship is not observed for a ferromagnetic 
body shows that the method of averaging to find 
H is not the obvious one, and we must try to find 
what method is used. To do this, we return to the 
case of the toroid on an anchor ring, for which we 
defined H as the amp-turns/m of the primary. 

For this special case, H is the same whether 
the anchor ring is made of ferromagnetic material 
or of air. Since the effect of a ferromagnetic body 
is due to the atomic CCs, the above fact suggests 
that in obtaining H from the lines of flux, we 
should count only selected lines of flux, which 
pass where the effect of the atomic CCs is least, 
which is in the gaps between atoms. The infer- 
ence is that in a block of ferromagnetic material 
of arbitrary shape, the value of H at any point 
is obtained by counting the number of lines of 
flux passing, in a certain small area, through gaps 
between the atoms in that area, and dividing by 
the area and by wo. 

The next simplest case is that of the short 
magnetized bar of method (b), Sec. IV, which we 
shall at first consider to be permanently mag- 
netized (no external currents). Assuming for 
simplicity that successive cross sections all have 
their atoms in line with those of Fig. 1, each such 
line of atoms will be similar to a line of coaxial 
coils. Some lines of flux will thread through the 
CC of each atom in turn along the line, giving 
continuous lines of flux in the general direction 
of B. These will be lines of B. But each atomic 
CC will also have some lines of flux which form 
closed loops round the current. These closed 
loops have the same general direction as B inside 
the atomic CCs, but the opposite direction out- 
side them, between the lines of atoms. As lines 
of H are observed between the atoms, they must 
be identified with these portions of closed loops 
between the lines of atoms. Thus, our assump- 
tion, i.e., that in obtaining H only lines of flux 
passing between atoms are to be counted, leads 
immediately to an explanation of the fact that 
H inside the material is in the opposite direc- 
tion to B. 

It must be realized that at this point in the 
discussion the ideas are only qualitative. As there 
is no means of knowing how many closed loops 
of flux should be associated with each atomic CC, 
no numerical value can be given for H even when 
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B and the geometry are fixed. The quantitative 
aspects of the problem are discussed in the next 
section. For the moment we shall assume that 
the total number of lines of flux due only to the 
atomic CCs is fixed and invariant, whether a 
current is flowing in an external solenoid or not. 
This is a reasonable assumption as long as the 
external current is not sufficient to cause realign- 
ment of the atomic CCs. 

If now a current in an external solenoid is 
applied, and gradually increased in the direction 
to assist B, more and more of the lines of flux 
will flow continuously through successive atoms, 
reducing the number of closed loops of flux and 
hence the reverse field intensity, until all the 
closed loops are absorbed (H=0).° Further in- 
crease of the current will force some continuous 
flux lines between the lines of atoms, giving H in 
the same direction as B. Figure 2 gives sketches 
of these three conditions, illustrating the way in 
which B and H are viewed on the ideas advanced 
above. 

It is interesting to compare the behavior of 
lines of B and of lines of H for such a bar when 
permanently magnetized. A line of B, followed 
from an arbitrary point in air, will proceed until 
it enters the bar; then, according to the above 
description, it will continue through the bar, 
threading its way through successive atoms, until it 
emerges from the bar and continues through air 
to its starting point.® Starting from the same 
point in air, the line of H follows the same flux 
line as the line of B, until it reaches the bar. The 
line of flux continues into the bar, but it is no 
longer called a line of H; the line of H stops short 
on the surface of the bar. Inside the bar are some 
closed lines of flux linking atomic CCs, and the 
portions of these between lines of atoms, where 
they are in the opposite direction to B, are 
counted as lines of H. In other words, lines of H 
outside the bar stop on the face of the bar and 
have no connection with the lines of H inside the 
bar, which in turn are not continuous but are 
only selected portions of complete, closed lines 
of flux. 


5 This statement relies on the assumption that the total 
number of lines of flux due to an atomic CC is fixed and 
independent of the current in the external solenoid. 

6 Or to near its starting point, as J. Slepian (Am. J. Phys. 
19, 87 (1951)) has shown. The difference does not affect our 
argument in any way. 


Fic. 2. Axial cross section through bar at various stages 
of magnetization. (a) No current in magnetizing solenoid; 
(b) solenoid current adjusted to give H=0 in bar; and 
(c) solenoid current greater than in (b). Solenoid current 
flows clockwise when viewed in direction of B. 

Six atomic currents are shown in cross section; a dot and 
a cross represent current flowing out of and into the paper, 
respectively. Full lines, with arrowheads, are lines of flux 
which are counted as lines of B; dotted lines are parts of 
lines of flux which are not counted for either H or B; and 
dashed lines are parts of lines of flux which are counted as 
lines of H. 
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Fic. 3. Axial cross section through cylindrical bar mag- 
net, showing symbols used in Sec. VI, Part 2. The principal 
circulating current, of (N72) amp-turns, is flowing like a 
current in a solenoid round the bar, counterclockwise when 
viewed in the direction of Hy. 


This picture explains the sudden change of 
direction of H at the face of the bar. In so doing, 
it makes it clear that H has little or no physical 
significance, but is introduced only because of 
the way we choose to specify the magnetic prop- 
erties of ferromagnetic materials. 


VI. A SUGGESTED ATOMIC INTERPRETATION OF 
B AND H: MATHEMATICAL RESULTS 


Since the portions of B and H due to external 
fields are readily calculated, in the following only 
By and Hy, the portions of B and H due to 
magnetism of the material, are considered. 
Earlier Hy was called the demagnetizing field. 
As already mentioned, when we speak of B and 
H in a solid body we refer to averages (taken in 
a special manner as explained above) over an 
area which contains thousands of atoms, but 
which is infinitesimal compared to the dimen- 
sions of the body. The next step is to progress 
from the atomic model to a picture with no 
smaller unit than these infinitesimal areas. Dur- 
ing this step a portion of the body with the cross 
section shown in Fig. 1 will be considered, and 
we shall continue to assume that we can see the 
lines of flux, and hence can count their number. 

For By, counting the lines of flux crossing an 
area is the same as counting the lines threading 
continuously through lines of atoms, since the 
closed loops of flux can add nothing to the count. 
In method (b) of Sec. IV the atomic CCs (light 
lines of Fig. 1) were replaced by a single CC 
round the outside of the bar (heavy line of 
Fig. 1), hereafter called the principal CC. This 
is equivalent to allowing the orbits of the atomic 
CCs to expand and distort until they touch at 
all points to form a complete mesh, while the 
value of the current decreases to maintain the 


number of continuous lines of flux constant. The 
number of amp-turns in the principal CC, which 
is chosen to give agreement with the experimen- 
tally observed B, is equal to the value for these 
expanded and distorted CCs. 

But the replacement of the atomic CCs by the 
principal CC gives a picture which excludes Hy. 
There should also be included reverse CCs, 
shown dotted in Fig. 1, in which the current is 
exactly equal and opposite to that in the princi- 
pal CC. In obtaining Hy by counting lines within 
an infinitesimal area, only lines passing through 
these reverse CCs must be counted. If now the 
orbits of the atomic CCs are allowed to expand 
and distort as explained above, with consequent 
decrease of the value of the principal CC and the 
reverse CCs, then the reverse CCs within the 
infinitesimal area must also expand (and possibly 
distort) to keep the number of lines through them 
constant. As a working hypothesis we shall as- 
sume that they meet to give a complete mesh, 
so that they can be replaced by a single reverse 
CC round the infinitesimal area. Consequently, 
Hy» is given by the number of lines passing 
through an infinitesimal area, divided by the 
area and by yo, when the principal CC flows round 
the outside of the bar and an equal reverse CC flows 
round the infinitesimal area. 

Since for both principal CC and reverse CC, 
adjacent turns are separated by the interatomic 
distance while their cross sections have dimen- 
sions vastly greater than this, both principal CC 
and reverse CC can be considered as flowing in 
tightly wound solenoids. The amp-turns for both 
these currents will be denoted by (Nz). With this 
picture we can now consider some mathematical 
results. 


1. Anchor Ring Sample, as in B-H 
Measurement 


Both the principal CC and reverse CC flow in 
toroids, so that each produces a field intensity of 
(Ni)/l, where J is the mean length of the toroid. 
As these intensities are in opposite directions, 
Hy=0. 

The starting point of the present discussion 
was the fact that H in such a sample was just 
that due to external currents, so that Hy=0. 
The fact that our modified picture of CCs has 
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produced this result is a justification of the work- 
ing hypothesis given above. 


2. Short, Uniformly Magnetized Bar 


An axial cross section, explaining the symbols, 
is given in Fig. 3. The principal CC is a solenoid 
of finite length, giving a field intensity Hp, where 


Hp=(Ni/21)(cosa+cosy). 


The reverse CC flows in a solenoid of infinitesimal 


cross section, and so gives a field intensity 
Her=(N2)/l. Hence, 


Huy =Her—Hp=(Ni/21){(1—cosa)+(1—cosy)}. 


This result is precisely that which would be ob- 
tained by assuming a uniform density of pole 
strength on the two faces, so that it agrees with 
the result given by the older method. 


3. The Circuital Relation 


The circuital relation, im air, can be stated in 
two forms: 


FSH-ds= i, (a) FB-ds=modi, (b) 


where 77 is the sum of all currents linked by the 
path of integration. If magnetic media are pres- 
ent anywhere on the path of integration, both 
results must still hold, but }°>7 must include the 
CCs linked by the path. For (a), the path passes 
inside both the principal CC and the reverse CC, 
which are equal and opposite, so the magnetic 
material adds nothing to gH-ds. But for (b), 
the path links only the principal CC, so the 
magnetic material makes a big contribution to 
£B-ds. Thus, (a) is correct even when }7 is 
redefined to be only the sum of currents flowing 
in actual circuits linked by the path of integra- 
tion (i.e., omitting CCs), but (b) is not true 
under these circumstances if magnetic media are 
present. As }°7 is usually given this definition, 
(a) is the correct form for the circuital relation. 


4. Laplace’s Equation (divB =0) 


The discussion in Sec. V makes it clear that 
lines of B are always continuous, while lines of 
H may be discontinuous at the surface of and 
within a ferromagnetic medium. Consequently, 


it is clear that divB=0, but that divH may not 
be zero. 


5. Boundary Conditions 


Once Laplace’s equation and the circuital rela- 
tion have been established, the boundary condi- 


tions at an interface between media follow by 
standard methods. 


6. Intensity of Magnetization 


It is possible to give a definition of M entirely 
in terms of CCs and to deduce Eq. (1). It is also 
possible to modify the idea of principal CC (re- 
placing it by a mesh of CCs, each unit in the 
mesh being an infinitesimal area but containing 
thousands of atoms) to cover the case when M 
is not constant. However, it seems convenient at 
this point to introduce poles as a convenient 
mathematical device, without physical significance, 
and to continue the theory in terms of poles. 

So far we have assumed we could ‘‘see’’ and 
count lines of flux. This is convenient for de- 
ducing theoretical results, but in practice we 
must rely on instruments for determining B and 
H. Methods of defining B and H, at a point in a 
ferromagnetic body, precisely and without any 
assumptions about visibility of lines, are con- 
sidered in the next section. 


VII. MACROSCOPIC DEFINITIONS OF B AND H 


All three methods of finding B in air require 
movement of the test apparatus. As this is im- 
possible in a solid medium, the only alternative 
is to cut a cavity in the medium. The idea of test 
cavities has fallen into disrepute, but it is surely 
better to accept definitions involving cavities 
than to try to manage without definitions. It 
must be emphasized that these cavities are actual 
holes cut in the material, not theoretical “gaps 
between atoms’ such as have been considered 
earlier. As the cavities contain air, the methods 
of Secs. III and IV for determining B and H in 
air can be used, and the remaining problem is to 
ensure that B or H in the cavity is the same as it 
was in the material. 

The first choice for the cavity is a cylinder with 
its axis parallel to B or H. Removal of the ma- 
terial of the cavity will cause a reverse CC to flow 
on the cylindrical wall.7 To measure B, the effect 


7 An alternative view is to say the material now has two 
boundaries, the outside and the cylindrical wall of the 
cavity, and the principal CC flows on both. It is of little 
importance which view is taken, as long as it is realized 
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of this reverse CC must be made as small as 
possible, which is achieved by making the length 
of the cavity infinitesimal, so that the number of 
amp-turns of the reverse CC is infinitesimal. To 
measure H, the reverse CC must be allowed to 
have its full effect, which is achieved by making 
the diameter of the cavity infinitesimal, so that 
the reverse CC flows in an effectively infinite 
solenoid. 

Since a cavity of (relatively) large diameter 
can be used for B, there seems to be little objec- 
tion to extending it sideways till it reaches the 
edge of the medium. It thus becomes a crevasse, 
and measuring apparatus can (at least in theory) 
be inserted and manipulated. A coil placed in 
the crevasse is in the wrong orientation to experi- 
ence a couple, so that method (B2) cannot be 
used, but methods (B1) and (B3) are feasible, 
and of these (B3) is preferred for reasons already 
given. Hence, the definition of B inside a solid 
medium could be worded as follows: 


The flux density at a point in a magnetic medium is a 
constant times the charge circulated when a small coil is 
moved from a crevasse cut at the point, the walls of the 
crevasse being normal to the flux density and an infinitesi- 
mal distance apart, to a point of zero flux density. The 
constant is the total resistance of the coil circuit divided 
by the effective area of the coil. 


For measuring H, method (H1) is used. The 
cylindrical cavity of infinitesimal radius allows 
no room for manipulating test apparatus. But a 
solenoid whose cross section is a rectangle with 
one side infinitesimal also gives a field intensity 
at the center of (Ni)//, so the cylindrical cavity 
can be replaced by one of rectangular section, 
and this can then be extended along its long 
dimension till it reaches the surface of the me- 
dium and becomes a crevasse. Owing to the 
orientation, the only suitable method of measur- 
ing B in this crevasse is (B2). The definition of H 
inside a solid medium could thus be worded: 


The field intensity at a point in a magnetic medium is a 
constant times the torque on a small current-carrying coil 


that the current on the cylindrical wall is exactly equal 
and opposite to that on the outside, so that for practical 
purposes it fulfils the same role as the reverse CC we have 
considered earlier. 
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placed in a crevasse at the point, the walls and floor of the 
crevasse being parallel to the field intensity and the walls 
being an infinitesimal distance apart. The constant is the 
reciprocal of the product of yo, the current in the coil, and 
the effective area of the coil. 


In both the above definitions, the coil is as- 
sumed to have finite diameter but infinitesimal 
thickness, so that it can fit into the crevasse in 
only one way. If the conditions have to be obeyed 
rigorously the definitions can only be theoretical, 
but in some cases a crevasse wide enough to 
accept a flat spiral coil of fine wire could be 
tolerated. The main criterion would be whether 
or not a crevasse of this width introduced a re- 
luctance appreciable compared to that of the rest 
of the magnetic circuit. 


VIII. CONCLUSION 


It is not suggested that the methods outlined 
in Secs. V—-VII will replace the concept of mag- 
netic poles, or dipoles, for purposes of calculation. 
The advantages claimed for the methods out- 
lined are the following: (a) A much clearer mental 
picture of the meanings of B and H is obtained. 
The way in which a demagnetizing field can be 
produced by atomic CCs is shown, and in so 
doing the artificial nature of the concept of field 
intensity is stressed. (b) The reasons why H is 
used in the circuital theorem; why B is used in 
Laplace’s equation; and why the boundary con- 
ditions use B for normal components and H for 
tangential components, are made clear. (c) 
Rigorous definitions of B and H at a point inside 
a solid magnetic medium are provided. (d) Poles 
are not mentioned in obtaining any of these 
results. They can be introduced when needed, 
but it is clear that they are only a convenient 
mathematical device, not physical entities. 

A further advantage of the methods outlined 
is that the early introduction of the atomic pic- 
ture, with the idea of lines of flux threading suc- 
cessively through lines of atoms, and with local 
distortions between atoms, lays a sound founda- 
tion for explanations of hysteresis and the idea 
of domains. 





A Classroom Antenna Experiment 
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This article describes equipment for studying antenna patterns without the use of cumber- 
some outdoor apparatus. Patterns for numerous arrays, in good agreement with theory, are 
readily obtained. Construction of arrays observed by students is encouraged. 





HIS article describes an apparatus designed 

to conduct antenna pattern experiments in 

college laboratory space. The equipment is not 

cumbersome nor does it use frequencies which 
might interfere with present TV channels. 

To place the receiver in the ‘‘far zone’’ use 
is made of a 12-cm wavelength. Lighthouse tube 
oscillators operating at this wavelength are fre- 
quently used in microwave optics experiments 
and may be purchased from the Central Scientific 
Company. Such equipment makes desirable the 
use of antennas as transmitting units. By re- 
ciprocity (except possibly for some of the “V”’ 
antennas) a receiving antenna pattern is the 
same as that of an identical transmitting 
antenna. 

The transmitting unit shown in Fig. 1 is 
easily constructed at small expense. The turn- 
table, 18 inches in diameter, is cut from ply- 
wood and mounted on the square plywood table 
by means of a bearing (New Departure No. 3203) 
and drawer rollers. Over-all height of the unit is 
about 24 inches. The scale is cut from bristol 
board, marked at one-degree intervals, and is 
glued to the turntable. To measure the angle of 
rotation a movable marker (not shown) is used. 
A two-prong microphone connector, centered on 
the axis of rotation, provides a termination into 
which the antenna can be slipped. A section of 
300-ohm transmission line an odd number of 
quarter wavelengths long is used to couple the 
antenna to the oscillator. 

The coupling loop is a piece of spring wire 
about 23 inches long bent into a “U.”’ The two 
wires of the feed line are soldered to the ends of 
the loop. Proper insertion of the loop into the 
oscillator cavity is necessary to obtain maximum 


* Now at the United States Air Force Institute of Tech- 
nology, Wright-Patterson Air Force Base, Dayton, Ohio. 
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output but no effort is made to match the line 
to the antenna. 

The receiver is an intensity meter consisting 
essentially of a receiving dipole, reflector plate, 
crystal rectifier, and a microammeter with a 
scale range of 0-20 microamperes. This receiver 
is available from the Central Scientific Company 
as a part of their microwave optics equipment. 

To make sure that radiation reaching the re- 
ceiver is due only to the antenna, a framed 
copper wire screen fastened to the base and ex- 
tended to a sufficient height effectively shields 
the receiving meter from transmission line and 
oscillator cavity radiation. 

A Yagi antenna mounted in the transmitter is 
shown in Fig. 1. (The darker background around 
the antenna is merely for contrast.) A pattern 
obtained by students from this antenna with the 
receiver mounted on a wooden tripod about 35 


Fic. 1, Photograph of the transmitter. 
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Fic. 2. Radiation pattern of a Yagi array. The radial 
coordinate is in db, that is 20 logZ/Zo, where J» is the 
maximum current reading of the receiving microammeter, 
I the current reading at arbitrary azimuth. The asymmetry 
may be due to reflection from a wall. 


wavelengths distant is shown in Fig. 2. A minor 
lobe in the 270° position is too small to show in 
the pattern on the scale selected. Reflections from 
objects and walls can affect the meter readings, 
and some care is necessary in the placement of 
apparatus and observers. 

Antennas made of copper wire about 5 inch 
in diameter, pulled to straighten and harden it, 
are mounted in polystyrene rods. The parasitic 
elements are inserted into holes, and the driven 
element is fitted into a slot cut in the mount. 
Details for the Yagi mount are shown in Fig. 3. 

For arrays with a single vertical driven ele- 
ment, 75-ohm coaxial cable is used. The antenna 
is essentially an extension of the central con- 
ductor. Vertical parasites are placed in the de- 


sired position in a thin sheet of polystyrene. This 
sheet slips over the transmitting antenna and has 
a pattern of holes suitable for a variety of arrays.! 

Often the most effective length of the parasitic 
elements is best determined by trial and error. 
Particularly is this true for the single parasitic 
director or reflector associated with a half-wave 
transmitting antenna. A tight fitting sleeve on 
the vertical antenna provides a means for chang- 


Fic. 3. (a) Details of the Yagi mount with folded dipole 
as fed element. (b) End view. The vertical segments are 
spaced to slip into the microphone connector. This figure 
is not drawn to scale. 


ing its length. This is a convenience in deter- 
mining the resonance length for the antenna. 

The author is grateful to George Guthrie, a 
sophomore physics major at Wabash College, for 
constructing the transmitting unit and many of 
the antenna models. 

1 Bailey, TV and other Receiving Antennas (John F. 
Rider Publisher, Inc., New York, 1950). This book gives 


geometrical details of a great variety of arrays and is an 
excellent book to use in conjunction with this experiment. 


Iron rusts from disuse, stagnant water loses its purity and in cold weather becomes frozen; even so 
does inaction sap the vigors of the mind.—LEONARDO DA VINCI. 
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Application of Dirac’s s-Function to Some Problems in Classical Physics 


L. S. KorHari* 
Fitzwilliam House, Cambridge, England 
(Received June 27, 1952) 


The Dirac 6-function is used in solving the problem of bending of a beam and also some 
problems of vibration. It is seen that fairly general results can be obtained without complicated 


calculations. 


HE function 6(x), introduced by Dirac,! has 
been extensively used in quantum me- 
chanics. On account of its queer properties it 
has, so far, found little application in classical 
physics. The purpose of the present paper is to 
show that the 6-function can be used with ad- 
vantage, in solving specially those problems of 
classical mechanics in which forces are concen- 
trated at a point. As examples we shall consider 
(1) the problem of bending of a beam and (II) 
some problems of vibration. It will be noticed 
that fairly general results can be deduced in a 
very simple way with the use of the 5-function.? 
As is well known, the 6-function is defined by 
the following two relations: 


i 5(x)dx=1, 


—o 


(1a) 


5(x)=0 for x0, (1b) 


i.e., 5(x) is zero everywhere except at the origin 
(x=0), where it tends to infinity in such a way 
that Eq. (1a) is satisfied. From this definition one 
can deduce a number of interesting properties of 
5(x), but we state here only the two which we will 
use later in the calculations. They are 


Jf sea@e-aae= so), (2) 


i(x—a)=0 f ¥a, 
(x—a) or xa 3) 


—o for x=a. 


I. We first consider the bending of a beam of 
length / hinged at both ends and acted upon by 
a force P at a point x=a, with (J>a>0). (We 


* On Government of India Scholarship. 

1P. A.M. Dirac, Quantum Mechanics (Oxford University 
Press, London, 1947), p. 58. 

2 Jointly with Dr. B. K. Agarwala, a preliminary note 


(195 _ been published in Physics Quarterly 2, 18 
1951). 


take one end of the beam as the origin.) Then the 
load density can be expressed as 


W(x) = Pé(x—a). (4) 


If Y is Young’s modulus of the material of the 
beam and J is the moment of inertia of the cross 
section, then the displacement of the beam from 


its unperturbed position is given by the differ- 
ential equation 


da‘ 


is — W(x) = —Pi(x—a). (5) 
dx* 


To solve this equation we make a Fourier 
analysis of 6(x—a). We write 


© nex 
é(x—a) => ay - 


n=l 


(6a) 
where 
2 7 nTx 2 nna 
a,=-— f 5(x—a) sin——dx =— sin—-. (6b) 
lL Jo l l l 
Substituting for 6(x—a) from Eqs. (6) in Eq. (5) 
and integrating four times successively, we get— 
for the present boundary conditions—the solution 


2PP 2» 1 ma naxx 7 

y = x i sin ; sin (7) 

It will be noticed that the point a (where the 

force is applied) and the point x (where the dis- 

placement is observed) can be interchanged with- 

out changing the magnitude of the observed 
displacement. 

On account of the factor 1/n*, the above series 
converges very rapidly; and even if we retain the 
first term only, the error involved is not large. 
We therefore take 


2PRh na xx (8) 
y= sin— sin—. 
wYVI l 
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For all values of a (J>a>0) the maximum dis- 
placement of the beam occurs at the center, 
x=1/2. Taking a=1/2, we get 


2PI3 
VI 


Vmax = 
Apart from a small numerical factor, this agrees 
with the result usually deduced. (The small dis- 
crepancy is removed if instead of working with 
Eq. (8) we work with Eq. (7) and note that 
> (1/n*) =2*/90.) With the help of Eq. (7) one 
can also study the bending of a uniformly loaded 
beam. We replace P by Wda in Eq. (8) and inte- 
grate with respect to a from 0 to /, getting 


4Wis | 
y= — sin—, 

wr VI n=1,3,5--- 15 
where W is the load per unit length of the beam. 

II. Asa second example we consider the vibra- 
tions of a struck string. The differential equation 
of a vibrating string is 


oy 


Ox? 


1 d’y 
_——— ’ (9) 
ce ot 


where c?=T/p; p is the mass per unit length of 
the string; and T is the tension applied to it. The 
general solution of the above equation, for the 
case of a string of length / fixed at both ends, is 


© NTX nirct 


niwct 
y= > - a, sin +b, cos). (10) 
n=l 

If the string is excited by striking it at a point 
x =a, and if the blow imparts an “‘average”’ veloc- 
ity ¥ to the string, we take as the initial condi- 
tions 

y=0 forallxat ¢=0, 
dy 
—=y=lv5(x—a) at t=0. 
ot 


These give 
21 
a,=— 
NTC Ho 
2vol nna 
=— sin—, 
nrc 


b,=0. 


nTx 
V0d(x _ a) oe 


KOTHARI 


Substituting these in Eq. (10) we have, for the 
displacement, 


2uol =~ 1 na max — noct 
y=— > —sin—— sin—— sin—, 
TC n=1 17 l l l 


(12a) 


x nmd nmx nC 
y= 2v) > sin—— sin—— cos—— 
l 


n=1 


(12b) 


gives the instantaneous velocity at any point x 
of the string. 

Forced oscillations of a string, when the force 
is concentrated at a point, can also be studied by 
making use of the 6-function.® 

As a further example of the application of 
6-function we consider the vibrations of a canti- 
lever of length /, excited by striking it at any 
point x=a. The vibrations of a cantilever are 
represented by the differential equation 


dty p dy 
corel igen, (13) 
Ox? YI ot? 


where p is mass per unit length of the beam. The 
general solution of this equation is* 


2rnct 2rnct 
y= ¥a(2)| B, cos +C 


n Sin 
n=1 


| (14a) 


where 


TBnX 


¥n(x) = ax| cosh 


TB nX 
— cos 
l 


TBnX TB nX 


—sin 


+b, sinh | (14b) 


2 l 
B, — J Yo(x)Wn(x)dx, 
(14c) 


1 l 
Giieos f (y)odalx) dx, 


TVnl 


and an, bn, and 8, are certain constants, while », 
gives the allowed frequencies of vibration of the 
cantilever. If the cantilever is struck at a point 
x =a, so that the “‘average’’ velocity imparted to 

3P. M. Morse, Vibration and Sound (McGraw-Hill Book 


Company, Inc., New York, 1936), p. 99. 
4 See reference 3, Chapter IV. 
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it is Yo, we have 


Vo 
C, _ —¥rn (a) . 


TVn 


The case when the cantilever is struck at the 
end is the one usually considered, and it is 
simple to show that here 


Vn(a =1) _ (= 4)>*vi. 


The above result is, however, more general, in 
that we can study the vibration of the bar excited 
by striking it at any point. If, for example, the 
cantilever is struck in the middle, a=//2, then 
apart from the first no other odd overtones 
(n=3,5,7, ---) are excited. 

As a final example we consider the vibration 
of a rectangular membrane of dimensions a Xb. 
A vibrating membrane is represented by the 
differential equation 


0°n Om 1 On 


Ox? dy? it 2 af’ 


(15) 
where n(x, y, ¢) is the displacement of the point 
(x, y) at the time ¢ from its original unperturbed 
position. The solution of this equation is of 
the form 


2  mnrx  nry 
n= > sin— =o COS2 TV mnt 


m,n=1 a 


+Cmmn SiN2TVmnt}, (16a) 


where 


4 MTX nTy 
eke f f no(x, 9) sin—— sin——dxdy, (16b) 
ab a b 


DIRAC’S s-FUNCTION 


2 
Cnn = J free y) 
TVmndb 
m 


X sin 


™x  mnTry 

sin——dxdy, (16c) 
a b 
and 


1/T\*( sm? /n\?)3 
maT) 1C)+G) | 

2\ p a b 

m=1,2,3--- 
. (16d) 
n=1,2,3--- 

p is the mass per unit area of the membrane, and 
Ymn represents the allowed frequencies of vibra- 
tion. If the membrane is excited by striking it at 


a point (xo, yo), and if the impact imparts an 
average velocity vo to it then 


Bice — 0, 

200 | MTXo 
sin 

TVmn a 


_ UTYo 
sin : 


mn — 
and hence, 


2u » 1 MnXo | NTYo | Max 
n=— DY —sin sin sin—— 
T m,n=l Vn a b a 


ary | 
x ~~ sin27vmnl. 


(17) 


Any particular results can easily be deduced from 
the general expression (17). 

The above few examples illustrate the use of 
the 6-function in solving certain problems of 
classical mechanics. Its use can be extended to 
more of such problems, where the forces are 
concentrated at a point. Further, as is clear from 
the problem of bending of a beam, even problems 
with uniform distribution of force can be treated, 
indirectly, by this method. 


For example, tests show that voltages as high as 5000 volts can be developed by such common- 
everyday-occurrences as getting up from a plastic-covered chair. Oscillographic records, however, 
indicate that if the resistance between the person and the chair is less than 20 megohms, the voltages 
will not exceed 300 volts, This is less than the minimum sparking voltage in air and is thus not 
hazardous.—Hazards from Static Electricity, U. S. Bureau of Standards Technical Report. 





The Eckart Conditions for a Polyatomic Molecule 


SALVADOR M. FERIGLE AND ALFONS WEBER 
Illinois Institute of Technology, Chicago 16, Illinois 


(Received May 13, 1952) 


The Eckart conditions for a polyatomic molecule are discussed. They are introduced to insure 
the smallness of the interaction between vibration and rotation and define a system of coordi- 
nates in which the vibration of the atoms around their equilibrium positions is described. 


ROM a completely general point of view a 

polyatomic molecule is a mechanical system 
formed by N nuclei and m electrons. To a high 
order of approximation, the motions of the elec- 
trons are independent of the motions of the 
atomic nuclei.! The nuclei, then, move in the 
electronic potential field and form a mechanical 
system with 3N degrees of freedom. This system 
is considered to be isolated, conservative, and 
mechanically stable, i.e., the interactions with 
other molecules are neglected, the forces derive 
from a potential, and there is an equilibrium 
configuration. There are a number of energy 
levels and corresponding states of this mechanical 
system which will be found by the solution of the 
corresponding Schrédinger equation and the 
transition between these states gives rise to a 
spectrum. The theory of molecular spectra deals 
with these two problems: the calculation of the 
energy levels or the frequencies, and the calcula- 
tion of the transition probabilities. We will be 
interested only in the first problem here. The 
first step in this problem is finding the classical 
Hamiltonian of the polyatomic molecule to ob- 
tain later the quantum-mechanical Hamiltonian. 
The potential energy of the molecule depends 
only on the relative positions of the particles and 
is therefore invariant under a rotation or transla- 
tion of the molecule as a whole. In the earlier 
work, the total energy of a molecule was con- 
sidered to be the sum of an energy of vibration 
and an energy of rotation; but refined mathe- 
matical methods as well as empirical data indi- 
cated that this is not the case, and interactions 
between vibration and rotation had to be taken 
into account. In the mathematical treatment, the 
so-called Eckart conditions? are applied in such 


1M. Born and J. R. Oppenheimer, Ann. Physik 84, 457 
(1927). 


?C. Eckart, Phys. Rev. 47, 552 (1935). 


a way that this interaction will be small so that 
a perturbation treatment is possible in the 
quantum-mechanical formulation of the problem. 
The aim of this paper is to give a detailed deriva- 
tion of the Eckart conditions following the ap- 
proach of Sayvetz’ and also to arrive at some of 
the formulas derived by Nielsen, Shaffer, and 
others‘ by a straightforward method and in 
greater detail than is given in the original papers. 

If an atom is displaced in such a way that its 
relative position with respect to the other atoms 
is changed, restoring forces ate called into play 
and thus the displacement from the equilibrium 
configuration will be small. Thus, the molecule 
behaves like a semirigid body which can translate 
and rotate freely while the nuclei move only 
slightly with respect to each other. In some cases, 
however, the relative motion of a part of the mole- 
cule with respect to other parts of the molecule 
may be considerable, as in the internal tor- 
sion of ethane or the inversion of nitrogen in am- 
monia; therefore, in general, we will have semi- 
rigid parts moving with respect to each other 
with their components moving only slightly from 
the equilibrium configuration and the whole 
structure still free to rotate and translate in 
space. 

If the slight motion from the equilibrium con- 
figuration would be frozen, the parts of one 
molecule would be perfectly rigid, the potential 
energy would be identically zero, and the system 
would be left with only M degrees of freedom of 
which 6 (or 5 if the molecule is linear) describe 
the over-all rotation and translation in space, and 
the other M—6 (or M—5) describe the gross 
internal motions (like the torsion of ethane or 
the inversion of ammonia) when they are present. 

3 A. Sayvetz, J. Chem. Phys. 7, 383 (1939). 


4H. H. Nielsen, Revs. Modern Phys. 23, 90 (1951) and 
references given therein. 
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Then the Cartesian coordinates in space-fixed 
coordinate axes are functions of M coordinates 
a,(p=1, 2, ---M), where a, for p=1, 2, --- de- 
scribe the translation and rotation of the mole- 
cule as a whole and a, for p=7, ---M describe 


the gross internal motions. The equations of 
transformation will be 


t= F,(ap), 
-3N. 


(1) 
where k= 1, 2, -- 


The kinetic energy will be simply given by 


OF, OF; 
=2 mx)? = i a 


—Gpkp’ 
kpp’ = Op, OA’ 


= ya A pp’ php’, 
pp’ 


or, in matrix form 
=&,A°a=P,A°—'P, 


where @ is a column matrix formed by the ay, a; 
is its transpose, P is a column matrix formed by 
the momenta P, conjugate to a,, and A®™ is 
the reciprocal of A°. Since the potential energy is 
identically zero, we can write the Hamiltonian as 


2H=2T=P,A°™P. 


If the atoms are now allowed to move slightly 
from their equilibrium configuration, the 3V — M 
degrees of freedom will reappear; let us call 
qi(j=1, 2, ---3N—M) the generalized coordi- 
nates corresponding to this internal motion. 
Since the configuration departs only slightly from 
the rigid configuration described above, the 
Cartesian coordinates are considered to be linear 
functions of the g; which are to be added to Eq. 
(1), and so 


Xe = Fe(ap) + D5 fei(ap)qi- (2) 


Let us now compute the Hamiltonian. The 
velocities are 


OF, Ofij 
=> et 5 eet E Sesdis 


Qp Pi Op 
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and the expression for the kinetic energy becomes 


OF), OF, 
2T= ~ mie = > me— —pity’ 


pr'k =a, Jap: 


pp'ji'k 0a» 


aF x Ofe 
+X m{— ; 


pp! ik Oap IAy’ 


OF, Of; ; 
> asin 


Oay IAy 


OF, 
+20 mifeifes didi +2 “ sts —faglogh 
ji'k 0ap 
Ofc; 
+2 LD mi—fuidiraerdi, 
pivk Ody 


which can be written in the form 


2T= u A pp'Giptp +2 ~ Brides t X Ci3Gi03, (3) 


pp’ 


where 


A pp =A pp’? + oj App igi t a A pp’ isis 
q2 


Byj;= Bp? +L Bois’ 
- 
Cie =D mufiifis, 
k 


OF, OF; 
Asp? =>, m— 


’ 
k OAy OAp: 
P P 


OF, Of; OF, Of; 


Oap IApy 


A pp - =2 m(— 


Oap IAy’ 


Of; Ofx;" 
A pp'i3' = LY m ’ 
k Oap Oap’ 


OF, 
Bp; = i ——firj, 
k Ap 


and 


Of; 
Boi = X mx kje 


Oap 


Equation (3) can be written in matrix form as 


2T =a,Aa+a,Bq+ 4.Cq, (3a) 
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Fic. 1. Nuclear coordinates of the ath atom 
of a polyatomic molecule. 


where & and g are column matrices and the sub- 
script ¢ is used to indicate the transpose. 
The potential energy 


U (xx) = U(Fi +X; feigi) 


may be expanded in a Taylor series about the 
equilibrium configuration; thus 


eu 
2U=> ( ) Saif Gis 
kk’ ii’ \OX,OX 7 zor 


=) bi9i9%, (5) 
ii’ 








where 
bv =L é ) Suifers. (6) 
kk’ \ OX,0X 4/7 2oP 
In matrix form Eq. (5) becomes 
2U =q.bq. (Sa) 


In the equation for the kinetic energy, the first 
group of terms represents the energy of the sys- 
tem made rigid (q¢;=constant) and includes the 
energy of translation and rotation; the third sum 
represents the kinetic energy of vibration (a, 
=constant); the second represents the possible 
interactions between vibration and the other 
types of motion. It is desirable to set the equa- 
tions in such a way that the interaction terms 
will be small enough to be considered as a per- 
turbation; or in other words, that B,; be small. 
From Eq. (4), since the g; are small, we can insure 
the smallness of B,; by setting 


B,;°=90, (7) 
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which was first formulated by Casimir® and, as 











we will see later, gives the so-called Eckart 
conditions. 


NORMAL MOLECULES 


A molecule is called normal when there are no 
gross internal motions and when in the vibra- 
tionless state it behaves like a rigid body, (zero- 
point energy not being considered). If the mole- 
cule is nonlinear, there are six a,(M=6), the six 
variables describing the classical motion of a 
rigid body. These are the three coordinates of a 
fixed point (ordinarily the center of mass) and 
three variables describing the orientation of a 
body-fixed system of axes with respect to the 
space-fixed axes. These are usually taken to be 
the Eulerian angles or the angles describing three 
infinitesimal rotations. 

Let the index k be split into two indices and a, 
where i=x, y, 2, and a=1, 2, ---N, N being the 
number of atoms in the molecule. Thus x;q is the 
ith coordinate of the ath atom in the space-fixed 
system. The following nomenclature will be used 
(see Fig. 1): R: position vector of the origin of 
the moving system O’ with respect to the space- 
fixed origin O, ra: position vector of the ath atom 
with respect to O’, rq’: position vector of the 
ath atom in the equilibrium configuration, 
d,=f.a—f.°: displacement vector of the ath atom, 
€;: unit vector along the ith space-fixed axis, 
e;: unit vector along the ith moving axis. 

Since the displacements are linearly dependent 
on the qj, we can write 


do=X dest (8) 


and 


i= to°+ 2X da;9;; 
7 
so that 


Xia=(R+1a)-e:=(R+10") ei +2, daj-eig;; (9) 


and from Eqs. (2) and (9) we obtain 


Fig= (R+r,°) *e; 
and (10) 


Ff iaj =a; Ci. 
By use of Eqs. (10), the elements of A, B, and 


5H. B. G. Casimir, Rotation of a Rigid Body in Quantum 
Mechanics (J. B. Wolters, The Hague, 1931). 
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C given by Eqs. (4) become 
a(R+r.°) 0(R+1r,°) 
A pp ———— ene 
Oa» Oap’ 


O(R+1.°) dda; 


Oap 


A pp'i= =2 mel 


Oap’ 
O(R+1.’) dda; 
+), 
Oa 0ap 
Ade; ddq;- 
Appi = Le Ma ’ 
@ 0a, Jay 


a(R+1r.°) 


Oay 


Bpi°= DL Ma 


pm lg, 


Ci = Lo madaj* a3. 


Let us obtain now the derivatives with respect 
to a, which appear in Eqs. (11). For a, corre- 
sponding to a translation, 


aa Cy 
Oay 


OF,° 


0a» 


0dq; 
0ap F 


and for a, corresponding to a rotation, 


dap 


where 1, is a unit vector along the axis of rotation 
corresponding to a,. If a, are taken to be in- 
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finitesimal rotations about the moving axes, 


1,=€p. 


(13a) 


Let us write explicitly Eqs. (7) by taking B,;° 
from Eqs. (11). For a, corresponding to a trans- 
lation, by use of Eqs. (12), 


o(R-+r,°) 


By? =L Ma *daj=Cp* D> mada;; 
a a 


Ap 
and Eqs. (7) require that 
Z. Mada; = 0, 


DL medai9i= 


a) 


> m.d.=0, (14) 
which is the first Eckart condition. From the 
definition of d,, Eq. (14) is equivalent to 


L mMlo= X ma’, 


which means that if O’ is taken at the center of 
mass of the molecule at the equilibrium con- 
figuration, then O’ will always be at the center of 
mass of the molecule. 

For a», corresponding to one of the rotational 
coordinates (Euler angle or infinitesimal rota- 
tion), by use of Eqs. (13), 


o 


daj= D0 ma(1pXTa°) - aj 


or, 
= z. Me 
0a 


a | Zz ML a° X da;, 
and Eqs. (7) require that 
Zz, Me¥a’ X dai = 0, 


(15) 
DY ma? X dajQi = 


az 


> mt.°Xd.=0, 


which is the second Eckart condition. The two 
Eckart conditions define the moving system. We 
have to describe the vibrations of the atoms in 
this particular system to insure the smallness of 
the B,;- 

By using Eqs. (12)—(15), the expressions (11) 
can be further simplified. For a, and ay: corre- 
sponding to translation, 


A pp? = XL Ma€p C= Spr DL) Ma, 
a a 
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where 


1 wf p=r, 
0 if p#¥r, 
A pj =A pp ii = 9. 


Spr = 


and 


For a», corresponding to translation and a,: 
corresponding to rotation, 


A pp’? =A pp j =A pp’ =O. 
For both a, and a, corresponding to rotation, 


A pp’? = DL Mal aS pp’ — (La°* Ep) (Ta’* Ep’) |, 
a 


A pp'i= _ Mala dajSpp’ 
—(ra°: fp’) (da; : fp) 


sa (z.° : fp) (d.j ? tp’) |, 
A pp’ ii’ = Zz mal da;: da jb pp’ 


(16) 


— (d,;: fp’) (dajr- a) i} 
since 
(aXb)-(c Xd) =a-bX(c Xd) 
=a-[c(b-d)—d(b-c) ] 
=(a-c)(b-d) —(a-d)(b-c). 
For a, corresponding to translation, 

Bpii =9; 

and for a, corresponding to rotation, 


Boi = = Ma(&p X daj) - da; 


=ey°), Maa; X Gaj’. (17) 


The above results show that—as we should 
expect—the energy of translation is separable as 


2T trans. = R? z. Ma. 


a 


If we leave out this term in Eq. (3a), since it 
is of no interest in connection with band spectra, 


2T =6Aw+oBG+4.CG, (3b) 


where w is used for & and the tilde is also used 
to indicate the transpose. To write the equation 
for the energy in Hamiltonian form, we have to 
define the momenta corresponding to the coordi- 
nates. Writing symbolically @7T/dw for the 
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column matrix whose elements are the deriva- 
tives with respect to the elements of w, 


oT 
P=—=Aw+Bq, (18) 
Ow 


or in components, 


P,= ZL A ppp + Boidi: 
p’ 7 


and 
p=0T/dgG=Bw+ CG, 


or in components, 
Pi= Li Bywj +L CisGy. 
P ¥ 


Equation (19) gives 
g=C—p—C—Ba, 
and by substitution in Eq. (18) we obtain 
P=Aw+BC'p—BC"'Buw 
=(A—BC™B)w+BC-p, (21) 


which shows that the angular momentum con- 
tains a term, which will be designated by L, 
which is due to vibration alone. The effective 
inertial matrix will be designated by J. Thus, 


L=BC-", (22) 


I=A—BC—B. (23) 


In terms of P, L, and , the velocities are 
w=yu(P—L) 
(24) 
q=Cp—CBy(P-L), 
where 
p=I-'=(A—BC"B)—. (25) 
By use of Eq. (24) the Hamiltonian form of 
Eq. (3b) is now written as 
2T = Pw+pq=Pu(P-L) 
+pC—'p—pC“'Bu(P—L) 
=(P—L)u(P—L)+pCp, (26) 
which represents the final form of the Hamil- 


tonian expression for the kinetic energy of a 
polyatomic molecule. 
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The normal coordinates are defined so that 
C5 = 855" 


and 


(27) 
bj = dj0ji". 


If we consider the g; to be the normal coordi- 
nates, Eq. (26) takes the form 


2T=(P—L)y(P—L)+2p, (26a) 


where P, L, I, and yw are given by Eqs. (19), (21), 
(22), (23), and (25) by taking C to be the unit 
matrix. Equation (26a) is given by Wilson and 
Howard® and by Herman and Shaffer’ as a step 
to form the quantum-mechanical Hamiltonian. 


If the atoms in the equilibrium configuration 
6 E. B. Wilson, Jr., and J. B. Howard, J. Chem. Phys. 
4, 260 (1936). 


7R. C. Herman and W. H. Shaffer, J. Chem. Phys. 16, 
453 (1948). 


Reproductions of Prints, Drawings, and Paintings of Interest in the History of Physics 
51. Caricature of Sir John Leslie 


E. C. Watson 
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are in a straight line, only five coordinates are 
needed to define the rigid motion of the linear 
structure, namely, the coordinates of the center 
of mass and two angles. In consequence, there 
are in this case 3N—5 degrees of vibrational 
freedom and therefore 3N—5 coordinates qj. 
Equation (7) gives only five conditions to define 
the system of moving axes and so, these will not 
be uniquely determined. Sayvetz* considers that 
the transformation is made as before, in which 
case one of the angles will enter as an arbitrary 
parameter. Setting this angle equal to zero, com- 
pletes the definition of the moving axes. 

If the molecule possesses gross internal motions 
which cannot be considered as small vibrations, 
as, for example, the ammonia and ethane mole- 
cules referred to above, Eq. (7) gives additional 
conditions which have to be considered sepa- 
rately for each case. 


California Institute of Technology, Pasadena 4, California 


(Received June 2, 1952) 


An excellent caricature of Sir John Leslie, Professor of Natural Philosophy in the University 


of Edinburgh and inventor of the differential thermometer, was made by John Kay shortly 


Str JoHN LESLIE (1766-1832), Professor of 
Natural Philosophy in the University of Edin- 
burgh, is remembered nowadays primarily for 
his invention and application of the differential 
thermometer and his discoveries in the field of 
heat radiation. Physics teachers should think of 
him, however, whenever they repeat for their 
classes the experiment of freezing water with a 
vacuum pump using sulfuric acid to remove 
the water vapor, for it is to him that this well- 
known experiment was due. The other main 
facts of his life, together with a discussion of 
his principal publications, will be found in the 
seventh and eleventh editions of the Encyclo- 
paedia Britannica, in the Dictionary of National 
Biography, and elsewhere, and so will not be 
repeated here. 


before his death. It is here reproduced together with some amusing biographical notes. 


The article in the Dictionary of National 
Biography states that a bust of SiR JOHN was 
executed by SAMUEL JOSEPH, of which a copy by 
Jonn RHINpD is in the National Portrait Gallery, 
Edinburgh, but it makes no mentien of the 
amusing caricature sketched by JOHN Kay 
(1742-1826), the self-taught Scottish etcher, who 
from 1785 until his death produced nearly 900 
plates of the oddities and celebrities of Edin- 
burgh. Almost every notable Scotsman of the 
period, with the exception of Burns, was repre- 
sented, and the etchings, while not entitled to 
high rank as works of art, possess a certain 
quaint originality as well as considerable fidelity 
as likenesses. 

Kay’s work was collected by HuGH PATON 
and published in two volumes under the title 
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DEMONSTRATION 
OF 
CAUSE & EFFECT. 


Fic. 1. Str JonNn LEsLtie (from the caricature 
by John Kay). 


A Series of Original Portraits and Caricature 
Etchings, by the late John Kay, minature painter, 
Edinburgh; with biographical sketches and illustra- 
tive anecdotes (Edinburgh, 1837 and 1838). The 
curious biographical notes that accompany the 
etchings were started by Kay, but were com- 
pleted after his death by JAMES PATERSON, 
author of The History of the County of Ayr, 
aided by Davin LAING, ALEXANDER SMELLIE, 
and other antiquarians. 

Several pages of these biographical notes 
accompany the plate which is here reproduced. 
They read in part as follows: 


“The character of Sir John has been subject 
to some little stricture. All have admired the 
inventive fertility of his genius—his extensive 


knowledge and vigorous mind. As a writer, 
however, his style has been criticised; and he 
has been accused as somewhat illiberal in his 
estimate of kindred merit, while he is represented 
to have been credulous in matters of common 
life, and sceptical in science. ‘His faults,’ says 
his biographer, ‘were far more than compensated 
by his many good qualities—by his constant 
equanimity, his cheerfulness, his simplicity of 
character, almost infantine, his straight-forward- 
ness, his perfect freedom from affectation’ and, 
above all, his unconquerable good nature. He 
was, indeed, one of the most placable of human 
beings; and if, as has been thought, he generally 
had a steady eye, in his wordly course, to his 
own interest, it cannot be denied that he was, 
notwithstanding, a warm and good friend, and 
a relation on whose affectionate assistance a 
firm reliance could ever be placed.’ In ‘this 
character we are disposed to concur. One slight 
blemish, however, has been overlooked—personal 
vanity; for, strange to say, although in the 
eyes of others the worthy knight was very far 
from an Adonis, yet in his own estimation he 
was a perfect model of male beauty. 

“The general appearance of Sir John is well 
represented in the Print which precedes this 
notice. He was short and corpulent—of a florid 
complexion*—and his front teeth projected 
considerably. In later life his corpulence in- 
creased;t he walked with difficulty; and he 
became rather slovenly in his mode of dress—a 
circumstance the more surprising, as his anxiety 
to be thought young and engaging continued 
undiminished.” 


* “What the natural color of his hair may have been we 
cannot say; but in consequence of the use of some tincture 
—Tyrian dye it is said—it generally appeared somewhat 
of a purple hue.” 

t “When unbending his mind from severer labours, the 
knight resorted to Apicius; and to his success in reducing 
to practice the gastronomical propositions of that interest- 
ing writer has been ascribed his somewhat remarkable 
exuberance of abdomen. A legal friend, now, alas! no more, 
once witnessed an amicable contest between Sir John and 
an eminent individual, celebrated for his taste in re 
culinaria. The latter was invincible in the turtle soup and 
cold punch, but the former carried all before him when the 
‘sweets’ were placed on the table. To show how easily the 
victory was won, besides other fruits produced with the 
dessert, the knight, without any effort, devoured nearly a 
couple of pounds of alm@nds and raisins.” 





















Demonstrations to Provide Data for Student Computation in Electricity* 


Tuomas J. BLISARD AND BERNARD A. GREENBAUM 
Newark College of Engineering, Newark, New Jersey 


(Received March 24, 1952) 


Two panel boards were devised for use in obtaining data for student computations in elec- 


tricity, the first for study of Ohm’s law and Kirchhoff’s laws applied to networks, and the second 
for study of line loss in power transmission. The boards are effective teaching aids because of 
their good visibility and the ease with which circuits can be altered. The data are sufficiently 


URING the past two years our staff has 

developed demonstrations to provide data 
for computations in general physics for engineer- 
ing students. These computations parallel the 
content of the classwork and alternate with the 
usual laboratory experiments during the two-year 
physics course. The staff has found that compu- 
tations based on data from experimental investi- 
gations carried on by the group result in 
significant learning that is not obtained through 
the solution of problems prepared by the 
instructor. 

Apparatus available from scientific houses 
was not sufficiently large to enable groups of 
thirty students to observe the procedures, hence 
it was necessary to plan and build the equipment. 
Two of these assemblies will be described. The 
first, a Vari-circuit board, was used to provide 
data for computations involving Ohm’s law, 
equivalent resistance, and Kirchhoff’s laws; the 
second was for computing line loss in power 
transmission. 


VARI-CIRCUIT BOARD 


A section of five-eighth inch plywood 4X3 ft 
was supported in a vertical plane by two short 
slotted feet. Visibility was increased by painting 
the board yellow, and the circuit was made with 
No. 14 black insulated stranded wire stapled to 
the board. (See Fig. 1). Nuway snap connectors! 
were used as terminals for resistors and meters. 
These connectors provided a large degree of 
flexibility and convenience as meters and can 
be connected to any two terminals by merely 


* Based upon a paper read at the twenty-first annual 
meeting of the AAPT, Barnard College, New York, 
January 31, 1952. 


1Nuway Snap Connectors may be purchased from 


Hathaway Manufacturing Company, 4347 West Twenty- 
fourth Street, New York, New York. 





precise to provide for results in agreement with physical laws. 





pressing the connectors on the meter leads onto 
the terminals. Thus, in measuring potential 
drops satisfactory connections can be made 
without disturbing the circuit. Additional compo- 
nents may be quickly inserted in the circuit by 
pulling off leads and snapping the new leads on 
the terminals. The connectors do not become 
loose with long use and continue to provide 
good contact with insignificant resistance. 
Resistor blocks were built by mounting low 
cost radio resistors and snap connector terminals 
on blocks or wood 2X4 xX in. The resistance of 
each resistor was measured with a high precision 
Wheatstone bridge and the values stenciled on 
the edge of the blocks. The blocks were suspended 
on the board by two round-headed screws 
which fitted loosely into holes drilled through the 
board. Thus, it was possible to move or replace 
















































Fic. 1. Vari-circuit board. (Meters shown were 
replaced with 7-in. wide view meters.) 
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_ TABLE I. Sample computation. To find the total re- 
sistance of a circuit by (1) applying Ohm’s law and (2) by 
computing the equivalent resistance. 


R 
= oe Standard 
I (ohms) (ohms) 


8.00 
3.19 
2.99 


(amperes) 


0.216 
0.536 
0.574 
1.33 


V 
R (volts) 


Ri, 
Rig 


—- 
~~ 
NN 


— 
~~ 
Nh 


Rig 

R; (Exp) 
R; (Cal) 
R 


1.26 
0.650 
0.420 
0.200 
1.26 


NNNnu 
pmo et pee bem 
coocoo 


Rs (Cal) 
Rr=1.29+4.35+ 1.67 =7.31 ohms 


resistors easily and rapidly. Meters were mounted 
on larger blocks and suspended on the board in a 
similar manner. This made changes in meter 
positions convenient and the meters were 
available for other experiments. Seven-inch wide 


TABLE II. Sample computation. To compute the total 
resistance from measurements of the potential drop and 
current for the complete circuit. 


E 
(E=Vi 
Vi V2 V3 +V2+Vs;) 
(volts) (volts) (volts) (volts) (volts) 


1.72 5.50 2.10 9.32 9.36 


Rr=E/I=—— =7.44 ohms. Percent difference for the two 
1.26 values of Rr=1.7. 


E 
Measured 


view meters were used. During the operation all 
the students could see the approximate needle 
positions and one or two students near the 
board read the meters precisely for the group. 
The board is set up for a series-parallel circuit 


r 


‘o-0) 
8} 
38 


LOAD (lamp bank) 


Fic. 2. Power transmission board circuit. 
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in Fig. 1. By changing the snap leads at R; and 
R; all the resistances can be connected in series 
within a few seconds. By changing other connec- 
tions various series-parallel combinations may 
be produced, and resistor blocks may be inter- 
changed or replaced with ones of different values. 
Resistor R2 is frequently replaced with a battery; 
thus, a second emf is introduced into the circuit 
for computations based on Kirchhoff’s laws. 
The variable resistor at the source S was used 
to vary the emf of the supply and data were 
taken for several current values. Quite compli- 
cated networks may be set up rapidly by 
snapping additional leads, with or without 
resistances, across various points on the board. 
A sample set of data and computations for 
the series-parallel circuit of Fig. 1 is included 
below. Students computed the resistance of each 
resistor from measurements of potential drops 


TABLE III. A sample student-computation. 


Ps Per- 
(Source) P (Loss) centage 
VsxlI Effi- 


(watts) 


PL 
Vi (Load) Vs 
(Load) VzXI (Source) 
(volts) (amps) (watts) (volts) 


80 
160 
240 
320 
400 


Per- 
Ps—Pt centage 


(watts) Loss’ ciency 


80 
160 
240 
320 
400 


PWNS CPW 
SCKHRAH CHRON DA 
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and currents. The total resistance Rr found by 
calculating the equivalent resistance was in 
close agreement with the value computed by 
using the potential drop across the entire circuit 
and the current at J. The measured emf was 
found to equal the sum of the potential drops 
within reasonable experimental error. (See Tables 


] and II.) 


LINE LOSS IN POWER TRANSMISSION 


The purpose of this demonstration was to 
indicate that the efficiency with which power is 
transmitted increases with increase in voltage. 
This means that the power loss for the same load 
increases as the voltage decreases. 

The apparatus consisted of a Variac, ammeter, 
voltmeter, variable load (1 lamp-blank) and a 
simulated line resistance, a coil of approximately 





DATA FOR STUDENT COMPUTATION 


2 ohms rated at 100 watts. All of these items 
except the Variac were mounted on a panel to 
afford better visibility to a large number of 
students. The circuit is shown in Fig. 2. 

The procedure followed was to maintain the 
output voltage at a constant value of 50 volts 
and to increase gradually the current through 
the load. Currents of 1.6, 3.2, 4.8, 6.4, and 8 
amperes were obtained by carefully adjusting 
the Variac and the lamp load. At these values 
the power delivered to the load was 80, 160, 
240, 320, and 400 watts. Also, at each setting 
of the Variac the potential drop across the 
Variac was taken. 

This procedure was again followed using 100 
volts as the constant potential across the load. 
However, this time the current readings were 
halved, and, therefore, the power delivered to 
the load remained at the same values. The 
corresponding voltage drops across the input 
were recorded. 

From the above data the students calculated 
the power input, (Vs J=watts) at the source, 
and the power delivered (V_xXJ=watts) to the 
load. The difference of these two values was the 
power loss in the transmission lines. The percent 
loss in the line and the efficiency of the transmis- 
sion were calculated also. (See Table III.) 

To permit a better understanding of these 
figures and a clearer picture of the conclusions 
to be reached, two pairs of curves were plotted. 
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Fic. 3. Percent line loss as a function of power at the load. 
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Fic. 4. Percent efficiency as a function of power at the load. 


In Fig. 3, a curve representing percent loss 
versus power at the load was drawn when the 
potential drop across the load was maintained 
at 50 volts, and a similar curve for 100 volts 
was drawn. Similarly, percentage efficiency 
versus power at the load was also drawn and is 
shown in Fig. 4. The graphs clearly indicate 
the advantage of transmitting power at a high 
voltage since the efficiency increases and the 
loss decreases with increased voltage. 

The two demonstrations described have been 
used in freshmen physics classes during the past 
two years. We have found that the students 
became intensely interested in the problems 
proposed by the demonstrations, and they under- 
stood the significance of the computations as 
they applied them to the operating apparatus. 
Principles and laws became meaningful to 
students as they used them in dealing with 
concrete situations. The teachers have found 
the apparatus convenient to operate and useful 
as a teaching aid because of its visibility and 
the degree of precision of the data obtained. 

The authors wish to acknowledge the helpful 
suggestions and advice of Paul O. Hoffmann, 
Chairman of the Technological Division and 
Physics Department at Newark College of 
Engineering. 





The Closed Differential Pulley: Analysis of Its Rotational and Linear Motions* 


LAURENCE ELLSwoRTH Dopp 
University of California, Los Angeles, California 
(Received February 25, 1952) 


A mechanical system of special design, primarily for instructional purposes, is proposed, 
mathematically analyzed, and tested experimentally. It may be regarded as a hybrid form 
between the compound wheel-and-axle and the differential pulley or chain hoist. Although the 
design is simple, the analysis is not too obvious at the start, and is somewhat complicated, ex- 
hibiting coordinated angular and linear motions with constant accelerations due entirely to the 
weight of the movable pulley. While the analysis assumes no friction, the actual heat loss can 
be readily calculated from time and displacement measurements and the equivalent friction 
force determined. Numerous equations and numerical data are included. The system offers 
various challenges to student interest in classical mechanics, both on the lower-division and 
advanced levels. Particular experimental work is suggested and analytical exercises given for 


use by teachers. 


OTATIONAL mechanics, aside from any 

lack of attention on programs, is easily 
underemphasized in the teaching of general 
physics for two reasons. A proportionately large 
amount of the student’s time at the start is 
needed in the mechanics of linear motion, such 
as of falling bodies, including projectiles. Then, 
having learned the basic formulas of translatory 
motion, he finds it simple to apply the same 
algebra to the corresponding problems of rota- 
tion, and to get the right answers without spend- 
ing nearly the time thinking about the physical 
situations as was spent with the problems in 
translation. 


PULLEY SYSTEMS FOR STUDYING ROTATIONAL 
AND LINEAR MOTIONS UNDER NON- 
EQUILIBRIUM CONDITIONS 


Pulley systems are dealt with in general 
physics as cases of equilibrium, the term ‘‘me- 
chanical advantage”’ being tacitly so defined. An 
exception occurs in the Atwood machine, with 
its single pulley and unbalanced weights giving 
acceleration to the system of masses. It appeared 
challenging to study a more complex pulley 
system in which unbalanced forces act to pro- 
duce accelerations. 

The older textbook of physics edited by Duff! 


* Paper presented at the meeting of the American Associ- 
ation of Physics Teachers, Southern California Section, 
held at Pomona College, April 21, 1951. 

1A. W. Duff, A Text-Book of Physics (Blakiston, 1916), 


fourth revised edition, p. 165, problem 18. This problem 
does not appear in the sixth revised edition, 1926, with the 
changed title, Physics for Students of Science and Engineer- 
ing, nor, I believe, in any subsequent editions. 


gave a problem in a nonequilibrium system of 
three coplanar pulleys, as follows: ‘A cord 
passes over two fixed pulleys and through a third 
pulley suspended between them. A mass of 10 
grams is attached to one end of the cord, a mass 
of 5 grams to the other end, and the suspended 
pulley and the attached weights weigh 2 grams. 
The cords being all vertical, what are the ac- 
celerations of the three masses?”’ According to 
the answers given in the book (809, 639, and 724 
cm/sec”, the outside weights moving downward, 
the movable pulley upward), the solution as- 
sumed no friction, and also neglected moments 
of inertia. Thus, the tension in the cord could 
be regarded as everywhere the same, simplifying 
the problem. In such a situation the vertical 
linear acceleration of the movable pulley is the 
average of the accelerations of the two outside 
weights, with algebraic signs reversed. 


THE “CLOSED” DIFFERENTIAL PULLEY 


We now propose to consider a particular 
modification of one of the more complex classical 
pulley systems. This modification is called the 
“‘closed’’ differential pulley, illustrated in Fig. 1. 
It may be regarded as a hybrid between the com- 
pound wheel-and-axle and the chain hoist or dif- 
ferential pulley. Also, it may be thought of either 
(1) asa compound wheel-and-axle with the wheel 
removed, or (2) as a differential pulley in which 
the characteristic free, or ‘“‘open,’’ loop is missing, 
which suggested the term ‘‘closed.’’ In case (1) 
the wheel might be retained but without any 
torque acting on it, the wheel then merely con- 
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CLOSED DIFFERENTIAL PULLEY 


tributing to the total moment of inertia of a 
rotating rigid body. A glance at Fig. 1 shows 
how the device works. 

If an actuating torque from outside were 
applied in the present system, it would be of the 
ordinary differential-pulley type. But here there 
is no torque from outside. The operation is 
therefore due entirely to the weight of the 
movable pulley and any additional mass attached 
to it. 


THE ANALYSIS 


In its mechanical design the apparatus is quite 
simple, rather disarmingly so when one starts 
analyzing. Eventually, however, one finds a 
minimum of twenty-four different magnitudes in 
the equations developed. In addition to these 
twenty-four there are others, related to one or 
more of the twenty-four. Although analytically 
complicated, the mechanical system is adaptable 
to the laboratory in general physics, with the 
analysis abbreviated. For the more technical 
pre-engineering physics, there can be less ab- 
breviation. For a still more detailed study, it 
might be used in advanced classical mechanics. 


ASSUMPTIONS 


While no apparatus is free from friction, it is 
acceptable for analysis to assume this freedom, 
as is often done. Seemingly inconsistent with the 
supposed absence of friction is the further as- 
sumption that no slippage of the one string occurs 
anywhere. This condition is met experimentally 
at the axle by having the two ends of the string 
attached permanently to it, and at the pulley 
by admitting enough friction to prevent slippage. 
For analysis, it can be imagined that the string 
is a frictionless chain, the links of which engage 
cogs on the pulley, as in the chain hoist. 

In the present case we further “‘idealize’’ by 
assuming that the moving elements all lie in the 
same plane. The experimental justification for 
this assumption of coplanarity in getting usable 
“working equations” will appear later. 

The principal assumptions are: (1) coplanarity 
of all the moving elements; (2) no friction, ex- 
cept at rim of pulley; (3) no slippage of string; 
and (4) thickness and weight of string negligible. 


Fixed Axis 


Fic. 1. The closed 
differential pulley. Cer- 
tain vector quantities, “ 
the velocities and ac- 
celerations both angu- 
lar and linear, are in- 
dicated, but not dis- 
placements. 


Movable Pulley 


QUANTITIES APPEARING 


For brevity, the term ‘axle’ will be ‘used 
instead of “‘compound axle.”’ Subscripts . and » 
refer to axle and pulley, respectively. Subscripts 
3 and , mean the same thing, both referring to 


the movable pulley. We have: #, time; 11, 72, 
radii of large and small drums of axle; m, mass 
of pulley and any attached weights; J,, J», mo- 
ments of inertia; g, gravity acceleration; fi, fe, 
tensions in cord at left and right (Fig. 1); 71, 72, 
torques applied by cord at left and right, to 
large and small drums of axle; 7,, net torque on 
axle; 7», torque on pulley; ai, @2, linear accelera- 
tions of left and right cords; a3, linear accelera- 
tion downward of pulley unit; ag, a,, angular 
accelerations; /;, he, hs, linear displacements, 
with respect to the ‘‘ground,” of left cord, right 
cord, and pulley unit, in time ?#, 'time being 
counted from instant of release of system; 2, v2, 
v3, corresponding linear velocities with respect to 
the ground, attained in time ¢; wa, w,, angular 
velocities; $a, ¢p, angular displacements; v,, vp, 
T., Ty, instantaneous frequencies and periods; 
and finally 2,, mp, total revolutions of axle and 
pulley—all in time ¢. 


GROUPING OF THE QUANTITIES 


The algebraic quantities total 24 as a mini- 
mum, as follows: the independent variable ¢; 6 
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independent parameters, 71, 72, Ia, Ip, g, m; 9 de- 
pendent parameters (dependent on the inde- 
pendent parameters but not on the time), f:, fe, 
T1, T2y (Ta), 73) @1, @2, 23, Qa, (ay); and 8 variables 
changing with time, hi, he, hs, v1, v2, V3, wa, 
(wp, Vay Vp» Tay Tp), bay (Gp) May Mp). Adding to the 
24 the quantities in parentheses, we have an 
over-all total of 34 algebraic symbols. 


OUTLINE OF THE DERIVATIONS 


With the beginning student in mind, we give 
a resumé of the mathematical analysis. The par- 
ticular algebraic route followed is probably vari- 
able, and might conceivably be shortened. While 
derivational steps will be interspersed with dis- 
cussion, the results will later be summarized. 
Downward direction is positive. 


RADII AND LINEAR DISPLACEMENTS 


Let us assume that 


r3=(ritre)/2. (1) 


Then since the linear displacements hy, hz of the 
cord at left and right are equal to the lengths of 
arc traversed by any fixed points on the pe- 
ripheries of the two drums, and since these are 
coaxial and fastened together, the corresponding 
angle turned through by the axle is 


oa= hi/ri= —h2/ro, whence hy/he2= —1;/fr2, (2) 


h, being negative because the cord on the right 
moves upward. This is a relation among the 
independent parameters 71, 72, which are inde- 
pendent of the time, and the variables fy, he, 
dependent on the time. 

Now /h; is the length of string that in time ¢ un- 
wraps from the large drum, while hz is the length 
which in the same time wraps up on the small 
drum. Remembering that hz is negative, we have 
hi+he as the total additional length of ‘‘free’’ 
string in time ¢. The free string consists at any in- 
stant of the two vertical segments not in contact 
with either axle or pulley. These remain equal, 
although each increases in length as the pulley 
descends, so that the additional free length is 
divided equally between the cord at the left and 
that at the right (Fig. 2). Therefore, the pulley 
must have descended a distance 


h3=(Ay+h2)/2. (3) 


LAURENCE ELLSWORTH DODD 


This is a geometrical and mechanical fact we 
can accept. Substituting in Eq. (3) from (2) 
gives 


h3=ch,/(2r1) = —ch2/(2re), 


introducing for brevity the symbol c=r,—re. 


(4) 


LINEAR VELOCITIES AND ACCELERATIONS 


It is unnecessary to resort to calculus, any 
more than in problems of freely falling bodies 
where acceleration is assumed constant. How- 
ever, we shall use some calculus. Differentiation 
of Eq. (3) with respect to time gives 


V3= (v1 +02) /2; 
similarly from Eq. (5), 
a3= (a\+a2)/2. 
Differentiating Eqs. (2) and (4), 
01/V2=hy/he= —11/r2, 


(5) 


(6) 


(7) 
(8) 


and 
v;/03= hy/h3= 2ri/c, 


whence by Eqs. (7) and (8), or by Eq. (4), 
V2/03=he/h3= —2r2/c. (9) 


Further, differentiation of Eqs. (7) and (8) 
gives 
Q1/d2=03/02=hi/h2= —r1/r2, (10) 
and 


Q;/a3=0;/03=h1/h3=2r1/c, 
whence from Eqs. (10) and (11) 
Q2/a3=V2/03=h2/h3= —2r2/c. 


(11) 


(12) 


Thus the ratios of the v’s and the h’s are the 
same as those of the a’s, all taken in corre- 
sponding pairs; but whereas the a’s are constant, 
both the v’s and the h’s are dependent on time, 
the v’s varying as the first power of #, the h’s 
as its square, as in ordinary free fall. 


ANGULAR VELOCITIES AND ACCELERATIONS, 
AND TORQUES 


Since the axis of the compound axle remains 
stationary, the v’s being referred to the “ground,” 
we have 

We = v/7r1 = —V2/7r2, 


(13) 


which is the same relation among 4, v2, 71, 72 
















as Eq. (7), while 


Aga /11 = —Q2/Tro, 


(14) 


agreeing similarly with Eq. (10). 

Expressions are now wanted for the corre- 
sponding angular velocity and angular accelera- 
tion of the pulley, which both rotates and has 
translatory motion as a whole. Here we must 
look at the physical picture closely. Referring 
to Fig. 2, we see that point F on the descending 
string at the left, which at any instant during 
the motion has just made contact with the 
pulley, will in any time ¢ thereafter have moved 
from its initial position, measured along the 
string in its instantaneous position at time ?¢, a 
distance h,. However, in this same time the 
pulley center has descended a distance h3, so 
that the length of string that has made contact 
with the pulley rim in this time is hi—h3=s, the 
length of arc through which any point on the 
perimeter must have moved. Therefore, the 
pulley’s angular displacement is 


op=5/r3= (hi—hs)/rs. (15) 
As before, differentiation gives 
Wp = (¥1—03)/73, (16) 
while Eq. (16) yields 
p= (41—43)/75. (17) 
The net torque acting on the axle is 
Ta=Ti1t7T2= firi— fore = aca. (18) 


The component torque 72 on the small drum is 
negative, or clockwise. But the cord tensions, 
fi and fe, are both regarded throughout the 
analysis as positive. The torque on the pulley 
must be 


Tp= (fo—fi)r3= pap. 


The pulley accelerates counterclockwise, and the 
torque has that direction. This means that 
fe>fi, although from Eq. (18), (fer2) < (fini). 


(19) 


AN IMPORTANT PRINCIPLE IN ROTA- 
TIONAL MECHANICS 


Next, connecting the linear acceleration a; to 
the weight mg, we meet a desirably simple situa- 
tion where the working of an important prin- 
ciple in the mechanics of rigid bodies is illus- 
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Fic. 2. Spread diagram as an aid in analyzing the per- 
formance of the closed differential pulley. Diagram shows 
the relation 4:=h,+s, where s/r3=¢, the angular dis- 
placement of pulley. Not to exact scale throughout: 
actually, |h3| <|he|. 


trated. The situation is simple because the 
forces, only three in number, are not only co- 
planar but also parallel. The principle states 
that ‘‘all the forces applied to a rigid body may 
be replaced by a single force and a single couple.’” 

Suppose that the beginning student reasons 
wrongly in the following manner: Of the total 
linear force, fi+f2, applied through the string 
and acting upward on the pulley, Fig. 1, the 
amount (f2—f1) is ‘used up” in producing, 
according to Eq. (19), angular acceleration of the 
pulley. Therefore, there is left, as the total ex- 
ternally applied force opposing the weight mg, the 
amount (f1+ fe) — (fe—f1) =2f1, so that the equa- 
tion for the linear acceleration of the pulley unit 
would in that event be mg—2f1=ma;. But if the 
student proceeds on this wrong basis, he will 
not find the resulting equations agree with the 
principle of the conservation of energy, mghs 
= $1 .wa?+ 31,w,?+3mv;?. But applying the prin- 
ciple above to get the resultant force on the 
pulley, we subtract from mg the sum (f:+f2), 
while the resultant couple is 7, as in Eq. (19). 


2A. G. Webster, The Dynamics of Particles and of Rigid, 
Elastic and Fluid Bodies (B. G. Teubner, Leipzig, 1912), 
second edition, p. 209. 
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The results will now agree with the energy 
equation. No part of (fi+f2) is ‘used up” in 
producing angular acceleration. This important 
relation as it affects the pulley is then 


mg — (fit fe) =mas. (20) 


THE DEPENDENT PARAMETERS 


By choosing from the foregoing, we seek ex- 
pressions for the dependent parameters—forces, 
torques, accelerations—which are independent 
of time, in terms entirely of the independent 
parameters, 71, 72, Ia, Ip, g, and m. 

Eliminating a, from Eqs. (14) and (18) leaves, 
@1=71(firi— fore)/Ta (21); and a, from Eqs. (17) 
and (19), ads=a:—r37(fe—f1)/Ip (22). Eliminate 
a; from Eqs. (6) and (22), @2=a1—2r3?(fe—fi)/Ip 
(23); and a2 from Eqs. (14) and (23), a1=rib(fe 
—f,)/(2Ip) (24) where b=7,+ 72. Finally, eliminate 
a; from Eqs. (21) and (24), fe/fi=(Jab+2Ipn)/ 
(Iab+2Ipr2) =A/B (25), where A and B are con- 
venient substitution symbols, for functions of the 
independent parameters only. Accelerations are 
constant by Eqs. (21), (30), and (31), a conclu- 
sion reachable also by a physical reductio-ad- 
absurdum argument. 

Equation (25) is of a symmetrical type, the 
only difference between numerator and denomi- 
nator being in the second term, where 7; appears 
in the numerator and 72 in the denominator. 
Since r1>7r2, we have f2>/f1; that is, tension in 
the cord on the right is greater than that on the 
left. This indicates that the analysis is correct 
thus far, since the pulley accelerates counter- 
clockwise. 

Eliminate a; from Eqs. (6) and (20), a2=2g—a, 
—2(fitfe)/m (26); and a2 from Eqs. (14) and 
(26), a@1=2r;[mg—(fitfe)]/(mc) (27); then a; 
from Eqs. (24) and (27), fe=[(mbc—4I,)fi 
+4I mg |/(mbc+4I,) (28). 

For brevity write I for J,+J,, then let Q=4I 
+mc (29); then eliminate f2 from Eqs. (25) and 
(28), f1=2Bmg/(bQ) (30); which is the desired 
expression for f;, in terms entirely of the inde- 
pendent parameters. From Eqs. (30) and (25), 
f2=2Amg/(bQ) (31), which is the desired ex- 
pression for f2. 

For further brevity let S=(ri—re)mg/Q 
=cmg/Q (32). Then by Eq. (30), f:=2AS/(bc) 
(33); while by Eq. (31), fe=2BS/(bc) (34), 
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where A = (JI,b+2J,r1), and B=(I,b+2Ipr2), as 
in Eq. (25). From Eqs. (30), (31), and (32), 
(fe—fr) =41pS/b = 2IpS/rz (35). Similarly, fit fe 
=4IS/c (36). Substitute Eq. (35) in Eq. (24), 
a,=2r,S (37). From Eqs. (23), (35), and (37), 
a2= —2r2S (38); and by Eqs. (6), (37), and (38), 
a3=cS (39). Equations (37), (38), and (39) are 
the desired expressions for the three linear 
accelerations. 

From Eggs. (18), (33), and (34), r2=2I4S (40); 
and by Eggs. (19) and (35), r,=2I,S (41). The 
net torques on the axle and pulley are counter- 
clockwise. Incidentally, from Eqs. (18) and (19), 
Ta/Tp=Tatta/(Ipap) (42); while from Eqs. (40), 
(41), and (42), ta/tp=Ia/Ip (43). By Eqs. (42) 
and (43), aa=a,, which means that axle and 
pulley always agree in phase. To check this re- 
sult, by Eqs. (14) and (37), aa=2S (44); and 
from Eqs. (17), (37), and (39), a,=2S=as=a 
(45). This result is obtainable directly from Eqs. 
(11), (12), and (17). 

Now 71= f171, whence by Eq. (30), 71: =71Amg/ 
(6bQ) =27r,A S/(bc) (46). Also r2= fore, whence by 
Eq. (31), 72= —2r2Bmg/(bQ) = —2r2BS/(bc) (47). 

To check Eq. (40) by Eqs. (46) and (47), 
since ta=7T1+72, T2 being negative, we get 7a 
=2I,S, the same as Eq. (40). Further, by Eqs. 
(46) and (47), we have 7;/72= —Ar;/(Bre) (48), 
which shows that since A>B, and r1>re, |71| 
>|72|, although, as noted above, f1< fe. 

Also, since by Eq. (19), tp=(fe—fi)rs, by 
Eq. (35), tp=2I,S (49). Incidentally, by Eqs. 
(46) and (49), 71/tp=71A/(Ipbc) (50); and by 
Eqs. (47) and (49), 72/tp= —7r2B/(Ipbc) (51). 

The desired expressions have thus been ob- 
tained for all the dependent parameters in terms 
of the independent ones. 


DERIVATIONS FOR THE VARIABLES 
DEPENDENT ON TIME 


We next want expressions, in terms of inde- 
pendent parameters, for the variables dependent 
on time. These include the velocities and dis- 
placements, both linear and angular. 


(A) The Pulley Displacement h; as 
Independent Variable 


The distance hz; through which the pulley 
has descended is a function of time, 43=cSf?/2 
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(68). Because this distance is readily measurable, 
expressions in terms of it are desirable. Accelera- 
tions are constant. 

032 = 2ayhy (52); v2? = 2Zaehe (53); v3? = 2ashs (54). 
From Eqs. (2) and (3), 41=2rih3/c, (55); and 
similarly, he= —2reh3/c (56). By Eqs. (25), (36), 
and (39), 1=2riw/(2Sh3/c) (57); by Eqs. (26), 
(37), and (40), ve= —2re,/(2Sh3/c) (58); and by 
Eqs. (3), (27), and (38), v3=+/(2cSh3) (59). To 
check, note that Eq. (58) could have been got 
directly by Eqs. (7) and (57), and likewise Eq. 
(59) either by Eqs. (8) and (57), or by Eggs. 
(9) and (58). Equations (55) to (59) are the 
desired expressions for linear displacements and 
velocities in terms of h3, the pulley displacement, 
a function of time. 

Now, by Eqs. (13) and (57), wa=v1/ri 
= 2(2Sh3/c)* (60); while by Eqs. (16), (57), and 
(59), wp=(¥1—v3)/r3=2(2Sh3/c)t=wa=w (61). 
Since t=0 at the instant of the system’s re- 
lease, da=a?/(2aa), whence by Eqs. (44) and 
(60), ¢4=2h3/c (62); and similarly by Eqs. (45) 
and (61), ¢,=2h;/c=¢.=¢ (63). Equations (60) 
to (63) give in terms of h3 angular velocities and 
angular displacements of axle and pulley, and 
they are the same for both. 

In view of these results, we have for the com- 
mon frequency of rotation attained in time f, 
v=w/(2r) =(2Sh3/c)*/m=1/T (64), where T is 
their common period. Also, the total revolu- 
tions are the same for both pulley and axle, 
n = o/(2mr) = h3/(mc) = hy/(24r1) = — ho/ (2282) 
(65). 


(B) The Time ¢ as Independent Variable 


Since in general h=ai?/2, we have by Eggs. 
(37), (38), and (39), respectively, h1=1r,St? (66); 
he= —r2S? (67); hs=cS/2 (68). Also similarly, 
V1=2r,St (69); ve=—2reSt (70); vg=cSt (71). 
For the corresponding angular quantities w and 
¢, since w=at and ¢=al?/2, by Eqs. (44) and 
(45), wa=2St=wpy=w (72). Similarly ¢,.=Sf 
=¢,=¢ (73). Equations (66) to (73) are the 
desired expressions for linear and angular ve- 
locities and displacements in terms of time. 
Frequency v=w/(2r)=St/r=1/T, (74); 
total revolutions n = ¢/(27) = St?/(27) (75). 


and 


SUMMARY OF EQUATIONS FOR THE 
DEPENDENT PARAMETERS 


For convenient reference, in addition to sec- 
tions (A) and (B) just above, we assemble 
“working equations’ from the expressions pre- 
ceding those, remembering the mechanical con- 
dition that the pulley diameter is the average of 
the two axle-diameters, r3=(ri+-7r2)/2=6/2 (1) 
bis. Wherever the factor b=(ri1+-72) appears, it 
is easy to introduce r; into the equation. As to 
the factor c=(r1—r2), we shall comment later 
under “Equivalent Axle.’’ Conciseness, as has 
been seen, is served by the abbreviating func- 
tions Q and S. These are parametric in character, 
since they are functions entirely of the inde- 
pendent parameters, S being a function of Q. 
They also will be discussed presently. We have, 
Q=4I+mce?, (29); S=cmg/Q, (32). 

Cord tensions: f1:=2AS/(bc) (33); fe=2BS/ 
(bc) (34); where A=J,b+2I,r1, and B=I,b 
+2Ipr2, as in Eq. (25). Linear accelerations: 
descending cord, @;=2r,S (37); ascending cord, 
ad2= —2re2S (38); pulley unit, ag=cS (39). The 
torques: on large drum of axle, 7; =27,Amg/(bQ) 
= 2r,A.S/(bc) (46) ;0n small drum, rz = — 27r2Bmg/ 
(bQ) = —2reBS/(bc) (47); net torque on axle, 
Ta=T1+7T2=fiti— for2=2I,S (40); net torque on 
pulley, tp=(fe—fi)rs=2I,S (41). Angular ac- 
celeration common to both axle and pulley, 
a= 2S (45). 


QUANTITIES OTHER THAN A; OR ¢t AS 
INDEPENDENT VARIABLES 


Obviously with time considered as a parameter 
any of the dependent parameters, themselves in- 
dependent of time, can be treated as independent 
variables, some other dependent parameter being 
regarded as dependent variable. By plotting 
these two the locus of the equation can be 
studied. Two examples follow: 


(a) Phase angle $ as independent variable. From 
Eq. (73), 2=¢/S, whence by Eq. (66), 4i=1719; 
by Eq. (67), ho=—red; by Eq. (68), hs=cd/2; 
by Eq. (69), 1=2riv/(S¢); by Eq. (70), v 
= —2re/(Sd); by Eq. (71), vs=e/(Sd); by 
Eq. (72), o=2/(S¢); by Eq. (74), »p=V/(S)/7; 
and by Eq. (75), »=¢/(27). This group of equa- 
tions is collectively designated (76). 

(b) Total rotations n as independent variable. 
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Since ¢=22n, we have by Eq. (73), 2=2n/S, 
so that substitution again in Eqs. (66) to (75) 
yields the obvious relations: hiy=2mrin; he 
=—2nrn; h3=men; v3=2r;(27Sn)!; ve= —2re 
<(2rSn)!; v3=c(2rSn)!; w=2(2rSn)'; and vp 
= (2Sn)!/x. These forms together are designated 
(77). 
CHECKING THE ANALYTICAL RESULTS 


Accuracy of the derivations was checked by 
dimensional inspection, interchecking among the 
various results, substitution in the energy equa- 
tion, and experimentally as described below. 


COMPOUND AXLE 
(Center section) 


Scale: % 


All dimensions 
given in mm. 


Connecting 
Strings 


fl 
| 
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OV: 
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Fic. 3. Vertical axial section of fixed compound axle 
and movable pulley, to scale. The string is shown un- 
wound from larger drum of axle, so that pulley is at its 
lowest level. Although pulley is shown here with its plane 
perpendicular to plane of paper, actually in this lowest 
position its plane as turned about a vertical axis through 
some 30° from its initial approximate coplanarity with 

‘plane of axle,”” here perpendicular to plane of paper. 
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APPARATUS 


The compound axle was lathed from a single 
piece of solid aluminum, and similarly the pulley. 
The vertical section, Fig. 3, is to scale. The 
string, over 9 ft long, had its ends fixed to the 
drums by drilling a hole from the periphery of 
each drum next the flange obliquely through the 
axle, to emerge through the outside flange-plane 
of the larger drum, at the right (Fig. 3), so that 
the two emergence openings are diametrically 
opposite on this face. Each opening is approached 
by the wound-up cord in a gradually deepening 
gouge in the flange, the fishline string keeping its 
helical path until it disappears into the flange. 
The end segments of the string were threaded 
through the holes from the drum peripheries and 
knotted. In effect, the cord was attached to each 
drum on its periphery. Each successive turn has 
close contact with the preceding and the follow- 
ing turn. Strictly, the 1-mm diameter of the 
cord should be added to diameters of drums and 
pulley, for calculating values in checking the 
equations, but the error is small. In fact it 
affects 7; +72 but not 73—7r2, so that neither of 
the functions Q and S is affected. 

Figure 4 gives comparison photographs of the 
apparatus, made at four different levels of the 
pulley during its descent. The agreement in 
phase is evident. 


NUMERICAL VALUES OF PARAMETERS 


Referring to Fig. 3, r;=2.65, r2=2.25 cm; 
b=4.9, c=0.4 cm; m=65 gm; g=980 cm/sec’. 
Calculated: J,~=705, I,=69 g-cm?, treating 
the axle as an assemblage of solid cylinders. 
From these values, J,=m(r1—r2)?=10.8, while 
Q=3107, S=8.2008. These precise values for Q, 
S assume similar precision in measuring the in- 
dependent parameters. Further abbreviating 
functions are: A=IJ,b+2I/,r2=3765.0, B=I,b 
+2I,r1= 3820.2 absolute cgsu. 

Minimum vertical distance center to center of 
pulley and axle, 9.7 cm; maximum distance after 
16 turns, 9.7+20.1=29.8 cm. Because of the 
2-mm thickness of the middle flange, Fig. 3, the 
plane of the pulley at its upper limiting position 
lacks by calculation and measurement about 3° 
of being parallel with the ‘“‘plane of the axle,” 
or plane of the analysis. At the lower limiting 












position reached after the maximum 16 turns, 
the pulley plane is inclined about 33°, since it 
veers steadily during the descent. However, even 
here, as calculation shows, the free lengths of the 
string do not depart from the vertical by as 
much as 1°, which accounts for the pulley’s not 
being perceptibly out of phase with the axle 
even after reaching the lowest level, after 16 
turns of axle and pulley, Fig. 4. In Fig. 3, the 
vertical section of the pulley perpendicular to its 
plane is shown as though it remained parallel 
with the axle throughout the descent. 

Applying some of the derived equations, which 
assume no friction, we have by Eqs. (32) and 
(39), a3/g=0.003347, or the linear acceleration 
of the pulley is only about 1/300 that of free 
fall. Also by Eq. (20), the sum of the string ten- 
sions fi+f2 is 99.665 percent of the pulley’s 
weight. This value is consistent with the pulley’s 
low linear acceleration. Whereas its weight is 
65.000 g, the total upward pull of the string on 
the pulley is 64.782 g. By Eq. (25), ratio fe/fi 
=B/A=1.0147, so that the tension fe is only 
about 1.5 percent greater than f;, while by Eq. 
(48), the ratio of torques on the axle is 7/72 
= 1.1607, 7; being 16 percent greater than 7, in 
spite of the force f; on that side, the left, being 
less than fe on the right. Tension f,;=32.155, 
while fe= 32.627 g. (For further data see section 
“Numerical Table” and Table I near end of 
paper.) 

EXPERIMENTAL CHECKING 


The calculated 3; for 15 turns was 18.8 cm, 
Eq. (65), which checked with direct measure- 
ment. With a stop watch, not a precise method 
considering shortness of the measured interval, 
the pulley’s descent through 15 turns was timed. 
However, several successive measurements agreed 
well and the measured time was taken as 4.5 sec. 
The calculated time, Eq. (68), was 3.38 sec. An 
error in this direction was expected, because 
friction slows the descent. 

Here two questions may be raised: (1) If the 
apparatus had no friction how much would the 
total moment of inertia J=I,+J, have to be 
increased to give the observed time of descent of 
4.5 sec, the system remaining otherwise the same? 
By Eq. (68), the new value of S would be 2h3/ (ct?) 
=2 X18.849/(0.4 X20.25) =4.6541. Then by Eq. 
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(a) (b) (c) (d) 


Fic. 4. Photographs of apparatus made at four different 
levels of pulley during its descent under gravity. Radial 
markers show that pulley and compound axle remain in 
phase. Also, visibility of the light-colored string is in- 
creased by intermittent inking. In (d), with pulley at 
lowest level, inclination of its plane from the initial ori- 
entation in (a) is apparent. This obliquity causes in the 

hotograph an apparently much greater departure from 
Coleman parallelism of the two radial markers than 
actually exists. 


(32), Q=cmg/S=0.4 X65 X980/4.6541 = 5474.7 
=4]I+mc*, by Eq. (29); therefore, the required 
new value would be J=1366. The actual value 
calculated from measurements was 774. There 
would have to be an increase of 76.5 percent in 
the combined moment of inertia to slow the 
descent as much as was actually caused by 
friction. 

(2) What is the equivalent force of friction? 
That is, assuming no friction what external 
force F applied vertically upward on the pulley 
would slow the descent down from the ideal 
calculated value of 3.38 sec to the observed 
time, 4.50 sec? Since distribution of friction in 
the system is unknown, f; and fe cannot be 
calculated and we resort to energy relations. In 
the descent of the pulley through h;=18.8 cm, 
for 15 turns from rest, loss of PE is mgh;=65 
980 X 18.8 = 1,200,000 ergs. The actual KE of 
the whole system is Iw?/2+mv,?/2. By Eq. (16), 
w=(v;—v3)/r3, and since v1=w71, we have 23 
=we/2. Then KE=([+mc?/4)w*/2, where the 
second term is small. To find w=at=a;t/r;, we 
have from Eqs. (1) and (10), a1=271a3/c, where 
a3=2h;3/t?, t being the observed time. There 
results @3=1.85, w=2a3t/c=2 X1.85 X4.5/0.4 
=43.9. Then KE=(774+65 X.16/4) X1927/2 
= 748,000 ergs, which is 62.3 percent of the loss 
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of PE. The difference 452,000 ergs=0.011 cal, 
or 37.7 percent, is friction loss. The equivalent 
friction force is then F=452,000/18.8 = 24,000 
dynes = 24.6 g, 38 percent of the pulley weight. 


THE FUNCTIONS Q AND S 


Let us attribute physical meanings to the 
abbreviating functions: Q=4(J.+],)+m(ri—re)? 
(29); and S=(r1—7r2)mg/Q (32), the latter only 
containing the gravity factor. Dimensionally Q is 
moment of inertia, while S, being torque divided 
by moment of inertia, is angular acceleration. 
Physical interpretation for S is left to the reader. 
As we shall see, Q is a measure of the total in- 
ertial reaction of the system, while S is affected 
by the actuating force mg. 

Although S is proportional to g, it is not pro- 
portional to m, or mg, because m appears also 
in one term of function Q. The tensions in the 
string, fi, f2, and all accelerations, velocities, and 
displacements, both linear and angular, are pro- 
portional to S. Thus all these are proportional to 
g but not to m, although their values are affected 
by it. 

This absence of direct proportionality of cord 
tensions and accelerations to the driving force 
is met in the analogous simpler system of two 
masses connected by a string passing over a 
pulley, the larger mass resting on a smooth table 
top. The whole system accelerates according to 
the second law, a=mg/(M-+m). If M is very 
large compared with m we have a good compari- 
son with the present case, where m(r1—7r2)? in 
function Q is only about 3 percent of Q. 

Suppose further a=a; (=3.3 cm/sec’), or 
mg/(M-+m) =(ri—re)S, by Eq. (39). Then from 
Eqs. (29) and (32), M=Q/(r1—r2)?—m=3107/ 
0.16—65=19 354 g, or about 42.5 lb being 
accelerated by a weight of 65 g. That is, of 
the whole inertial reaction, the larger part due 
to the relatively large moments of inertia of axle 
and pulley is equivalent in the analogous system 
to M=42.5 lb on the table top, while a weight 
of 65 g accelerates both masses. Ratio m/(M 
+m) equals ratio mc?/Q=a;/g in the present 
system. 

To discuss Q from another point of view we 


note that Q/4=(I.+I,)+m[(ri—r2)/2P=I+1., 


where I,=m[(ri—7r2)/2 P=mr2. Then Q/4 is 
the sum, total moment of inertia J plus that of 
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an equivalent particle of mass m at distance r, 
=(r1—r2)/2 from an axis. Suppose this particle 
has an angular acceleration a equal to that of the 
axle and pulley in common. Then a=a,/r.=25S, 
by Eq. (45), so that a,.=(r1—72)S=cS, where a, 
is the tangential acceleration of the equivalent 
particle. Then by Eq. (39), a.=a3, which also is 
obviously the tangential acceleration of any 
point on the axle located at distance r. from 
rotation axis. Reversing this reasoning, from 
Eq. (39) a3/g=ml[(r—12)/2 P/(Q/4) =1-/(I+ 1). 
Thus Q/4=IJ+I,., the relation started with 
above. 

Function Q is a measure of the total inertial 
reaction in the system, in terms exclusively of three 
moments of inertia—those of axle and pulley 
together with that of the equivalent particle. 


EQUIVALENT AXLE 


Referring again to the equivalent particle, 
since @,.=4a; is the tangential acceleration of any 
point located on the axle at distance r.=(r; 
—r»)/2 from its axis, imagine next a third drum 
attached coaxially to the compound axle, with 
radius r, but negligible mass, so that it has no 
effect on the motions in the system. Let a weight- 
less cord be wrapped around it, with one end 
fastened to the drum, and a small weight, 
negligible but just enough to keep the string 
taut, attached to the free end. If initially this 
small weight is hung at the same level as the 
pulley axis, but not touching it, then after re- 
lease of the system the weight, accelerating at 
rate @.=4d3, will remain on the same level as the 
pulley. For the present apparatus, the diameter 
of this imaginary drum is c=7,—7r2=4.0 mm, 
which contrasts with diameters of the actual 
drums, 26.5 and 22.5 mm. 

Now suppose this “equivalent axle’ is an 
actual drum replacing the whole system. Let 
the mass hanging from the string be m, the same 
as that of the replaced pulley. We ask, what mo- 
ment of inertia 7, would be necessary for a,=4d3. 
We have 7t.=mg-r.=mg(ri—r2)/2=Ia=I,-2S 
by Eq. (45) so that J.=mgr./(2S)=Q/4, the 
total inertial reaction of the replaced system, 
in terms exclusively of moment of inertia. Thus, 
wherever we meet the factor c in the equations, 
we can be reminded of this equivalent axle, 
the torque on which is r,=QS/2=cmg/2. 
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TWO ROLES OF THE PULLEY MASS 


In the closed differential pulley when operat- 
ing, the mass m of the pulley unit plays two 
separate roles: (1) in the force factor mg its role 
is merely to supply a material object the weight 
of which accelerates the system; (2) in the m- 
term, m(ri—re)?, present in the function Q— 
this m-term being 4 times the moment of inertia 
of the equivalent particle—its role is to con- 
tribute to the total inertial reaction of the system. 

A similar situation exists in the analogous case 
where a=mg/(M-+m). Here the role of m in the 
numerator, where it determines the applied force, 
is different from that in the denominator, where 
it presents inertial reaction. Thus, in general, 
any mass has the two roles mentioned, whether 
or not it is part of a complicated system. The 
present apparatus could be laid on a smooth 
table top, and a force F, applied to the pulley 
along the line of centers of pulley and axle that 
would be entirely independent of its mass, thus 
eliminating the first role of m while retaining the 
second. This means substituting F, for mg in Eq. 
(20). If from this point, further derivation pro- 
ceeds on this basis, the same expressions are 


obtained throughout, except that F, now appears 
in place of mg. The inertial term m(r,— 12)? for the 
linear motion of the pulley remains. 


ADAPTATION TO STUDENT LABORATORY 


By suitable selection from the working equa- 
tions, the closed differential pulley is available for 
student laboratory. Procedure can be simplified 
for lower division classes, or made more elaborate 
for an advanced course in mechanics. For the 
latter, derivations can be required. 

The vertical distance of descent of pulley h; 
is measurable quite precisely. If a stop watch is 
used a longer path of descent is desirable. This 
can be supplied with apparatus of the present di- 
mensions by using finer fishline or strong thread, 
giving a correspondingly larger number of turns 
on the axle, or by increasing somewhat the widths 
of the drums, while retaining the 1-mm cord. 

To prevent jerking when the pulley reaches 
its lowest position, where the cord starts abruptly 
to wind again on the large drum, a simple brake 
can be quickly applied to the axle just before the 
limiting instant, and also a small platform can 
be placed in the pulley path. 
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By weights one can vary the mass of the 
pulley unit to control the descent time over a 
chosen distance. Thus Eq. (68), where S depends 
on m, can be the working formula. If h; is to be 
treated as a variable parameter, the relative 
values of the independent parameters can be 
changed from the particular ones here, in order 
to give a wider spread in values of t. 

The experiment can include measuring para- 
metric quantities 71, 72, Ia, J», and m. Moment of 
inertia J, can be calculated, as was done here, on 
the basis of the compound axle’s being an as- 
semblage of coaxial solid cylinders, while if the 
pulley is a solid cylinder, J, is calculated simi- 
larly. Thus, the student himself finds the neces- 
sary numerical parameters. For advanced labora- 
tory the experiment has rather extensive pos- 
sibilities. 

SUGGESTED EXERCISES 


The descending pulley offers practice with 
vectors, if attention is fixed on a point P at the 
pulley rim. Particularly choose the point which 
at t=¢=0 was at the right end of a horizontal 
diameter. This point revolves counterclockwise 
according to Fig. 1, about the center of the pulley 
as it descends, its position on the circle at any 
instant being given by the phase angle ¢. Let the 
origin of Cartesian coordinates be at the position 
of the pulley center when at the highest level. 
The origin is thus stationary with respect to the 
ground. If the vector is drawn full scale its tail 
can be at the origin while the head remains with 
the moving point P. The motion of point P, 
particularly with respect to the ground, is some- 
what complicated. 

Exercise 1. In time t what is the instantaneous 
displacement of point P from the origin? The dis- 
placement vector changes length while oscillating 
through the y axis. The x component is 73 cos¢, 
while y=h3;—r3 sing, whence by Eq. (63), y 
=(c@—b sing)/2, downward direction positive. 
From these the resulting vector has length R 
= (b?—2bco sin¢d+c*¢’)!/2. The direction angle y 
with the y-axis is given by tany=x/y. 

Exercise 2. What is the instantaneous velocity 
of point P with respect to the ground, that is with 
respect to the origin? It is the vector sum of the 
tangential velocity vr of P with respect to the 
pulley center and the velocity v; of the pulley 
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TaBLE I. Calculated values of variables dependent 
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center with respect to ground. Basically, this is 
the problem of a boat on a river that moves with 
respect to the bank. Here, however, both compo- 
nent vectors are variable. By Eqs. (72) and (73), 
vr = 2r3(S¢)!=b(S¢)', and by Eq. (71), v3=c(S¢)?. 

Exercise 3. What is the instantaneous linear ac- 
celeration of point P with respect to the ground? 
Here three vectors are added: (1) radial accelera- 
tion, toward the pulley center, a, =7;w? = 2bS¢, by 
Eqs. (72) and (73); (2) tangential acceleration, 
ar=r3;a=bS by Eq. (45); and (3) acceleration of 
the pulley center itself, aa=cS by Eq. (39). This 
resultant vector presents a challenge, requiring 
both length and direction as functions of ¢. The 
length of the first of the three components de- 
pends on time, that is on ¢, while the lengths of 
the others do not. However, directions of both the 
first two depend on ¢, while that of the third 
does not, since that component is constant. 

Exercise 4. The actual path traversed by point 
P with respect to the ground is followed by the 
head of the displacement vector, if it is drawn 
full scale and treated as above. What is the 
equation of this path? Procedure calls for elimi- 
nation of @ between expressions for components 
x, y, as given under Exercise 1. 

Exercise 5. What is the kinetic energy of twist 
of the pulley about a vertical axis, for a given h3? 
This is obviously small, but how small? 


NUMERICAL TABLE 


Assuming no friction in the present apparatus, 
Table I gives numerical values dependent on 
time, calculated in terms of total rotations 
from the start. Formulas will be found in the 
group designated Eq. (77) above. The time rela- 
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tion is that of Eq. (75). (See ‘‘Numerical Values 
of Parameters’ above.) 

Values of dependent parametric quantities are: 
(a) string tensions, by Eqs. (33) and (34): fi 
= 31,507, fo=31,969 dynes, (note that sum f1+ fe 
= 64.771 g is only slightly less (0.35 percent) 
than weight of pulley, 65.000 g, while f.2—/f,, 
Eq. (35), is only about 0.73 percent of mg); 
(b) torques, by Eqs. (40), (41), (46), and (47): 
on axle, 71= 83,493, ro= —71,931, and re=71+ 72 
= 11,563, which checks with 2/7,S; on pulley, 
Tp =21,S = 1160.0, checking with (fe—f1)r3, Eq. 
(41); and finally, (c) accelerations, by Eqs. (37), 
(38), (39), and (45): a,;=43.5, a2= —36.9, az 
= 3.28 cm/sec’; while a= 16.40 rad/sec?. 


POINTS OF INTEREST 


(1) Of primary importance to teaching is the 
fact that the apparatus is mechanically simple 
enough for the student to follow the working 
out, in a rather complicated way, of Newton’s 
second law, operating here in simultaneous, re- 
lated linear and angular motions. A consider- 
ably larger number of magnitudes is met than 
in the Atwood’s machine. 

(2) Mathematical relations include three dif- 
ferent-valued linear motions, having to do with 
accelerations, velocities, and displacements, and 
also corresponding to these the angular quan- 
tities of two different rotary units, one fixed, the 
other having translation. Although many equa- 
tions appear during the derivations, the working 
formulas are not too numerous and are quite 
simply expressed. 

(3) The abbreviating functions Q and S have 
physical meaning as well as mathematical utility. 

(4) Interpretation can be made in terms of an 
“equivalent axle,” the diameter of which equals 
the difference in radii of the drums of the com- 
pound axle. 

(5) The distance the pulley descends in one 
rotation is equal to the difference of radii of axle 
drums multiplied by z. In the present apparatus 
the pulley rotates nearly four times to descend a 
distance equal to its own diameter. 

(6) Although the tangential force on the small 
drum of the axle, which is the tension in the right 
string, Fig. 1, is greater than that on the large 
drum, which is the tension in the left string, the 
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torque on the small drum is Jess than that on the 
large one. The ratio of the two torques is equal 
to the ratio of the two tensions multiplied by the 
ratio of the radii. This product here is —0.986 
X1.18= —1.16, the tensions differing from each 
other by only 1.5 percent, while the radii differ 
by 18 percent, causing the torques to differ by 
15 percent, oppositely from the tensions. 

(7) Two visible characteristics of the me- 
chanical action: (a) the cord winds up on one 
side at the same time it unwinds on the other, 
giving somewhat the impression of a “‘self- 
acting’? mechanism, but an “externally applied 
force” is necessary; (b) the descending pulley 
remains in phase with the fixed rotating axle, 
analysis showing that pulley and axle have the 
same angular acceleration and velocity at every 
instant. 

(8) Along ~vith the Atwood’s machine and the 
inclined plane, this system offers a mechanical 
means of reducing acceleration of gravity, so 
that a/g can be varied. Here the acceleration of 
the pulley is about 3 percent of g. A smooth 
inclined plane giving this acceleration to a mass 
sliding down it, would need to be tilted from the 


horizontal only 11.5 mins of arc. Correspondingly 
in the Atwood’s machine, neglecting friction and 
moment of inertia, the mass of the rider would 


be 0.66 percent of the mass of one of the equal 
weights. 


CONCLUSION 


In general physics, “simple machines” and 
mechanical advantage are treated as elementary 
problems in equilibrium, any accelerations being 
usually omitted. These systems, especially for 
the college student, have taken on something of 
a “‘cut-and-dried”’ nature. Their basic importance, 
however, needs continued emphasizing. This 
paper proposes that the classical systems of 
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pulleys be dealt with from a nonequilibrium 
point of view. 

A particular pulley system is proposed and 
described. From the nonequilibrium standpoint 
it offers broad opportunities for instructional 
use in lecture, quiz, classroom, and laboratory. 
The paper analyzes the system mathematically 
while emphasizing the physics. The resulting 
working equations are summarized. The deriva- 
tions, in a particular sequence, keep the be- 
ginning student in mind. 

The student, finding himself here in what may 
seem at first a somewhat baffling experimental 
situation, by continued ‘looking beneath the 
surface” gradually recognizes an increasing num- 
ber of related mechanical facts, assuming he has 
had no information about the system beforehand. 
At the same time he observes in this example 
nature’s ability to ‘hang together,” avoiding 


entanglements and maintaining a “balanced” 


condition amid complexity, through an im- 
pressively consistent use of simply stated prin- 
ciples, especially that of Newton’s second law. 

Provocative questions, problems and exercises 
have been included for use by teachers. The 
plotting of graphs, like J,.=f(m), the numerous 
other quantities including time remaining the 
same, would add interest. Obviously, teaching 
possibilities have not been exhausted in the 
present paper. 

What spurred the present study was the re- 
quest made by the Secretary, Southern California 
Section of AAPT, for program papers preferably 
in the fields of alternating currents, Newton’s 
second law and rotation. The last named subject 
appealed particularly to this writer. 

Finally, thanks are due to Mr. Lewis H. 
Humason for his assistance in the mechanical 
design, and to the Physics Shop for construction 
of the apparatus. 


We, like the eagles, were born to be free. Yet we are obliged, in order to live at all, to make a cage 
of laws for ourselves and to stand on the perch. We are born as wasteful and unremorseful as tigers; 
we are obliged to be thrifty, or starve, or freeze. We are born to wander, and cursed to stay and dig. 


—WILLIAM BOLITHO. 
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The uniform model of the nucleus, proposed by Wigner in 1937, assumes that the motion of 
the nucleons within the nucleus are so complicated, and the interactions between them are so 
strong, that the details of the motion are of no great importance. This model, therefore, differs 
diametrically from the shell model of the nucleus, in which each nucleon is assumed to move 
independently of all others to a first approximation, the various nucleons filling up a regular 
sequence of possible orbits (shells). In this paper we re-examine the evidence for the uniform 
model to see whether, in the light of data accumulated since the time the model was proposed, 
the agreement between the predictions of the uniform model and experimental binding energies 
is sufficiently good to serve as an argument against the shell model. This is found not to be the 
case. The uniform model gives qualitative agreement with experiment, but the detailed quan- 
titative agreement is only moderately good. This rather indifferent record must be compared 
with the excellent evidence for nuclear shell structure in the same region of mass numbers. 
Nevertheless, we feel that the uniform model should not be abandoned altogether, but rather 
should be reinterpreted to apply to the “characteristic state” postulated by Hurwitz and 
Bethe, rather than to the ground state. No test of this hypothesis is possible until more is known 
experimentally about the level spectra of nuclei. In this connection a modified semiempirical 


form is suggested for the density of nuclear levels as a function of excitation energy. 


1. INTRODUCTION 


HE uniform model of the nucleus, proposed 
by Wigner,! assumes that the motion of 
the nucleons within the nucleus is so complicated, 
and the interactions between them are so strong, 
that the details of the motion are of no great 
importance. The interactions between nucleons 
in the nucleus are assumed to be given by ex- 
change forces of the Majorana (space exchange) 
type, with perhaps a small admixture of ordinary 
(Wigner) forces. Spin-dependent forces are as- 
sumed to be small enough so that they can be 
neglected in a first approximation. The contribu- 
tion of each pair of nucleons to the total poten- 
tial energy then depends only upon the type of 
the pair, i.e., upon whether the wave function is 
symmetric or antisymmetric under the exchange 
of the space coordinates of the two nucleons, and 
upon the nuclear radius. 

* Assisted by the joint program of the ONR and AEC. 

{+ Based upon a S.M. thesis submitted by one of us 
(I.G.W.) to the Massachusetts Institute of Technology. 

t Part of this work was done while both authors were at 
Massachusetts Institute of Technology. Additional re- 
search to complete the paper was performed while the 
second author was at the University of Illinois. 

1E. P. Wigner, Phys. Rev. 51, 106, 947 (1937); Bicen- 
tennial Symposium, University of Pennsylvania (1940). 


E. P. Wigner and E. Feenberg, Reports on Progress in 
Physics (The Physical Society, London, 1941), Vol. 8 
p. 244. 


Wigner! and Barkas? have shown that the 
theory based upon these assumptions gives ex- 
cellent agreement with experimental data on 
binding energies, especially for the differences 
between binding energies of neighboring isobaric 
nuclei. Margenau and Wigner* showed that the 
observed magnetic moments of nuclei could also 
be understood on the basis of the uniform model. 
Wigner‘ applied the theory to the calculation of 
the beta-decay constants (ft values) of ‘‘allowed 
and favored” beta-transitions in light nuclei, 
again obtaining excellent agreement. 

Thus, some ten years ago, the uniform model 
was considered very successful. However, various 
developments since that time have made a re- 
examination of the model desirable. First, Wig- 
ner! already called attention to the fact that the 
model underestimates the absolute amounts of 
the binding energies, even though it gives a good 
approximation to the difference between binding 
energies of neighboring isobars. Second, the dis- 
covery of the quadrupole moment of the deu- 
teron® and its interpretation in terms of a non- 


2 W. H. Barkas, Phys. Rev. 55, 691 (1939) 


3H. Margenau and E. P. Wigner, Phys. Rev. 58, 103 
(1940). 


4E. P. Wigner, Phys. Rev. 56, 519 (1939). 


5 Kellogg, Rabi, Ramsey, and Zacharias, Phys. Rev. 
57, 677 (1940). 
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central, spin-dependent tensor force® made the 
assumption of spin-independent forces somewhat 
doubtful. Recent calculations of Feingold’ on 
light nuclei indicate that the tensor forces make 
important contributions to the total binding 
energy. Third, the high energy neutron-proton 
and proton-proton scattering data obtained at 
Berkeley show that the basic law of force be- 
tween two nucleons is very much more compli- 
cated than had been believed. The scattering 
data seem to indicate that the nucleon-nucleon 
force does not have enough exchange character 
to lead to saturation in heavy nuclei.’ While 
these lines of evidence tend to weaken the basic 
assumptions of Wigner about the force law, the 
recent evidence for shell structure in nuclear 
ground states and in low-lying excited states® 
weakens the “uniformity”? assumption quite 
badly. The shell structure evidence indicates that 
the nucleons move, to a first approximation, in- 
dependently of each other in well-defined orbits. 
This is the exact opposite of the “uniform” 
assumption. 

Thus, it appears quite surprising that the uni- 
form model should be so extremely successful in 
predicting nuclear binding energies. A re-exami- 
nation of the evidence by Frisch, using more 
recent experimental data, did indeed indicate 
somewhat poorer agreement. The present investi- 
gation was undertaken to see to what extent the 
agreement of the uniform model with experiment 
can be considered proof of the ‘‘uniformity”’ as- 
sumption, and therefore as an argument against 
the usual explanation of nuclear shell structure in 
terms of an independent particle model of the 
nuclear motion. Since there are theoretical rea- 
sons for expecting the uniform model to break 
down for very light nuclei on the one hand (the 
uniformity assumption surely fails there) and for 
heavy nuclei on the other hand (Coulomb effects 


( ont) Rarita and J. S. Schwinger, Phys. Rev. 59, 436, 556 
1941). 

7A. M. Feingold and E. P. Wigner, Phys. Rev. 79, 
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8R. S. Christian and E. W. Hart, Phys. Rev. 77, 441 
(1950); R. S. Christian and H. P. Noyes, Phys. Rev. 79, 
85 (1950); E. Gerjuoy, Phys. Rev. 77, 568 (1950). 

°M. G. Mayer, Phys. Rev. 74, 235 (1948); 75, 1969 
(1949); 78, 16, 22 (1950). Haxel, Jensen, and Suess, Phys. 
Rev. 75, 1766 (1949). L. W. Nordheim, Phys. Rev. 75, 
1894 (1949); 78, 294 (1950). E. Feenberg and K. C. 
Hammack, Phys. Rev. 75, 1877 (1949). 

10D. H. Frisch, Phys. Rev. 84, 1169 (1951). 
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become very important), we have restricted our- 
selves to nuclei with mass numbers between 16 
and 65. The precise range chosen is somewhat 
arbitrary, but a line had to be drawn somewhere 
in order to avoid the temptation to include ‘‘just 
one more good case.’’ We have compared the 
observed binding energy differences between 
neighboring isobars and the total binding energies 
of nuclei within that range with the predictions 
of the uniform model. 

The result can be summed up as follows: The 
uniform model gives qualitative agreement with 
experiment, but a more quantitative comparison 
shows that the detailed agreement is only moder- 
ately good. In particular, the effects of magic 
numbers appear very clearly, and even for non- 
magic nuclei the uniform model has a tendency 
to underestimate the amount by which even-even 
nuclei are favored over odd-odd nuclei. The old 
disagreement with the total binding energies re- 
mains. This rather indifferent record of the uni- 
form model should be compared with the excel- 
lent evidence for nuclear shell structure in the 
same region of mass numbers. In our opinion, 
the comparison speaks against the uniformity 
assumption for the ground states of nuclei. 
Nevertheless, we feel that the uniform model 
should not be abandoned altogether, but rather 
should be reinterpreted to apply to the ‘‘charac- 
teristic level’’ postulated by Hurwitz and Bethe" 
rather than to the ground state. This suggestion 
is discussed in Sec. 5 of this paper. No definite 
conclusions can be reached until more is known 
experimentally about excited states of nuclei. 


2. THE BINDING ENERGY FORMULA OF 
THE UNIFORM MODEL 


In this section we compile the relevant for- 
mulas for easier reference. The total energy of a 
nucleus consists of three parts: the nuclear poten- 
tial energy V, kinetic energy 7, and electrostatic 
(Coulomb) potential energy C. The nucleus is 
specified by stating its mass number A, and its 
neutron excess T; defined by 


T;=}(N—Z), (2.1) 


where N is the number of neutrons and Z is the 
number of protons. The potential energy for a 


1H. Hurwitz and H. A. Bethe, Phys. Rev. 81, 898 
(1951). = : 
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pair of nucleons is assumed to be 
Vie= Vw (rie) + Va(riz) Pio™, (2.2) 


where Vy is the potential of the Wigner force, 
Vu gives the distance dependence of the Ma- 
jorana exchange force, and P™ is the Majorana 
operator which exchanges the positions in space 
of the two nucleons in question, without affecting 
the spin coordinates. 

We distinguish between ‘“‘symmetric’’ and 
“antisymmetric” pairs of nucleons (under the 
operation of space-exchange). We denote by Vyw* 
the average value of Vw(r;;) for symmetric pairs, 
and by Vyw* the corresponding average for anti- 
symmetric pairs. A similar notation is used for 
averages of the Majorana potential. Let n, be 
the number of symmetric pairs, n_ the number of 
antisymmetric pairs. Then the nuclear potential 
energy is given by 


V=n,Vw'+n_Vw*+n.Vut—n_Vy*. (2.3) 


This expression can be rewritten as 


V=—L(n,—n_)—L’(n,+n_), (2.4) 


where 
L=- 3( Vyri- Vw’) 
—3(Vau'+ Vu’) =L;/A, 
L!=—\(Vw'+ Vo") 
—3(Vau'—Vu*)=Li'/A. (2.6) 


The uniform model of Wigner assumes that the 
quantities L and L’ are smooth functions of the 
mass number A, independent of special properties 
of individual nuclei. Furthermore, both L and L’ 
can be estimated to be inversely proportional to 
A, so that the quantities L; and L,’ defined by 
Eqs. (2.5) and (2.6), respectively, are taken to 
be constants. While these constants may be esti- 
mated from our information about forces be- 
tween isolated nuclear particles, they are usually 
taken to be adjustable parameters in this theory. 
The signs in Eqs. (2.5) and (2.6) have been 
chosen in such a way that L and L’ are expected 
to be positive. 

The form of Eq. (2.4) is useful since the second 
term is proportional to the total number of pairs, 


n4+n_=3A(A—1), (2.7) 


and may hence be ignored in comparisons of bind- 
ing energies of isobaric nuclei. The crucial quan- 


(2.5) 
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tity in that connection is the difference n,—n_ 
between the numbers of symmetric and anti- 
symmetric pairs.” This quantity can be derived 
by group-theoretical methods. For any one nu- 
cleus it can assume many different values, the 
spectrum of possible values depending upon A 
and T;. According to Eq. (2.4) we are interested 
in the highest eigenvalue of »,—n_ consistent 
with a given mass number and neutron excess; 
for this gives the minimum potential energy V 
and, hence, is appropriate for the ground state 
of the nucleus. Physically this is reasonable since 
the Majorana force is attractive for symmetric 
pairs, repulsive for antisymmetric pairs, so that 
large values of n; —n_ give large binding energies. 
This highest eigenvalue of n,—n_" is given by 


n4—n_= —(A?/8)+2A+(5/2)—§ (2.8) 


where ~ depends upon the type (even-even, odd- 
odd, even-odd) of the nucleus in question: 


§=37T2+2|T;|+(13/4) odd A, 
&€=472+2|T;|+(5/2) even A, even Z, 


§=17242|T;| +4 even A, odd Z, 
except Z=N, (2.9c) 


even A, odd Z, 
with N=Z. (2.9d) 


We observe that ¢ is lowest for | T;| close to zero. 
This is in agreement with the empirical rule that 
nuclei prefer to have equal numbers of neutrons 
of protons, except for the Coulomb effect which 
favors a neutron excess. In a similar way, the 
following three empirical rules are also direct 
consequences of Eqs. (2.8) and (2.9): (1) The 
binding energy of odd-A isobaric nuclei is a 
smooth function of the neutron excess. T;; (2) 
even-even nuclei are more strongly bound than 
odd-odd nuclei; (3) odd-odd nuclei with N=Z 
form an exception to rule (2), several of them 
being beta-stable. Thus, Wigner was able to 
show that the assumption of spin- and charge- 
independent forces of predominantly Majorana 
exchange character together with a uniform 


(2.9a) 
(2.9b) 


g=5 


12 This quantity is called = in reference 1. The & of Eqs. 
(2.8) and (2.9) is Wigner’s =’. 

18 Under Wigner’s assumptions about the nuclear forces 
(charge- and spin-independent), and with the further as- 
sumption that the Coulomb repulsions between the protons 
are negligible in a first approximation, the operator n,—n_ 
commutes with the Hamiltonian H of the nucleus. 





UNIFORM MODEL OF THE NUCLEUS 


model for the nuclear motion provides a natural 
explanation not only of saturation, but also of 
the empirical rules for nuclear stability enumer- 
ated above. 

The relations (2.5)—(2.9) contain the estimate 
of the uniform model for the nuclear potential 
energy. The kinetic energy T is estimated by 
assuming that the nucleons form a Fermi gas 
confined to a sphere of the nuclear radius 


R=roAt=1.47A!X10-" cm. = (2.10) 


The kinetic energy of the Fermi gas also de- 
pends upon the neutron excess to some extent, 
since the height of the Fermi distribution is in- 
creased by an increased number of identical par- 


ticles. For |7T;|<«<A, we get the approximate 
relation 


T=T,A+(Ti/A)nt+::-, (2.11) 


where Jo, 71, and 7» are defined as follows 
(M =nucleon mass): 


To = (3/40) (92) #(h?/ Mr?) 13.5 Mev, 
T,= (40/9) T>=60 Mev, 
n=3Tf+i odd A, 
=43T; 
n=3T +3 


(2.12) 
(2.13) 
(2.14a) 
(2.14b) 


even A, even Z, 


even A, odd Z, 


except Z=N, (2.14c) 


(2.14d) 


The numerical estimates in Eqs. (2.10), (2.12), 
and (2.13) are derived from the study of the 
Coulomb energy difference between mirror nu- 
clei, which depends upon the nuclear radius. The 
value 1.47 X10-" cm for 7o is also in reasonable 
agreement with data from nuclear reactions. The 
terms neglected in Eq. (2.11) are of order 
(T;/A)? or smaller. They are negligible compared 
to the intrinsic error of the Fermi gas approxima- 
tion. Indeed, the estimate of the kinetic energy 
is probably one of the weaker points of the 
theory. One must remember that about half the 
potential energy of attraction in the nucleus is 
compensated by the kinetic energy, so that an 
error in the estimate of the kinetic energy is very 
serious indeed for the binding energy. 

Since T, in Eq. (2.11) is positive, nuclear sta- 
bility is enhanced by a low value of 7. Equations 


n= even A, odd Z, Z=N. 
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(2.14) show that » has properties which are in 
many ways similar to , Eq. (2.9). 

Finally, we need an expression for the Coulomb 
energy of the nucleus. In accordance with the 
“uniform’’ assumption, the Coulomb energy is 
estimated as that of a uniform distribution of Z 
protons over a sphere of the nuclear radius R, 
Eq. (2.10), i.e., 


C=}32(Z—1)(6e?/5R) =4u,Z(Z—1)A-3. (2.15) 


Collecting the expressions (2.4), (2.7), (2.8), 
(2.11), and (2.15), we can write the total energy 
E(A, T;) of a nucleus as the sum of two terms, 
the first of which depends only upon the mass 
number, while the second contains the stability 
relations between isobaric nuclei: 


E(A, T;)=Eo(A)+E\(A, T>), 
where E,(A) can be chosen to be 


E,(A) = L,(A/8—2—5/2A) 
—$L;'(A—1)+TA+u,A*8(1—2/A). (2.17) 


With this choice of Eo(A), Ei(A, Ty) becomes 
E\(A, T:) = (Lig+ Tn)/A 


—4u,A-'T;(A —1—T)). 


(2.16) 


(2.18) 


In this way of writing the nuclear energy, 
there are five adjustable parameters: Ly, Ly’, To, 
T,, and u,. Wigner! took the following approach: 
u- was taken from the Coulomb energy differ- 
ences between mirror nuclei ( and y in Eq. (2.18) 
are the same for two mirror nuclei, 7; and — Ty), 
and To, JT; were taken from the Fermi gas ap- 
proximation for the kinetic energy. This left only 
two parameters, L; and L,’, only one of which 
(namely, L;) was of importance for the deter- 
mination of the stability relations between iso- 
bars. We shall allow ourselves more ‘latitude in 
fitting the data: we also take u, from the mirror 
nuclei, but we use the Fermi gas approximation 
for the kinetic energy only qualitatively, i.e., 
in so far as it determines the form of the depend- 
ence of kinetic energy upon A and T;. The quan- 
tities Ty) and 7, will be treated as adjustable 
parameters. According to Eq. (2.18), once u, is 
fixed, we have two adjustable parameters (Li 
and 7;) to fit the observed stability relations 
between isobars. We shall show that the present 
data cannot be fitted very well by formula (2.18) 
no matter what values are assumed for L; and 7. 
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Fic. 1. Ordinate: The nuclear energy difference AE, 
between two neighboring isobars (corrected for the Cou- 
lomb effect and for the neutron-proton mass difference) 
in Mev. Abscissa: The inverse of the mass number. The 
figure contains only those isobaric pairs which have odd 
mass number A (g—u) and whose two members have 
neutron excesses Ty=} and Ty=}, respectively, (7, =). 
Circles represent measured values, semicircles with arrows 
represent upper or lower limits, the arrows pointing into 
the allowed region. Double circles represent magic pairs 
(i.e., one or both of the members of the isobaric pair is 
magic); the particular magic number involved is indicated 
next to the point. The shaded region is excluded by nuclear 
stability data. According to the uniform model, these 
points are expected to fall on a straight line passing 
through the origin, of slope 


s=L,At+TiAn, 


where LZ; and 7; are empirical parameters to be deter- 
mined, and Aé and Ay are given in Table I. 


If a fit to the formula is nevertheless attempted, 
the values used by Wigner (L,;=43 Mev, T,=60 
Mev) are not the best possible ones. Lower values 
of L; and higher values of T, seem indicated. 


3. ANALYSIS OF BINDING ENERGIES OF 
ISOBARIC NUCLEI 


The basic assumptions of the uniform model 
restrict the range of validity of Eqs. (2.16)- 
(2.18) to a rather narrow range of mass numbers. 
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The model does not apply to very light nuclei, 
since the ‘‘uniformity’’ assumption is almost cer- 
tainly wrong there. it also cannot be expected to 
apply to nuclei of mass number more than 50-70. 
For A> 60, the Coulomb repulsion between the 
protons becomes so significant that it can hardly 
be treated as a small perturbation. The basis of 
the group-theoretical calculation of n,—mn_ (the 
charge- and spin-independence of the main forces 
in the nucleus) then disappears, and the model 
fails to give correct answers. Thus, any experi- 
mental test of the uniform model must be re- 
stricted to the intermediate range of mass num- 
bers, say between 16 and 65. 

The experimental data on stability relations 
between isobaric nuclei are mostly derived from 
B-decay data. Since there are always two nuclei 
involved in a #-transition, we need a uniform 
notation. Let the 8-transition involve the two 
nuclei (A, T;) and (A, T;’). (The direction of the 
B-decay is of no consequence here.) We shall de- 
note the larger one of T;, T;’ by Ty. Any §-transi- 
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Fic. 2. Same as Fig. 1, for isobaric pairs of 
type g—u, T,=5/2. 
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tion is specified (except for the direction of the 
decay) by the mass number A and by T3. 

For any isobaric pair of nuclei, (A, T,) and 
(A, T.—1), we are interested in the energy differ- 
ence between their ground states. Thus, if the 
B-decay proceeds not to the ground state of the 
daughter nucleus but to an excited state, we need 
to know the energy of the subsequent y-ray (or 
y-rays) in addition to the maximum energy in 
the B-spectrum. 

We have tabulated these transitions in the 
Appendix. For each transition, we list the mass 
number and atomic number, the type of decay 
(8-, B+, K capture, y-rays), the kinetic energy 
release AE’ (=the maximum kinetic energy in 
the 8-spectrum plus the energy of the subsequent 
y-rays), the Coulomb energy difference between 
the two nuclei [computed from Eq. (2.15) with 
ue=0.147 Mev], and the binding energy differ- 
ence corrected for the Coulomb energy difference 
and the neutron-proton mass difference. This 
latter quantity is called AE, in the tables. It 
represents the difference in the nuclear binding 
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Fic. 3. Same as Fig. 1, for isobaric pairs of 
type g—u, T,=7/2. 











Fic. 4. Same as Fig. 1, for isobaric pairs of 
type g—g, T>=1. 


of the two isobars. According to Eq. (2.18), AE, 
is given by 
AE, = (LiAE+TAn)/A, (3.1) 


where Aé is the difference between (A, T,) and 
(A, T,—1), computed from Eq. (2.9), and An 
is the corresponding quantity for y, Eq. (2.14). 

We have classified the energy release data ac- 
cording to the ‘“‘type’’ of the isobaric pair. By 
‘“‘type’’ we mean the value of T,. together with 
the specification of the even-odd, even-even, or 
odd-odd character of the nucleus (A, T,).!* Iso- 
baric pairs with odd mass number A are classified 
as g—u; isobaric pairs with even mass number A 
are classified as g—g if the nucleus (A, T5) is 
even-even, and are classified as u—u if the 
nucleus (A, T;) is odd-odd. 

This classification is useful because all isobaric 
patrs of the same type have the same values of Ag 
and An. For example, isobaric pairs of type u—u 
with T,=1 have A§=13/2—5/2=4 and An=1 

We emphasize once more that the direction of the 
B-decay does not concern us here. An isobaric pair (A, T>) 
and (A, T,—1) is classified as g—g if the nucleus with the 


larger neutron excess, (A, T,), is an even-even nucleus, no 
matter whether this nucleus is #-stable or B-unstable. 
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Fic. 5. Same as Fig. 1, for isobaric pairs of 
type g—g, T>= 


—0Q=1, according to Eqs. (2.9) and (2.14), re- 
spectively. Returning to Eq. (3.1), we see that 
AE, is expected to be proportional to A-! for 
isobaric pairs of the same type, with proportion- 
ality constant equal to (L;AE+7,An). Thus, if we 
plot the experimental values of AE, against the 
inverse of the mass number, the points ought to 
fall on a straight line passing through the origin, 
of slope s= L,AE+7,An. 

Such plots are shown for the relevant types 
(nuclei in the mass number range 16<A <65) in 
Figs. 1-9. It is apparent from these figures that 
the assumptions of the uniform model are rather 
poorly satisfied: The points show considerable 
scatter about a straight line, indeed about any 
smooth curve. 

While this is perhaps unfortunate, it is by no 
means unexpected and can hardly be considered 
as a significant failure of the model. After all, 
any such over-all, global approach necessarily 
neglects many detailed features of individual 
nuclei which affect their binding energies to some 
extent. Fluctuations of the points about any 
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smooth curve must be expected for that reason. 
The most the uniform model can hope to ac- 
complish is to give the general trend correctly, 
neglecting the individual fluctuations. 

We have therefore adopted the following point 
of view: We use all the data on each figure to 
determine a value for the slope s of the “‘best”’ 
straight line (shown as a solid line on each figure), 
and we take into account the fluctuations by 
assigning an error to this slope. That is, we draw 
two additional lines on each figure, corresponding 
to our estimate of the highest and of the lowest 
slope still consistent with the main trend of the 
data, and consider their slopes as upper and lower 
limits of the ‘“‘best’’ value of s, respectively. 
These limiting lines are shown dashed on each 
figure. One should then expect that these ‘‘best”’ 
slopes are represented correctly by the theo- 
retical formula s=L,AE+7),An, within their ex- 
perimental errors. 

In determining the ‘‘best’’ straight line, it is 
not necessary to rely exclusively upon measured 
energy differences between isobaric pairs. Con- 
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Fic. 6. Same as Fig. 1, for isobaric pairs of 
type g—g, T>=3. 
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siderable information can be obtained from par- 
tial information about the decay scheme. For 
example, the information that a 6- decay with 
a certain maximum energy has been observed 
puts a lower limit on AE,, even if the subsequent 
y-radiation has not been analyzed. Finally, the 
very information that certain nuclei are stable 
puts upper or lower limits, as the case may be, 
upon A£,. On the figures, points which represent 
upper or lower limits are shown with arrows, 
while regions of the graph which are excluded by 
the known stable nuclei are shown shaded. These 
data have also been included in the tables in 
the Appendix. Liberal use has been made of this 
additional information in order to decide upon 
the ‘‘best’’ slope. 

In principle, we should exclude from these 
figures all isobaric pairs for which one or both of 
the nuclei involved are ‘‘magic.’’ Magic pairs are 
shown on the figures by double circles. We see 
that the density of ‘‘magic numbers” is so high 
that about half of all the pairs are magic. Thus, 
exclusion of the magic pairs would have de- 
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Fic. 7. Same as Fig. 1, for isobaric pairs of 
type u—u, T,=1. 
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Fic. 8. Same as Fig. 1, for isobaric pairs of 
type u—u, T,=2. 


creased the statistics very badly. However, the 
magic is a random effect in the sense that either 
of the two nuclei is as likely to be magic as the 
other one, other things being equal. Thus, on the 
average, for any magic pair with (A, T,) magic, 
there is another magic pair with (A, T,—1) 
magic. In the first case, AE, is unusually low, in 
the second it is unusually high. To the extent 
that there are equally many cases of each kind, 
the slope of the best straight line is not affected 
by the inclusion of magic pairs. The error in the 
slope is increased (i.e., the dotted lines have to 
be drawn farther away from the best line), but 
this is perhaps a healthy influence, since the 
usual tendency is to overestimate the accuracy 
to which the slope is known. 

We would like to point out that, as an inci- 
dental by-product of this investigation, we obtain 
a rough estimate of the energy contributed by 
the various magic numbers. We see that the 
magic number 20 leads to an increase of the 
binding energy of approximately 1 Mev, and the 





ISABELLA GOLDIN WEINBERG AND J. M. 




















| 

f 
r 
| 

| 
0 


u, T>*3, 


| 


el 
| 60 50 
015 02 








Fic. 9. Same as Fig. 1, for isobaric pairs of 
type u—u, T,=3. 


magic number 28 leads to a somewhat smaller 
increase (perhaps around 4 Mev), whereas the 
so-called magic number 14 does not show up 
at all. 

The best slopes together with their errors are 
compiled in Table I. In most cases it was suffi- 
cient to assign the same error to the slope in both 
directions. But for some of the figures the upper 
and lower limits on s are significantly unsym- 
metrical about the best value of s; in these cases 
the upper and lower limits are stated separately. 
Also given in Table I are the values of Aé and 
An for each type of transition. 

We must now determine how we are to test 
these values of s for consistency with the uniform 
model. One particular case is simple to test: the 
g—g transitions with T,=1. The agreement is 
excellent. Having disposed of this case, we have 
8 values of s which must be fitted by 2 inde- 
pendent parameters (ZL; and T;). Since there are 
only two parameters to be fitted, it is possible to 
plot the data in such a way that these parameters 
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are related to the slope and intercept of a straight 
line. To do this, we define the quantities x and y 
as follows: 


g—u transitions, 
g—g transitions, 
u—u transitions, 


= 3(T>+9) 
t= 6T, 
x=2(T>+3) y=2s. 

Then the following equation holds for all 
transitions (except the g—g, T.=1 type): 


y=(L1+T1)x—-6T}. 


y=3s; 


y=6s; (3.2) 


(3.3) 


Thus, a plot of y against x should be a straight 
line of slope L1+T7; and of intercept —67,. The 
points are shown in Fig. 10. We see that the 
points do not fall very well on any straight line. 
There is a tendency for the g—g points to lie 
lower than the u—wu points by appreciable 
amounts. This means that the binding energy 
formula of the uniform model underestimates the 
difference in binding energy between even-even 
and odd-odd nuclei. This difficulty was already 
emphasized by Wigner in connection with the 
observed existence of triple isobars which would 
not be expected according to the uniform model. 
Figure 10 shows that the contradiction with ex- 
periment is not absolute: It is still possible to 
draw a straight line (the solid line shown on the 
figure) which passes through every point within 
the assigned errors. However, the appearance of 
the figure makes it likely that more, and more 
accurate, data will destroy even that bare 
agreement. 

Let us now look at the intercept of the line on 
the vertical axis, which should be equal to — 67}. 
The solid line leads to the value 7,;=109 Mev, 
compared to the theoretical estimate of T,=60 
Mev. This is not very good agreement. To show 
that this is not just an accident, we have also 


TABLE I. Best values of s=L,AE+7,Ap. 


Type of transition An 


g—u T,=3/2 
Ty =5/2 


s (in Mev) 


188+18 
306+35 
490+30 


0+11 
103(+60, —20) 
265+55 


236(+60, —30) 
390443 
530+50 
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drawn in the best line with the theoretical value 
of 7; (i.e., with intercept —360 Mev). This 
(dashed) line gives an appreciably poorer fit to 
the points, missing two points altogether. 

The situation is improved somewhat if we take 
the attitude that the uniform model breaks down 
for A250 so that nuclei with A>50 should be 
excluded from this analysis. The points at x =15 
and x=18, as well as the triangle (u—x) point 
at x=12, depend predominantly upon data for 
mass numbers between 50 and 65. If these points 
are omitted, the dashed straight line (with the 
theoretical value of 7,) goes through all the 
remaining points. 

The best straight line on Fig. 10 corresponds 
to the parameters 


L,=30 Mev, 7T,=109 Mev (best line), (3.4) 


whereas the dashed straight line has the pa- 
rameters 


L,=47 Mev (with T,=60 Mev taken 


from the Fermi gas model). (3.5) 


As is to be expected, the value L;=47 is in 
good agreement with Wigner’s value (Li=43 
Mev). We see that an independent adjustment of 
L, and T; tends to raise TJ; and lower L, com- 
pared to the values with 7, fixed theoretically. 

This is an unfortunate tendency from the 
point of view of the uniform model. The main 
idea behind the model is to explain the differences 
between energies of isobars by exchange forces 
together with the Pauli principle (which deter- 
mines the numbers of symmetric and antisym- 
metric pairs). Thus, the model is constructed so 
as to place all the emphasis upon the potential 
energy. The kinetic energy is treated as an after- 
thought, as it were, by an admittedly poor ap- 
proximation. This procedure can be justified a 
posteriort provided it turns out that the main 
contribution to the energy differences between 
isobars really comes from the potential energy. 
But this would mean L;>T7,, a condition which 
is violated already for Wigner’s values, Eq. (3.5), 
and is badly violated for the best values, Eq. 
(3.4). In view of these results it appears unten- 
able, in an explanation of isobaric energy differ- 
ences with a uniform model, to place the main 
emphasis upon the potential energy. On the con- 
trary, the kinetic and potential energies should 
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Fic. 10. Test of the formula s=Z,AE+7 An for the 
slopes of the preceding figures. The ordinate y and abscissa 
x are defined by Eq. (3.2) and the points are expected to lie 
on a straight line given by Eq. (3.3). The fact that the 
points from g—g transitions (the squares on the figure) 
fall systematically lower than the best line, while the points 
from u—u transitions (the triangles on the figure) lie some- 
what high indicates that the uniform model underestimates 
systematically the energy difference between even-even 
and odd-odd nuclei. The solid straight line represents the 
best fit (obtained by a modified least-squares method), 
while the dashed straight line represents the best fit subject 
to the condition that 7,;=60 Mev, the theoretical estimate 
of that quantity. 


be treated on an equal footing, with equally good 
approximations. No such treatment is available 
at this time. 

One might object that we have completely 
ignored a great success of the uniform model: 
The zero slope for the transitions of type g—g, 
T;=1 in Table I. While this is undoubtedly a 
success for the uniform model, it should be 
realized that this equality of nuclear binding for 
A=4k+2 (k integral) nuclei with 7;=0 and 
T;=1 also follows from any shell model: These 
nuclei differ from each other merely through the 
fact that a proton in the top shell has been trans- 
formed into a neutron in the same shell. Thus, 
this success of the uniform model does not count 
in a comparison between the uniform model and 
the various shell structure (independent particle) 
models of the nucleus. This comparison, however, 
is our main point of interest, and it is for that 
reason that we have given little weight to this 
agreement between theory and experiment. 


4. THE TOTAL BINDING ENERGY AND 
NUCLEAR SATURATION 


Wigner already pointed out that the uniform 
model fails seriously in one respect: The model 
is unable to account for the magnitude of nuclear 
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binding energies if the parameters are adjusted 
to fit the binding energy differences between iso- 
bars. Referring back to Eq. (2.17), we see that 
the potential energy per nucleon is (in the limit 
of large mass numbers, A>1) 


V/A =L,/8—Ly,'/2+order A“. = (4.1) 


This quantity has to be negative (attractive), 
and must be sufficiently large in absolute value 
to overcompensate the kinetic energy and the 
Coulomb energy. The analysis of the energy 
differences between isobars gives an experimental 
value of L;. Wigner pointed out that Z,’ cannot 
exceed L, if the nuclear forces obey the assump- 
tions of the uniform model, i.e., are charge- and 
spin-independent and give saturation. If we use 
the theoretical values 7,;=60 Mev, 7 )=13.5 
Mev for the kinetic energy terms, we get [see 
Eq. (3.5) ] L1=47 Mev, and hence, L;’<47 Mev. 
The potential energy per nucleon in the nucleus 
is then less than (3/8)(47) =17.6 Mev. We sub- 
tract the kinetic energy per nucleon, T)=13.5 
Mev, to get a remainder of only 4.1 Mev. This 
would not be enough to account for the observed 
binding energies (of the order of 8 Mev per 
nucleon), even if there were no Coulomb repul- 
sion. Actually, for Z=20, A =40, the Coulomb 
energy amounts to about 1.6 Mev per nucleon, 
so that the binding energy per nucleon predicted 
by the uniform model turns out to be less than 
2.5 Mev. This is in clear contradiction of the 
experimental data. 

This discrepancy becomes even worse if we 
use the value L,=30 Mev obtained from the 
solid line of Fig. 10. With this value of LZ; and 
any reasonable estimate for the kinetic energy 
per nucleon To, we get no binding at all for the 
nucleus. 

Wigner (private communication) has shown 
that the difficulty with the binding energy can 
be overcome by taking the tensor force into ac- 
count. The tensor forces between the nucleons 
contribute an appreciable fraction of the binding 
energy, perhaps enough to explain the whole dis- 
crepancy. However, the inclusion of tensor forces 
into the calculation also has the effect of decreas- 
ing the expected energy difference between even- 
even and odd-odd nuclei. Since the conventional 
uniform model, using only central forces, already 
gives too small a value for this difference, the 
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inclusion of tensor forces does not really help 
very much. 


5. DISCUSSION 


While the uniform model of the nucleus could 
muster overwhelming theoretical and experi- 
mental support some 10 years ago (with the one 
exception of the discrepancy in the absolute 
binding energies) the situation is nearly reversed 
today. On the theoretical side, all the assump- 
tions underlying the uniform model are in doubt: 
The Wigner assumption about the nuclear forces 
is contradicted by the tensor force as well as by 
the two-nucleon scattering data and the conse- 
quent difficulties of explaining saturation by ex- 
change forces; the uniformity assumption about 
nuclear structure is contradicted by the mount- 
ing evidence for shell structure in nuclei. On the 
experimental side, this investigation, as well as 
others,!° has shown that the agreement of the 
binding energy formula with experiment is only 
moderate, and in some instances (total binding 
energy) rather poor. Thus, we feel that there is no 
reason to consider the ‘‘success’’ of the uniform 
model as an argument against shell structure 
theories based upon an independent particle model 
of low-lying nuclear states. 

Nevertheless, we do not think that the uni- 
form model should be abandoned altogether. In 
the remainder of the discussion, we would like 
to make some tentative suggestions for modifying 
the uniform model in such a way that agreement 
with experiment can be achieved. 

Weisskopf'® has pointed out that the evidence 
for nuclear shell structure is restricted entirely 
to low-lying (not very highly excited) nuclear 
states. The evidence is derived from binding 
energies, magnetic moments, spins, isomerism, 
and beta-decay, all of which are properties of 
low-lying states. Hurwitz and Bethe! have em- 
phasized that the fast neutron-capture cross 
sections are inconsistent with the assumption 
that the entire level system of the nucleus is 
affected by the ‘‘magic” of the shells. Rather, 
they postulate the existence of a ‘‘characteristic 
level’”’ whose binding energy is a smooth function 
of the mass number and shows no shell effects. 
Levels below the characteristic level are subject 
to shell effects, while more highly excited 


1 VY, F. Weisskopf, Helv. Phys. Acta 23, 187 (1950); 
Science 113, 101 (1951). 
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levels are assumed to show no evidence of shell 
structure. 

Hurwitz and Bethe propose that the binding 
energy of the characteristic level is given to a 
first approximation by the Weizsacker semiem- 
pirical formula, without the ‘‘odd-even’’ term. 
We would like to suggest that the energy of the 
characteristic level is instead given by the 
“‘sround-state’”’ energy formula (2.16)—(2.18) of 
the uniform model, which does include an ‘‘odd- 
even” term. The argument for this proposal is 
purely qualitative and very tentative only. We 
imagine that the highly excited states of a 
nucleus are so complicated in their structure that 
the uniform model gives a better approximation 
than any ordered model (e.g., independent- 
particle or alpha-particle model). The spectrum 
of these ‘‘uniform’”’ states extends down to the 
‘characteristic’ level, which is therefore effec- 
tively the ground state as far as the uniform 
model is concerned. The few states below the 
characteristic level are highly ordered states for 
which the independent-particle model gives a 
good first approximation. Shell effects are impor- 
tant for those states, and the uniform model is 
not applicable to them. 

This suggestion presupposes that the charac- 
teristic level, which separates these two regions, 
can be identified with some degree of accuracy. 
This could be done provided that there is a sharp 
break in the density of levels as the excitation 
energy increases above the characteristic level. 
Usually it is assumed that the density of levels 
w(E) is a smoothly increasing function of the 
excitation energy E, of type 


w(E) =C exp(2{aE}?), (5.1) 


where the constants C and @ are adjusted so as 
to reproduce the observed level densities near 
the ground state and near the separation energy 
of the first neutron (where the level density is 
known experimentally from a study of slow- 
neutron resonances). Expression (5.1) is based 
upon a statistical, thermodynamic treatment of 
the nucleus. This treatment may well be valid 
for levels which can be described by a uniform 
model. The treatment is surely invalid for states 
whose motion is so highly ordered that it can be 
described by an independent-particle (shell) 
model. These levels are effectively at ‘zero en- 


LEVEL DENSITY in (Mev) 


ENERGY in Mev 
a a 2 6 


Fic. 11. Schematic picture of the level density w(EZ) as 
a function of the excitation energy E. Curve 1 shows the 
usually assumed behavior, formula (5.1), with constants 
fitted at E=1 Mev and at E=8 Mev. The very sharp slope 
of w(EZ) near E=0 is usually considered spurious and is 
therefore ignored; the extrapolated value of w(Z) at E=0 
would be w(0)=3.5X10-? (Mev), an impossible value. 
Curve 2 shows the result of formula (5.2); the energy of 
the characteristic level was assumed to be E,=3 Mev; the 
sharp discontinuity in slope at this energy is undoubtedly 
not realized in nature. Curve 3 shows the result of formula 
(5.3); this is qualitatively the same as formula (5.2), except 
that the discontinuity has been smoothed out arbitrarily. 
This figure shows that even though the constants in for- 
mulas (5.1) and (5.3) have been fitted so as to give the 
same level densities near the ground state and at 8 Mev, 
these two formulas predict quite different behavior of the 
level density in the intermediate region. 


tropy”’ and should be excluded from thermo- 
dynamic considerations. While we have no clear 
idea of the behavior of the level density in this 
zero entropy region, it is surely not too far from 
a uniform density of levels, at least qualitatively. 
We are therefore led to consider the following 
level density formula: 


w(E)=C for 0<E<E,, 
w(E) =C exp[2{a(E—E.)}*] for E>E., 


where E, is the excitation energy of the character- 
istic level and C and a are adjustable constants 


(5.2) 
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which are to be used to fit the observed level 
densities near the ground state and near the 
separation energy of the first neutron. 

Formula (5.2) has one qualitative defect: 
There is an extremely sharp break in the level 
density at the energy of the characteristic level. 
Indeed, the derivative dw/dE is zero if E ap- 
proaches E, from below, infinite if E approaches 
E. from above. There is no reason to believe that 
the distinction between the two regions (above 
and below £,) is that sharp. We therefore need a 
formula which approaches a constant value for 
E<E,, and approaches the value given by Eq. 
(5.2) for E>E,. The following is suggested as a 
plausible interpolation formula: 


a(E—E,) i 
—~——___| . (5.3) 
1—exp[—a(E—E,) ] 


w(E)=C exp | 


The differences between formulas (5.1), (5.2), 
and (5.3) are illustrated in Fig. 11. In this figure 
the level densities are adjusted to the same value 
at E=8 Mev, and Eggs. (5.2) and (5.3) are ad- 
justed to a common value at E=0. Since it has 
been the conventional practice to adjust the 
constants in Eq. (5.1) to fit the density of low- 
lying levels at E=1 Mev rather than at E=0, we 
have adjusted the constants in Eq. (5.1) so that 
this function has the same value at E=1 Mev 
as the other two functions assume at E=0. The 
energy of the characteristic level was taken arbi- 
trarily to be E.=3 Mev. 

We remark that the conventional practice of 
ignoring the sharp energy dependence of Eq. 
(5.1) near E=0 conceals a difficulty: accord- 
ing to its derivation, the exponent 2{aE}} 
represents an entropy, and its energy derivative 
{a/E}* represents the inverse of the tempera- 
ture of the nucleus. Thus, if we replace the theo- 
retical level density Eq. (5.1) near E=0 by a 
more slowly varying function (which is in better 
agreement with experiment), we are decreasing 
the energy derivative of the entropy drastically. 
Since experimentally the level density near zero 
energy is a slowly varying function of E, the 
thermodynamic interpretation of the nuclear 
level structure leads to the paradoxical result 
that zero energy of the nucleus corresponds to a 
very high temperature. The paradox is resolved 
of course by remarking that thermodynamic con- 
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siderations cannot be applied to these low-lying 
levels, so that the exponent should not be inter- 
preted as a nuclear entropy in that region. 

There are various indications that Eq. (5.1) is 
not adequate to describe the level spectrum of 
actual nuclei. The most direct experiments on 
level densities are probably those of Gugelot!* on 
the spectra of secondary neutrons from nuclear 
reactions. Gugelot adjusts both C and a in Eq. 
(5.1) to fit his data near the separation energy of 
the first neutron. Then he observes that the 
extrapolated value of the level density near E=0 
is much too low, i.e., the extrapolated energy 
spacing of the low-lying levels is very much 
larger than the observed one. It is apparent that 
a level density formula such as Eq. (5.3) would 
give a better fit. However, we do not believe that 
we are justified in fitting Eq. (5.3) to Gugelot’s 
data at this time, since it has not been demon- 
strated experimentally that the sharp break in 
w(E) at E=E, predicted by Eq. (5.3) really 
exists in nature. 

The excitation energy E, of the characteristic 
level is not expected to be a smooth function of 
N and Z. Energy E, is the difference between the 
binding energy of the ground state of the nucleus 
(which is subject to shell effects) and the binding 
energy of the characteristic level. Only the latter 
quantity is expected to be a smooth function of 
N and Z, given by Eqs. (2.16)—(2.18). In par- 
ticular, we expect E, to be abnormally large 
(perhaps 4 to 6 Mev) in ‘‘magic’’ nuclei, com- 
pared to its value (perhaps 2 to 4 Mev) in 
ordinary nuclei. 

Nothing can be said about either the validity 
of this whole approach, or the numerical values 
of E., until much more is known about level 
spectra of medium-weight and heavy nuclei. We 
have given Eq. (5.3) not because we are sure it 
is correct, but because we believe that Eq. (5.1) 
is incorrect and some alternative formula is 
necessary, if only as a stopgap until the level 
spectra become better known. It is perhaps not 
sufficiently realized just how important the level 
density formula is for many of the common esti- 
mates in the theory of nuclear reactions, such as 
the ratio of widths for emissions of different 
particles. In general, formula (5.3) favors emis- 
sion of neutrons from the compound nucleus over 
~ 16 P. C. Gugelot, Phys. Rev. 81, 51 (1951). 
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two-neutron emission, and emission of a neutron 
and a proton, for example, is extremely depen- 
dent upon what is assumed for the level density 
in the various residual nuclei. 

At the very worst, Eq. (5.3) is likely to give a 
better approximation to the truth than Eq. (5.1) 
merely because it contains one more adjustable 
parameter (the excitation energy of the charac- 
teristic level). 
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APPENDIX: TABLES OF TRANSITION ENERGIES 


These tables give the energy differences be- 
tween isobars, organized according to the ‘‘type”’ 
of the isobaric pair. Symbol Ty is the value of 
the neutron excess T;=3(N—Z) for the isobar 
with the Jarger neutron excess (independently of 
whether this is the radioactive or the stable mem- 
ber of the isobaric pair in question). Symbol g—u 
indicates that the mass number is odd; g—g 
indicates that the isobar (A, T,) is an even-even 
nucleus; u— indicates that the isobar (A, Ts) 
is an odd-odd nucleus. 

In each table, the first two columns give the 
mass number and atomic number of the nucleus 
from which the relevant information is derived 
(usually the radioactive nucleus). The next two 
columns give the type of radioactivity observed 

17Way, Fano, Scott, and Thew, “Nuclear Data,” 


Circular of the U. S. National Bureau of Standards 499, 
September 1, 1950. 
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(st =stable) and the magic numbers involved, if 
any. The fifth column gives the kinetic energy 
release AE’ in the radioactive decay. This quan- 
tity is defined as the maximum kinetic energy of 
the decay electron (or positron) plus the energy 
subsequently emitted in the farm of y-rays. In 
the case of K capture (where the information is 
derived from nuclear reaction data rather than 
from B-decay data) AE’ is the total energy differ- 
ence and is equal to the energy given to the neu- 
trino. The fifth column gives the Coulomb energy 
difference |AC| computed from Eq. (2.18) with 
4u,=0.590 Mev (this value is derived from the 
energy differences between mirror nuclei). The 
sixth column gives AE,, the nuclear energy 
change corrected for the Coulomb energy differ- 
ence and for the neutron-proton mass difference. 
According to the uniform model, AE, is given 
by Eq. (3.1) and is inversely proportional to A. 
Column 8 therefore contains the value of A-}, 
and Figs. 3.1 to 3.9 contain plots of AE, versus 
A-. A star is attached to energies AE, which are 
considerably (more than about 1 Mev) removed 
from the best straight line on the corresponding 
figure. Some of these deviations can be explained 
by magic numbers (see column 4), and in all 
these cases the deviation is in the right direction, 
i.e., the magic nucleus is unusually stable. We see 
that the magic number 14 is not very potent, but 
the magic numbers 20 and 28 lead to an increase 
in binding energy of between 1 and 2 Mev. There 
are some starred cases in which neither of the 
two nuclei is magic. Some of these deviations 
may be due to experimental errors; the others are 
not explained by either the uniform model or by 
the shell model. The last column contains those 
references which were not included in the Com- 
pilation of Nuclear Data by the U. S. Bureau of 
Standards.” Inclusion of all the relevant refer- 
ences would have lengthened the tables unduly. 


The adventurer is an individualist and an egotist, a truant from obligations. His road is solitary, 
there is no room for company on it. What he does, he does for himself —WILLIAM BOLITHO. 
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Momentum Thrust of a Rocket 


R. K. SHERBURNE AND W. L. WEEKS 
New Mexico College of Agriculture and Mechanic Arts, 


State College, New Mexico 

HE common methods for deriving an expression for 
the thrust developed by a rocket motor vary in 
complexity from simply writing down the result! to the 
derivation of an auxiliary general principle of mechanics.? 
Some writers utilize a theorem from hydrodynamics*‘ or 
an integration of the forces exerted over the exterior and 
interior of the rocket,’ including in the integration the 
reaction force due to the burned propellant leaving the 
rocket. The authors feel that it is desirable to derive the 
momentum thrust (that thrust due to the acceleration of 
the jet mass) directly from classical mechanics and to 
consider any other forces, such as pressure differences over 
the rocket surface, as external forces. This opinion does 
not imply that the derivations referred to above are in- 
correct, but rather that they are open to one or more of 

the following objections: 

(1) They require special study of an auxiliary theorem, 
such as the hydrodynamic theorem or the auxiliary prin- 
ciple of mechanics. 

(2) Some of the assumptions may be open to question; 
for example, the application of the hydrodynamic theorem 
assumes that the flow of the surrounding medium is not 
disturbed and that the jet mass does not mix with the 
surrounding medium. 

(3) They do not clearly present the mechanics of the 
problem; this is particularly evident if the result is written 
down immediately with little or no attempt at a deriva- 
tion. To be more explicit, consider the result as obtained 
in reference 1. It is concluded as being almost immediately 
obvious that the force on a rocket due to the jet is 


F=—mv,., 


where m is the rate at which the mass is leaving the rocket 
and v, is the velocity of the jet relative to the rocket. 
However, if an unwary student applies Newton’s second 
law to the rocket to obtain 


F=—2v.=d(Mv)/dt, 


where M is the mass of the rocket and v is its velocity, 
and, sensibly, attempts to account for the variation of 
rocket mass in the differentiation, he will obtain an incor- 
rect expression. 

(4) Velocities are sometimes determined with respect 
to more than one coordinate system; in particular, the 
rocket velocity v is determined relative to a fixed coordi- 
nate system, while the jet velocity v, is determined relative 
to the rocket, with no discussion presented. 

(5) The relative jet velocity is considered constant in 
most derivations, and some do not include a variation in 
absolute jet velocity as the rocket moves. 

It is true that some of the above methods have particular 
advantages. For example, the method which uses the 
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hydrodynamic theorem includes in the result the pressure 
thrust as well as the momentum thrust. On the other 
hand, the derivation which follows has the advantages that 
(1) it is direct; (2) it applies to the general case in that it 
considers variations in all quantities which can change; 
(3) it attempts to present a clear physical picture. 

We shall apply a principle from mechanics which states 
that the time rate of change of the total momentum of a 
system of particles is equal to the vector sum of the 
external forces on the system. It should be emphasized 
that this principle pertains to a system which always 
contains the same particles. We apply the theorem to the 
system which is defined to be that group of particles which 
comprise the total mass of the rocket and fuel at a time ¢. 
Letting G represent the momentum of this system, 


dG/dt=2F, 


in which 2F is the sum of the external forces. 

We define the following quantities: 1/=total mass of 
rocket at time ¢; v=velocity of rocket at time ¢, measured 
relative to a fixed coordinate system; v;=velocity of jet 
as it crosses the exit plane of the rocket nozzle at time ?, 
measured relative to a fixed coordinate system; Vv, = velocity 
of jet as it crosses the exit plane of the nozzle at time ¢, 
measured relative to the rocket nozzle; and Am=propellant 
mass which leaves the rocket during a time interval Af; 

At the time ¢, the total mass of the system is moving 
with the speed of the rocket, so 


G(t) = Mv. 
At the time ¢+At, a quantity of mass Am has left the 
rocket at an average velocity v;+ 3Av;, so that 
G(t+ At) = (M—Am)(v+Av)+Am(v;+ $4v;). 
By definition of the derivative, 
dG/dt= lim [G(t+At)—G(¢)]/At 
Ato 


= lim [((Mv+ MAv—vAm—AmAv-+Amv; 
Ato 


+ 4AmAv;) — Mv ]/At. 


It is convenient to express the final result for, momentum 
thrust in terms of v., the jet velocity measured relative 
to the rocket nozzle, since this quantity is usually known 
to the designer of rockets. Toward this end, since, vec- 
torially, 

V;=V+V. 
and 
Av;=Av+Av., 


it follows by a straightforward substitution that 
dG /dt= lim (MAv— 4AmAv+Amv,+ tAmAv,)/At. 
Ato 
Note that some of the terms cancel before proceeding to 
the limit, and the remaining terms can be interpreted as 


follows: The first term claims that the rocket has suffered 
a change in momentum Mav, while the second term 


ces” 
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counterclaims that not quite all of M was accelerated by 
the amount Av since some mass had left the system. The 
third term describes the momentum change of the jet 
mass, while the last term corrects this for a possible 
change in relative jet velocity. Proceeding to the limit, 


dG/dt = Mdv/dt+dm/dtv., 


in which the second and fourth terms have been discarded 
as terms which approach zero in the limit. 
The final result is, as is usually obtained; 


Mdv/dt = —dm/dtv.+2F, 


where dG/dt has been replaced by ZF. The minus sign 
indicates that the rocket acceleration is opposite to the 
direction of v.. The pressure thrust, A.(P3:—P.), which 
may exist because of a difference between the pressure P, 
of the jet at the nozzle exit, and the pressure Ps; of the 
surrounding medium, is included in ZF. v, is the actual 
relative jet velocity. A theoretical expression for v. can 
be obtained in terms of the conditions in the rocket motor 
combustion chamber and the pressure at the nozzle exit.® 

The expression is valid for all time, since no restriction 
was placed upon ¢. It should be borne in mind, however, 
that the total mass M is a function of time, so that the 
acceleration of the rocket body is not constant, even if m, 
ve, and =F are constant. 

1 Seifert, Mills, and Sommerfield, Am. J. Phys. 15, 1 (1947). 

2 Rosser, Newton, and Gross, Mathematical Theory of _— Flight 
(McGraw-Hill Book Company, Inc., New York, 1947), p. 

3M. J. Zucrow, Principles of Jet Propulsion and Gas Pusbines (John 
Wiley and Sons, Inc., New York, 1948), pp. 38, 491. 

4G. P. Sutton, Rocket Propulsion Elements (John Wiley and Sons, 
Inc., New York, 1949), . ~ 


5 Reference 3, pp. 38, 
® Reference 4, p. 45, Ee Xs 14). 


A Method for Deriving Various Formulas in 
Electrostatics and Electromagnetism Using 
Lagrange’s Trigonometric 
Identities 


Eppiz Ortiz Mufiz 
College of Agriculture and Mechanic Arts, Mayaguez, Puerto Rico 


N many textbooks in college physics authors very often 
do not derive formulas in electrostatics dealing with 
continuous distribution of charges and in electromagnetism 
dealing with the fields set up by electric currents because 
a summation process, integral calculus, is required. With 
calculus such derivations are easily carried out but without 
it they are rather difficult. 
The purpose of this paper is to present a method for 
computing various formulas in electrostatics and eens 


magnetism using the trigonometric summations 2 cosn@ 
1 


N 
and 2 sinv@ known as Lagrange’s trigonometric identities. 
1 


Throughout this paper all angles will be divided into a 
large finite number of equal units. As the size of the unit @ 
is arbitrary, as small as one pleases, and hence the number 
of units very large, the trigonometric summations can be 
approximated as shown below: 
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Fic. 1. Calculating the work required to bring a charge g from a to b 
in the field of a point charge Q. 


Summation Exact value 


Approximate value 


N : ’ 
2 cosné 1, i i _1 sinNe 
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2 2 sin}0 2 0 
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1—cosN@ 
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2 sin 30 6 
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Ne Ne M 


2 cosn@ = ~ cosn@— ~ cosn@ 
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Nez Nez M 
2 sinnd= =2 sinn@— 2 sinné, 
M 

the above equations can be approximated to 


Ni sin N20 —sinN\6 
> spelt ae 
M 6 


sinn@ = So — (2) 
M 6 

The application of the method will be illustrated by the 
following two examples. 

1. To calculate the work required to move a positive 
charge g from a point a to a point b in the region of a 
positive point charge Q. 

The force acting on g at a distance R from Q is 


=qQ/R?. 


The work done against the electrical force F in a displace- 
ment Al is 


’ (1) 


AW= — F cosdAl= —qQ cos¢Al/R?. 
It can be seen from Fig. 1 that 


Al=AS/sing = R@/cosné 
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Fic. 2. Calculating the magnetic field of a piece of wire carrying a 
9 I. The magnetic field is toward the reader and is represented 
y O. 


and 
R=d/cosné. 
Hence, the work becomes 
W=-—qQ0 sinnd/d. 


The total work required to move g from a to b is 


N2 
W=-—qQ0/d = sinné. 

Mi 
When Eq. (2) is used, 
Ww 


a qQ(cosB —cose) 
a paeete 
Since 

Ra=cosa/d, R,=cosB/d, 
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Since the quotient W/g is defined as the potential dif- 
ference between a and 3, the potential at a point at a 
distance R from Q relative to a point at infinity is V=Q/R. 

2. To calculate the magnetic field at a point P of a 
piece of straight wire carrying a current J. 

According to Ampére’s theorem, a current element JAl 
(Fig. 2) produces a field AH at P given by 


AH=IAlI sin¢/R?, 


then 


where J is the current, R the distance from the element 
to P, and ¢ the angle between JAl and R. 
It can be seen from the Fig. 2 that 


Al=R6/cosn@, R=d/cosné. 
Hence, the contribution AH due to the element becomes 
AH=I@ cosn6/d. 


Since the contribution of all current elements is in the 
same direction, the total field at P is 


Nz 
H=I6/d = cosné. 


M 
When Eq. (1) is used, 
H=I(sin8—sina) /d. 


If the wire is infinitely long, B=2/2 and a=—7/2, then 
H=2I/d. This formula is generally called the Biot-Savart 
law. 

It has thus been shown that for various field problems 
in electrostatics and electromagnetism it is possible to 
derive formulas by means of two trigonometric identities. 
It is advantageous to make derivations by a method 
closely paralleling the calculus approach, since it will mesh 
better with any work the student may take later. The 
derivations lend themselves to be used in a course of 
college physics, and in the light of the author’s experience 
it may prove to be a valuable teaching aid. 
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LETTERS TO THE EDITOR 


Complementarity in the Life Sciences 


R. POLLARD’S interesting article, “Significance 
of Complementarity in the Life Sciences,’! has 
prompted me to make a few remarks. 

I think we must be warned that verbal-interpretations 
such as “‘complementarity” and “indeterminacy” which 
are made directly from an abstract mathematics, because 
the phenomena which this mathematics describe are un- 
tenable to five-sense experience, really tell us more about 
ourselves than about the nature of the phenomena. Because 


words are rooted in five-sense experience they may not be 
appropriate or meaningful for phenomena beyond our 
direct perception, and any attempted interpretations may 
reveal certain predispositions in our own natures. In some 
respects the mathematical equation becomes a projective 
device, not unlike the psychologist’s Rorschach test 
(an abstract series of ink forms) from which verbal inter- 
pretations give insight into personality. 

Though verbal interpretation from abstract mathe- 
matics helps to satisfy our psyches more than the mathe- 
matics can (because the language of the psyche is and can 
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only be the verbal language), there is no necessity that the 
essential intelligence of the mathematics is verbally trans- 
latable. After all, we use mathematics precisely because of 
the fact that words are inevitably handcuffed to certain 
unconscious associations while mathematical symbols 
are not and hence are free to transcend the limitations of 
whatever predispositions are inherent in human nature. 

It should not be surprising to find that in other fields 
(history, biology, psychology) the organization of verbal 
statements divides along a cleavage similar to that found in 
physics. This merely says that the physicist, historian, 
biologist, psychologist, is the same kind of being. If it is 
true that the verbal concepts, ‘complementarity’? and 
“indeterminacy,”’ are applicable to these other fields, this 
informs us of something interesting about the nature of the 
mind, and it is possible that some work may be short- 
cutted in advancing the other fields by properly orienting 
the biologist, historian, etc. But this must suppose, of 
course, that the particular phrasing of ‘‘complementarity”’ 
and ‘‘indeterminacy”’ is general enough to apply elsewhere. 


E.tnu FEIN 
Fisk University 


Nashville, Tennessee 


1W. G. Pollard, Am. J. Phys. 20, 281 (1952). 


Complementarity in the Life Sciences 


ROFESSOR Fein’s comments arise out of a point of 
view which looks upon notions of ‘‘complementarity” 
and “indeterminacy” in physics as being abstract mathe- 
matical concepts. From this viewpoint it then appears that 
what I am suggesting in my paper is just another way to 
put the life sciences on a rigorous mathematical basis, 
or to express their content within the framework of a 
mathematical theory. This, however, was not at all the 
objective I had in mind. Actually, I think it is erroneous to 
regard either complementarity or indeterminacy in physics 
as interpretations of phenomena ‘‘made directly from an 
abstract mathematics’”’ as Professor Fein suggests. Both 
of these concepts have to do with the way in which our 
experience of physical phenomena must be handled if 
we are to achieve an adequate synthesis of it. They are 
most naturally expressed in experiential, not abstract, 
terms. They are statements about the way in which the 
physical world, as apprehended in human perception, is 
constituted. It is true that these concepts, once they were 
clearly recognized, were, in physics, expressed in abstract 
mathematical terms, e.g., the commutation relations, and 
used in this form to obtain the quantum mechanics. But 
if my suggestion that similar concepts might prove fruitful 
in the life sciences should be borne out, I would not expect 
that they would prove to be analogously expressible in 
mathematical terms, but rather that they would be applied 
in ways which correspond with the normal means of 
scientific investigation which already characterize the 
several fields involved. 
A more profound question raised by Professor Fein 
concerns the extent to which any given concept in science 
is characteristic of ourselves as perceiving subjects, and 
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the extent to which it is characteristic of the nature of 
phenomena apart from our mode of perception. This 
question, of course, goes well beyond considerations of 
complementarity and indeterminacy and ultimately applies 
to all knowledge, scientific or otherwise. The point of view 
from which my paper is written is a purely pragmatic one 
and is, as Bohr expresses it in one of the passages I quoted, 
“only an endeavor to clarify the conditions, in each field 
of knowledge, for the analysis and synthesis of experience.”’ 


WILLIAM G. POLLARD 
Oak Ridge Institute of Nuclear Studies 
Oak Ridge, Tennessee 


Negative Coefficient of Expansion of 
Stretched Rubber 


HILE picking up our children’s toys from the 

livingroom the other day, I threw an old, half- 
deflated balloon in the fireplace to watch it expand and 
burst. I reasoned that the side nearest the fire would be- 
come soft, expand, and break, but, for the moment, I had 
forgotten that the thermal coefficient of expansion of 
stretched rubber is negative. To my surprise, the side 
nearest the fire did not become soft enough to expand 
but violently contracted instead. This I could observe 
because the balloon was of the multicolored (marblized) 
variety. 


Linfield College 
McMinnville, Oregon 


GEORGE BARNES 


An International Language for Science? 


N the September number of this Journal, Forrest F. 

Cleveland discussed under the heading ‘‘Book Reviews”’ 
the necessity of a neutral and easy international language 
for publication of scientific results in general, and the 
merits of a certain linguistic project ‘Interlingua’ in 
particular.! I heartily agree with what F. F. Cleveland 
writes about the advantages of the introduction of such 
international language in general. Many of you who have 
attended international conferences in which French or some 
other foreign language was the official medium know how 
much simpler it would be to follow the scientific discus- 
sions if the language difficulty were less formidable. 

I may add another important argument. In the past 
centuries Latin has been used as the international language 
of science, later French, not so long ago German, and quite 
recently (since about 1933) English. In fact, since Latin, 
because of its difficulty, lost its position as a neutral means 
of communication, the national language of the country 
producing the greatest number of important scientific 
publications has always served as the international language 
of science. This, however, is not always the same country, 
and we therefore cannot expect English to remain the 
international language of science for all centuries to come. 
Rather than wait till we all have to learn Japanese, 
Spanish, Russian, Danish, Chinese, or who knows what 
language next, we should anticipate this change of inter- 
national language by returning to the idea of a neutral 
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language, avoiding the irregularities of grammer which 
made the now obsolete Latin inappropriate. Acceptance 
of a new easy and neutral language for science, however, 
can be achieved only gradually and only if enough scientists 


will join in promoting its use by publishing from time to 


time results in such language themselves. This is what I 
personally—like several others—have been doing with 
some of my publications, not only to set an example but 
also to show that it can be done.? 

On the other hand, I heartily disagree with the proposal 
that we use “Interlingua,’’ or any of the other hundreds 
of language projects that are not real languages. A language 
that is not spoken and thought by a group of people is not 
a living language, and nobody could readily express him- 
self in a ‘paper’ language during a debate in an inter- 
national conference. Also, it is practically impossible to 
interest large numbers of scientists in spending their time 
learning a language for which they have no use and in 
which not even some literature exists. Most of these so- 
called international ‘‘languages’’ have never grown nor 
do they show their ability to grow beyond the scope of a 
project on paper. A couple of them may have been learned 
by a few people, but never to the extent of becoming 
living languages with character of their own. 

Exceptions to this rule up to the present have been only 
two. Volapiik, when it first was published, caught the 
fancy of millions and for a short time was learned by many. 
However, its grammar was far too complicated, and it 
broke down completely when the first international con- 
gresses of Volapiikists were held and the pronunciation 
turned out to be too hard because of its eight different 
vowels. 

The second exception is Esperanto. The start of Espe- 
ranto in the wake of the failure of Volapiik was slow. 
Progress, however, was rapid from the time the first 
Esperanto congresses were held and fluent conversation 
in Esperanto was shown to be an easy matter.’ Esperanto, 
contrary to Volapiik, did not die a natural death but 
developed to become a spoken language. The author of 
this Letter has often spoken it himself with foreigners 
from many countries. During such conversations—or 
when writing a letter in Esperanto—one need not arti- 
ficially translate into Esperanto, but one can readily think 
directly in this language. Esperanto, therefore, cannot 
be put on one line with Occidental, Ido, Interlingua, or 
whatever those other projects may be called, for Esperanto, 
contrary to these projects, has become a real and live 
language. 

Remarks that this or that word, rule of grammar or 
of spelling in Esperanto might better have been chosen 
differently sound just as useless and a lot less convincing 
than similar remarks on the flaws of English, whose spelling 
and pronunciation are scarcely less than farcical. It is 
quite possible that some Interlingua people think that 
“lingua” looks neater (to them personally, that is!) than 
the Esperanto word “‘lingvo,’”’ but the introduction of an 
international language is not helped by trying to promote 
one’s own taste and making new projects (of which Inter- 
lingua is one of the newest). It is only helped by a joint 
endeavor of all of us to push ahead the one existing well- 
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established neutral international language, which, by the 
way, is as free from flaws as one can reasonably expect of 
any. scientifically man-made language—that is, a lot freer 
from flaws than man-made so-called ‘‘natural’’ languages. 
Moreover, it has a well-developed world literature, of 
which I personally have over a hundred volumes on my 
shelves—a fraction of the total. This is but one of the 
reasons why knowledge of Esperanto can be worth while 
right now. It also offers possibilities of much wider inter- 
national contacts than can easily be obtained otherwise. 
If these lines can convince anybody that our personal 
studying of Esperanto is at least as worth while and im- 
portant for the future of the world as campaigning for 
control of atomic energy and the like (and you'll find that 
it is much more fun, too), then I’ll gladly answer personal 
requests of information on Esperanto courses, diction- 
aries, textbooks, literature, magazines, journals, and 
societies. 
FREDERIK J. BELINFANTE 


Purdue University 
Lafayette, Indiana 


1F. F. Cleveland, Am. J. Phys. 20, 382 (1952). 

2 Text in Esperanto: Prog. Theoret. Phys. 3, 460 (1948); 4, 165 
(1949); 6, 202 (1951). Also, Scienca Revuo 1, 98, 100, 106 (1949); 
3, 61, 67 (1951). 

Abstract in Esperanto over English text: Physica 6, 849, 870, 887 
(1939); 7, 305 (1940); 12, 1, 17 (1946). 

3 There have been severe set-backs during each of the world wars and 
because of persecution of Esperantists by Hitler wherever he was in 
power, and (since about 1936) also in Russia. Recently, there has been 
quite a comeback of Esperanto in many countries. 


Further Demonstrations in Analytical Mechanics 


PURSUE in this note my ambition to demonstrate cer- 

tain of the classical problems of analytical mechanics, 
and I report here two which lend themselves elegantly 
to very close approximation. Indeed, these demonstrations 
possess a sensible beauty which can be discovered only by 
seeing them actually work. If the charge against. them 
(as demonstrations) is that they are simple, I must confess 
that I am always intrigued by the elegant simplicity of 
Nature. 

Problem 1: Find a point P on a vertical circle such that 
the time of descent along a radius to the center is the same 
as along a chord to the lowermost point of the circle. It 
is easily shown that the point P lies 60° from the top of 
the circle. : 

Demonstration: On a sheet of plywood 4 ft by 8 ft, 
standing on end, describe a semicircle with the mid-point 
of a long edge as center O. Draw a vertical diameter; call 
its lower end A. With a protractor locate point P. At 
points O, P, and A set wood screws or screw-eyes and at- 
tach piano wire along PO and PA, the wire carrying small 
riders. (A double electric lead connector makes a very good 
rider, being a metal cylinder about 2 cm long with a smooth 
hole along its length.) The wire can be drawn up taut by 
turning the screws or screw-eyes as turnbuckles. This 
arrangement gives a path of 4 ft along a radius and about 
7 ft along a chord. On simultaneous release of both riders 
at P, it is observed that both arrive simultaneously at O 
and A. Indeed, with one’s back on the experiment the 
simultaneity of impact is easily detected. Since the time of 
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descent is really too short for safe clocking, the simultaneity 
of descent as detected by the eye or ear is safe evidence. 

Problem 2: Show that the time taken by a particle to 
slide down any chord which begins at the highest point 
of a vertical circle is constant. It is easily shown that the 
time is given by 2(r/g)4, where r is the radius of the circle. 

Demonstration: The board used in Problem 1 can be 
fitted with screw-eyes and chords bearing riders. Too many 
at once, however, crowds the riders at the highest point. 
The simultaneous arrival at the circumference is an 
engaging thing to witness. 

JuLrus SUMNER MILLER 


Dillard University 
New Orleans, Louisiana 


Some Remarks on D. Bernoulli’s Formula 


HE paper of G. A. Lindsay on ‘‘Pressure Energy and 
Bernoulli’s Principle’! invites a commentary. First, 

we may observe that the name of Bernoulli’s principle 
is usually attributed to the very important theorem desig- 
nated otherwise as the principle of virtual displacement or 
work. This was first proposed intuitively by John Bernoulli 
(1667-1748) and afterwards elaborated with more rigor 
by Lagrange in his system of mechanics. But what G. A. 
Lindsay’s paper is concerned with is the remarkable 
formula demonstrated by Daniel Bernoulli (1700-1782) ,? 
the paramount importance of which has been realized 
only in our time. This formula can be obtained from two 
starting points which, for most physicists, may seem very 
far apart: (1) the principle of conservation of energy, and 
(2) Newton’s second law.* It would be, of course, faulty 
to consider the pressure as an energy. It is but a factor, 
i.e., the intensity, the volume being the extensity, of a form 
of energy called volume or spatial energy. But the demon- 
stration of D. Bernoulli’s formula from the principle of 
conservation of energy can be done quite correctly and has 
so many merits that it should be included into syllabuses. 
The point that it can be also deduced starting from 
Newton’s second law is very interesting and can be most 
easily explained on the lines of my system of energetics. 


Z. KLEMENSIEWICZ 
40 Foxgrove Road 


Beckenham, Kent, England 


1Am. J. Phys. 20, 86 (1952). 

2 At least three members of the Bernoulli family gave distinguished 
services to physics and mathematics. 

3 See, for example, E. Schmidt, [Thermodynamics, (Oxford University 
Press, London 1949), English translation, pp. 324, 332], who was rather 
startled by this fact. 


On an Aspect of Demonstration Experiments 


HE place which demonstration experiments occupy 

in physics instruction needs no discussion. As in- 
structional devices they are probably unexcelled. In the 
hands of a good demonstrator, accompanied by the ap- 
propriate “language,’’ even the most elementary experi- 
ment can be stirring. There are, of course, as many ways to 
demonstrate an experiment and as many ways to ‘“‘say 
the physics” as there are teachers, but one very important 
aspect exists of which, I have repeatedly observed, teachers 
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fail to take advantage. This aspect I might title ‘“‘Ques- 
tions to be Asked.’’ During the demonstration the teacher 
usually recites what he has, what he is going to do, and 
what will be observed to happen. And it usually happens. 
(If it does not, some good physics can still be taught.) 
Now the detail to which I refer explicitly is this: it is highly 
instructional to ask here and there, ‘‘Will this happen?’’, 
or, “Why will this not happen?” I can best illustrate the 
point by reference to an elementary demonstration. 

Consider a demonstration of linear momentum. I use 
a track made of angle-iron on which reside a string of 
billiard balls. One ball is taken up the incline and released. 
After impact, one ball runs off from those at rest. Two 
balls are taken up the incline; after impact, two balls run 
off from those at rest. Suppose now that more are taken up 
the incline than remain on the level track. What will 
happen? How many will go away? It is surprising to dis- 
cover how puzzling this question is. But let’s go on. After 
it is shown that as many go away as came down, I like 
to ask this: with two coming down, say, we found two 
going away after collision. The mv is taken care of. Now 
could not the mv be accounted for by one ball going away 
with a velocity 2v? And with this I repeat the experiment 
several times, but this does not happen! Why? (I once 
had a good student who had read about entropy. Regard- 
ing the billiard balls, he was of the opinion that it might 
eventually happen, just as the water in a glass might 
suddenly freeze on a hot day, or the original order in a 
shuffled deck of cards be gotten again.) 


JuLtius SUMNER MILLER 
Ford Foundation Fellow 


University of California 
Los Angeles, California 


Magnitude of the Newton 


UR teaching methods should utilize as many of the 

five senses as is practicable. Whenever we make use 

of the sense of touch, we add to the ability of the student to 

grasp the idea under consideration. To present the student 

with an appreciation of its magnitude, we have constructed 
a demonstration ‘‘newton.”’ 

The newton is defined as the force required to give 1 
kilogram an acceleration of 1 meter/sec?. The acceleration 
is approximately one-tenth that of gravity, so the earth 
pulls 0.1 kilogram with a force of approximately 1 newton. 
The error is 2 percent. The precise mass required for an 
earth-pull of one newton is 1/g where g is the local value 
of the acceleration of gravity. This mass is more nearly 
102 grams. 

To make our demonstration ‘‘newton,’’ we added the 
appropriate number of grams of solder to a brass 100- 
gram weight so as to obliterate the original 100 grams. 
Next we stamped on the solder the phrase ‘“‘1 newton at 
U. of Ala.’’ When we pass around our “‘newton”’ in elemen- 
tary classes, it becomes a “‘real’’ physical quantity to the 
student, not just another theoretical quantity to com- 


plicate life. 


Cart C. SARTAIN 
University of Alabama 
University, Alabama 
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Theory of Superconductivity. M. von Lave, translated 
by L. MEYER AND W. Bann. Pp. 140+x, Figs. 36, 
154X234 cm. Academic Press, New York, 1952. 
Price $4.00. 


The theory of which M. von Laue, Nobel laureate in 
physics, has so carefully written is the phenomenological 
electrodynamics which results when the equations of 
Maxwell are applied to the empirical equation for a super- 
conductor by Fritz and Heinz London. The author has 
collected a number of interesting problems on variously 
shaped superconductors about which the field distribution, 
current distribution, and related effects are worked out. 
The experimental physicst will find this 140-page book 
a valuable addition to his reference collection. The transla- 
tion was completed in 1949, but there are references and 
remarks about work published in 1950. The theory and 
mathematical steps are well written and easy to follow. 
Professor von Laue considered the quantum-theoretical 
efforts made on the detailed picture of superconductivity 
as premature and has confined his writing to only the 
classical theory. The ponderous tensor notation which 
the author presents to the reader in the early pages is 
avoided in most of the application to specific problems. 

The praise which this book received from Fritz Sauter 
[see Naturwissenschaften 37, 263 (1950) ] that the phenom- 
enological theory is ‘‘in ersten Linie mit dem Namen von 
Laues verbunden” is not wholly justified. In part this 
impression might be given because von Laue does not 
burden himself with giving prior credit. Thus, for example, 
the experiments on superconductors at very high frequency 
had been carried out by H. London, who also first gave the 
proper theory for the energy absorption. But in the book 
the inaccurate studies by McLennan and co-workers are 
evaluated with a theory claimed by von Laue. The last 
chapter of the monograph, “A Nonlinear Extension of the 
Theory,’”’ deals with some highly controversial theory 
which was stimulated by Heisenberg’s theory of super- 
conductivity. 

Chapter 19 deals with the phase transition into the 
superconducting state and discusses the intermediate state. 
No use is made of the concept of a surface energy between 
the superconducting and normal state. The interested 
reader will welcome a supplement to this chapter by 
reading von Laue’s article in Amnalen der Physik 10, 
296-316 (1952). The editors of the journal were cooperative 
in requesting Professor F. London to add some helpful 
comments in response to von Laue’s paper and in this way 
agreement on the surface energy has been reached. The 
unconvincing work by Stark, Steiner, and Schoeneck 
received some attention and yet such things as the classical 
work by Kikoin and Guber on the gyromagnetic ratio of 
a superconductor is not included in the monograph. 


C. F. Squire 
Rice Institute 


Proceedings of the London Conference on Optical Instru- 
ments 1950. Pp. 264+-xv, Figs. 103. John Wiley and 
Sons, Inc., New York, 1952. Price $7.00. 


The present volume is a compendium of papers resulting 
from a topical symposium on modern optical instruments 
held under the auspices of the International Optical 
Commission in London during July, 1950. This conference 
on optical instruments was designed to follow in the 
tradition of the great International Optical Conventions 
of 1905, 1912, and 1926. The objective of the conference 
was principally the design and application of modern optical 
instruments such as photographic objectives, reflecting 
microscopes, grating spectrometers, phase contrast micro- 
scopes, spectrophotometers, reflecting telescopes, and 
other modern developments. 

The introductory remarks by Sir T. Merton aptly 
point out that optics today is in a period of extended 
maximum activity and is actually becoming chiefly an 
applied science. The tremendous accumulation of instru- 
mental research during the past decade makes the need 
for a volume such as this one imperative in order to bring 
under one cover at a reasonable cost a condensation and 
integration of those developments which are often pub- 
lished in such obscure literature as patent descriptions or 
are otherwise scattered throughout the scientific literature 
of the world. 

Three experts contribute to the chapters on ‘‘Photo- 
graphic and Projection Lenses.” R. Kingslake of the 
United States writes comprehensively on ‘‘Some Recent 
Developments in Photographic Objectives.” Excellent 
descriptions are accompanied by numerous illustrations 
of various lens types. Conventional designs of double, 
triplet, four-element, Petzval-type lenses as well as wide- 
angle lenses, telephoto, and zoom lenses are discussed. 
Reference is made to the recent application to the photo- 
raphy of catadioptic systems involving a reflection element 
and a correction plate. H. H. Hopkins of Great Britain 
relates in detail the theory and properties of the zoom 
lenses in the section entitled ‘‘A Class of Symmetrical 
Systems of Variable Power.” In the chapter on “A New 
High Aperture Photographic Objective Having a Spherical 
Field,” A. Warmisham (Great Britain) points out the 
great flexibility and precision attainable with the Schmidt- 
optics type of lens. 

A truly recent development in optics is the reflecting 
microscope. A purely reflecting system is completely 
achromatic, and only the reflectivity of the surfaces limits 
the useful wavelength range. A Burch microscope has 
been used from 240 my in the uv to 14, in the infrared. 
These three chapters are ably contributed by R. Barer 
(Oxford), A. Bouwers (Netherlands), and D. S. Grey 
(United States), and emphasize so well how the biologist 
and the physicist can profitably cooperate in the field of 
biophysics. 

It is always a thrill in science when one witnesses the 
invention of an instrument which will improve the order 
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of measurement by a full magnitude. Such is the case for 
the echelle grating discussed by G. R. Harrison and C. L. 
Bausch of the U. S. A. in the first chapter under “Gratings 
and Grating Instruments.’’ ‘“‘A Grating Spectrograph 
with Interchangeable Receivers for Industrial Purposes’ 
is discussed by E. Ingelstam (Sweden), and ‘‘Some New 
Mountings of Diffraction Gratings and the Plane Mirror” 
by E. Hulthen and E. Lind of Stockholm. 

F. W. Cuckow of Great Britain gives an excellent de- 
scription of the advantages of the phase-contrast micro- 
scopic method applied to the study of surface structures 
in metallurgy. The advantages are clearly shown by com- 
paring an ordinary photomicrograph with an electron 
micrograph and with a phase-contrast photomicrograph 
of alloy steel, each at the same magnification. The de- 
scription of the apparatus and of the techniques necessary 
is very well illustrated. M. Francon of France presents the 
recent progress in phase-contrast microscopy of the Insti- 
tute of Optics in Paris. He describes a mirror reflection 
scheme which is absolutely stigmatic and aplanatic to be 
placed between the sample and the objective of the micro- 
scope to give an image in phase contrast without the 
use of a phase plate. 

A comprehensive review of spectrophotometry in the 
United States is presented by S. S. Ballard which would 
serve as an excellent guide to anyone attempting to reach 
a decision on the type of instrument best suited for a 
particular application. The photometric equipment of the 
Institute of Optics in Paris is thoroughly described by 
F. Desvignes. An excellent discussion of the slit-width 
dependence of the output is described in French and pre- 
sented graphically. 

A detailed discussion of the modern trends in reflecting 
telescope design is presented by E. H. Linfoot of Great 
Britain and recent Swiss developments by W. Lotmar. 
Nearly every modern type of instrument will be found 
described in these two chapters including the Baker- 
Schmidt type for the ADH-observing station at Blom- 
fontein and a similar instrument for the St. Andrews 
Observatory. 

Erik Bergstrand of the Swedish Geographical Survey 
reports in Chapter 17 on an apparatus intended for meas- 
urements of distance by high frequency light signaling. 
It is based on the principle used by Fizeau for terrestrial 
determinations of the velocity of light modified by making 
use of a Kerr cell and a photocell. Conversely, if the dis- 
tance is known, the velocity of light may be determined. 
The method and results are thoroughly described and 
illustrated. 

Chapter 18 is concerned with apparatus for the pho- 
tometry of optical instruments by J. Guild of Great Britain. 
He describes apparatus designed at the National Physical 
Laboratory for measuring such photometric properties 
of telescopes and photographic lenses as the paraxial 
transmission factor, the relative brightness of various 
zones of the field, and scattering and multiple reflection 
effects. 

The concluding chapters are written by S. S. Ballard 
of the United States on ‘Synthetic Optical Crystals,” 
by I. C. Gardner of the National Bureau of Standards on 
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“‘New Types of Optical Glasses Available in the United 
States,” and by H. C. Raine of Great Britain on ‘‘Plastic 
Glasses.’ These three chapters rather thoroughly sum- 
marize the developments during the last decade on syn- 
thetic optical materials. It will be noticed that Ballard’s 
article does not contain a bibliography, but he points out 
that the firm of Baird Associates, Inc., will shortly publish 
a handbook of the properties of synthetic crystals in the 
Journal of the Optical Society of America. At this writing 
(July, 1952) this handbook has not appeared in that 
journal. 

The book is exceedingly well written throughout and 
profusely illustrated. Each chapter contains an extensive 
bibliography which will be a great time-saver for the 
researcher. An unusual feature in the book is the presence 
of a blank page for notes at the end of a couple of the 
chapters. This is an excellent reference book for the teacher 
of optics whether teaching at the elementary or at the 
advanced level. The book reveals’ the tremendously rapid 
progress being made in optical instrument design today 
and should be inspiring reading to the student of optics 
as well as a challenge for him to contribute to these ex- 
citing events. Perhaps another International Optical 
Conference may be held in a few years to record these 
events in a treatise as well done as this one. 


Davip M. GaTEs 
University of Denver 


Low Temperature Physics. F. E. Simon, N. Kurt, J. F. 
ALLEN, AND K. MENDELSSOHN. Pp. 132, Figs. 59. 
Academic Press, Inc., New York, 1952. Price $3.50. 


Ten years ago there were fewer than a dozen labora- 
tories in the world equipped for research in low temperature 
physics—today there are more than fifty in the United 
States alone. Such institutes represent fairly formidable 
investments in time and money; why then has such an 
interest, almost unprecedented, arisen in so short a time? 
This small book, aimed principally at the nonspecialist, 
provides much of the answer. 

The four chapters, of roughly equal length, represent a 
series of lectures, one by each author, delivered before the 
Royal Institution in February and March, 1950. The 
writers are distinguished British physicists, each of whom 
has made substantial contributions to cryogenic physics, 
and no more authoritative group in this field could easily 
be assembled. 

The first chapter, by Simon, is a rather general survey 
(a sort of introduction to what follows) covering, in simple 
nonmathematical language, a diversity of low temperature 
phenomena. He brings out very clearly the predominant 
role of ‘‘macroscopic quantum effects’ in the phenomena 
occurring at liquid helium temperatures. A number of 
enticing suggestions for further experimentation in little- 
exploited branches of the subject are included. 

Chapter II, by Kurti, deals with ultra low tempera- 
tures—the technique of adiabatic demagnetization. This 
is a very readable account, again with a minimum of mathe- 
matics, of the principles and practices involved in the 
production and use of temperatures within 0.01 degree, 
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and less, of absolute zero. A particularly clear exposition 
is provided of the proposals involving production of 
temperatures around 10° °K by two-stage demagnetiza- 
tion, using electronic and nuclear spins, a technique which 
is currently under investigation in several laboratories. 

The third chapter, by Allen, discusses the strange 
phenomena associated with the ‘‘quantum-liquid” helium 
II. This, again, is a very clear, commendable, and critical 
account of most of the discoveries up to 1950. The author 
suffers the fate of all reviewers in a fast-moving subject— 
a lot of the exciting work on He’ was published after the 
lectures had been given. 

The book closes with Mendelssohn’s chapter on super- 
conductivity; again an excellent account considering the 
size of the subject and the space available. Mendelssohn 
closes the chapter with a rather moving appeal for a ‘‘new 
concept in physics’’ to explain the so diverse and yet so 
similar properties of superfluidity and superconductivity. 

Despite the multiplicity of authors, the book has a 
surprisingly coherent style. And the formal text is il- 
luminated (in the English manner) with many simple and 
intriguing observations. Thus, ‘‘at 0.03°K the vapor pres- 
sure would be so small that in a volume as big as our whole 
galaxy not a single atom would be in equilibrium with 
liquid helium’’ (Simon), or ‘‘the heat capacity of 10 gm 
of iron ammonium alum at 0.05°K is about the same as 
that of-a quarter-ton of copper . 


. . at the same tempera- 
ture” (Kurti). 
C. T. LANE 


Yale University 
° 


General Education in Science. Edited by I. BERNARD 
COHEN AND FLETCHER G. Watson. Pp. 217+ xviii, 
14X21 cm. Harvard University Press, Cambridge, 
1952. Price $4.00. 


Anyone who has followed the growth of general educa- 
tion in the sciences can hardly fail to be stimulated by a 
reading of this small volume. Fifteen scientists, presenting 
their views before a Workshop in Science in General 
Education, can hardly be expected to concur on many 
aspects of the subject, even if nine of them are from 
Harvard! They are unanimous, however, in believing that 
science has played and will play an important part in our 
civilization and that traditional introductory courses in 
the sciences fail to convey an adequate understanding of 
the nature and role of science to those students who will 
major in some nonscience field. This indictment of tradi- 
tional science courses is laid by some to the increased 
specialization in the sciences during the past half-century 
and the consequent narrowing of the training received by 
college teachers. René J. Dubos of the Rockefeller Insti- 
tute for Medical Research decries the widespread equation 
of science with technology by most people, and their 
total ignorance that science can enrich life with new 
esthetic, emotional, and intellectual satisfactions. He traces 
this attitude to the failure of science in the nineteenth 
century, through positivist philosophy, to solve the riddle 
of the universe and existence. 


Several of the contributors give their ideas on what are 
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the most important aspects of science that should be pre- 
sented to nonscience students, and the most effective 
methods. They are generous enough to acknowledge that 
no single approach has all the answers and that is it more 
important for a teacher to use one in which he feels a 
natural interest and some ability to prepare himself. The 
case method, the problem approach, critical thinking in 
science, the emphasis on history, or on philosophy, are all 
treated by teachers who have developed and used them. 
Anyone contemplating the organization of a general educ- 
ation science course should place this volume on his 
personal list of required reading. Worth-while thoughts 
about what such teachers should know about philosophy 
are presented. Three of the contributors make valuable 
suggestions on how a teacher may enlarge his background 
in historical materials, either for developing historical 
cases or merely to convey a sense of history. Even the 
worrisome matter of evaluation is treated helpfully. 

To strive for an understanding of what science is and 
what it does for a sense of its history or a glimpse into 
some of its philosophical implications—apart from mere 
knowledge of some of its facts, principles, and applications 
—is what is being attempted in these various courses for 
nonscience students. But this reviewer cannot help 
wondering about the remaining 31 percent (at Harvard) 
of “‘science’’ students. How do they acquire this liberalizing 
“understanding of science’ which these contributors feel 
is so important. Its absence from the standard intro- 
ductory courses must be presumed by the very develop- 
ment of general education courses. Perhaps it is acquired 
at advanced undergraduate or graduate levels. (The 
narrowness of modern graduate work in science raises 
some doubts on this point.) But very few of the 31 percent 
of freshman carry through to graduate study in science. 
A number follow engineering, others transfer to a non- 
science major before the bachelor’s degree. In short, a 
large block of students is still not receiving some of the 
benefits of general education in science. 

This is a difficult situation to deal with. Experience has 
shown that occasional footnotes or side remarks are of 
little value. A real impression must be made. But training 
in current science and technology is a time-consuming 
task, and, where a gap can be pried open, most people 
feel that a little more of the humanities or the social 
sciences should be wedged in. Two of the contributors to 
this volume have made brief suggestions to meet this 
problem. Paul B. Sears, a biologist at Yale, claims that 
the objectives of both general education and basic prepara- 
tion in the sciences could be met by a unified and reorgan- 
ized two-year sequence for science and nonscience students 
alike. Sidney J. French, a chemist at Colgate, has interest- 
ing ideas for reorgainzation of science courses at the upper 
levels. Both admit, along with this reviewer, that there 
would have to be sacrifice in a number of places. Judged 
by the current product, our science training programs are 
good. Is it too much to expect that, with careful planning, 
they could be improved by including some of the objec- 
tives of general education in science? 

CLEMENT L. HENSHAW 
Colgate University 





148 


Impact of Science on Society. United Nations Educational, 
Scientific and Cultural Organization, Paris, France. 
Quarterly journal. Annual subscription $2.50, single 
copies 75¢ (in U.S.A.). Distributed in various coun- 
tries by UNESCO sales agencies (see below). 


One of the most frustrating features of the life of today’s 
scientist or science teacher is the universally acknowledged 
difficulty of even scanning the material offered by the many 
journals in his own and related fields. It would therefore 
appear to require some courage to bring to the attention 
of the readers of the American Journal of Physics still 
another publication that may legitimately claim their 
notice. However, even apart from its inherent excellence 
and the unique function it was created to serve, there need 
be no hesitancy in recommending to scientists in all fields— 
and indeed to every thoughtful person who cares to under- 
stand the importance of science in the affairs of men— 
this new and valuable quarterly produced by the Depart- 
ment of Natural Sciences of UNESCO. For it is precisely 
the inclusion of such a journal in the busy scientist’s 
regimen of reading that will balance the many technical 
items he must peruse and that is capable of reminding 
him of the place and importance of his special field in the 
life of our time. 

As the name and sponsorship of the publication suggest, 
the purpose of IMPACT is to focus attention throughout 
the world on science as a powerful social force. Much 
excellent work is being done by social scientists in studying 
and appraising the consequences for society of past ac- 
complishments of science and technology. The editors of 
IMPACT feel, however, that the delay in the social study 
until after the consequent processes are already in motion 
is unfortunate. What they hope to accomplish through this 
publication is to stimulate social scientists, educators, 
economists, legislators, and leaders of thought to consider 
the possible future effects of present scientific discoveries. 
It must be granted that in many cases these consequences 
can be foreseen, at least in a general way, at the time when 
new scientific or technical advances are announced. 

These aims are realized through the medium of contri- 
buted articles on subjects concerned in some way with 
the interrelations of science and society; and a number of 
rather complete and detailed summaries of current books 
or articles on this field round out the text of each issue. 

Although IMPACT was founded two years ago, it 
now appears with completely revised content, style, and 
format. The Summer, 1952, issue which is at hand for 
review includes the following articles: ‘“‘The Education of 
the Public in Science,’’ by I. Bernard Cohen; “The Bulletin 
of the Atomic Scientists—A Magazine for Science and 
Public Affairs,” by Eugene Rabinowitch; and “A Basis 
for Social Physics,” by John Q. Stewart. Professor Cohen’s 
paper recognizes the disturbing problem arising from the 
fact that modern means of communication give the average 
person much information about science but contribute 
little to its understanding; and he discusses this question 
in the light of his own experience in teaching science to 
college students. This article is one that every teaching 
physicist will wish to read with care. 


ANNOUNCEMENTS AND NEWS 


Because the Bulletin of the Atomic Scientists and the 
organization which it represents have brought to the 
attention of American scientists their social responsi- 
bilities in a way never before attempted, the editors in- 
vited Dr. Rabinowitch to prepare his description of the 
Bulletin and its aims. This article will be read with great 
interest particularly by foreign scientists who are not 
yet acquainted with the Bulletin or who are unable to 
obtain it regularly. 

Readers of this Journal probably are already familiar 
with the name of the subject discussed in the article of 
Professor Stewart through his paper in the May, 1950 
issue. Sharing the conviction of many scientists, both 
social and physical, that the social sciences will eventually 
be developed toward the status of exact disciplines in- 
volving quantitative laws and generalizations, Dr. Stewart 
has himself succeeded in discovering several of these 
relationships. In the paper in IMPACT he sketches a 
program for working out a ‘“‘physics of society’’ and ex- 
amines some of the highly suggestive comparisons that can 
be drawn between certain physical and social phenomena. 

Subsequent issues of the magazine are to contain article 
on such subjects as ‘‘National Development Based on 
Science: The Case of Denmark,’’ “The Parliamentary 
and Scientific Committee in Britain,’ ‘‘Population Move- 
ments among Animals and People,” ‘‘The Origins of the 
Popularization of Science in France,’”’ and ‘‘Cybernetics 
and Civilization,” written by authorities from various 
countries. To judge from the contents of the current issue 
and the offerings promised in future ones, not the least 
of the uses of JMPACT to the’teacher of science will be 
to furnish material for meaningful classroom discussion, 
particularly in courses for nonscience students. 

IMPACT is published in French as well as in an English 
edition. In the United States, the English edition is avail- 
able through the Columbia University Press, 2960 Broad- 
way, New York City. 

IRA M. FREEMAN 
Rutgers University 


The Theory of Relativity. C. M@LLER. Pp. 386+-xii. Oxford 
University Press, London, 1952. Price $7.00. 


The publication of a book on the theory of relativity 
is a rare event these days. The appearance of the present 
book is to be welcomed, however, not only for this reason, 
but also because it provides a serious, carefully written, 
and unusually thorough account of the theory that can be 
highly recommended for the student of theoretical physics. 
For that matter, it can be recommended for anyone with 
an adequate knowledge of mathematics and physics who 
would like to become acquainted with the beautiful 
structure of the relativity theory and its applications. 

Well over half of the book is devoted to the special 
theory of relativity. The remainder deals with the general 
theory and includes a brief introduction to relativistic 
cosmology. Although the book concerns itself only with 
the classical relativity theory, it is evident that its author 
is interested in the quantum theory and the application 


. of relativity theory to the latter. Thus, the motion of waves 





ANNOUNCEMENTS AND NEWS 


receives considerable discussion, including a comparison 
of the phase velocity of a wave with the velocity of a 
particle. This paves the way for the introduction of the 
ideas of de Broglie concerning particles and waves which 
played such an important role in the development of wave 
mechanics. There is a derivation of the expression for the 
Thomas precession, the precession of a vector associated 
with a particle, such as the spin of an electron, when the 
particle undergoes an acceleration not parallel to its 
velocity. As an example of a field, the scalar meson field 
is briefly treated. 

As he points out in the preface, the author tries to make 
clear the physical meaning of the quantities occurring in 
the relativity theory by giving a three-dimensional formu- 
lation, in addition to the usual four-dimensional one. 
For example, the fundamental or metric tensor which 
determines the geometry (and gravitation) of the four- 
dimensional space-time continuum is split up into a three- 
dimensional metric tensor, a gravitational “scalar poten- 
tial,’”” and a gravitational ‘‘vector potential.” 

The author has performed a valuable service in bringing 
together many of the important results obtained during 
the early part of this century by Einstein, Poincaré, 
Minkowski, and others. Physics work tends to get buried 
under the work of subsequent generations. How many of 
the younger physicists are familiar with the disagreement 
between Abraham and Minkowski concerning the proper 
expression for electromagnetic momentum density in a 
material medium? Mller presents arguments supporting 
Minkowski’s standpoint. These might well be noted by 
future writers of books on electromagnetic theory. 

In the general theory of relativity both geometry and 
gravitation are described by the expression for the interval 
or line element between two neighboring events in terms 
of the coordinates and their differentials. In principle, 
this interval can be measured by a suitable procedure 
involving the use of clocks and measuring rods. In this 
connection, on p. 232, the author writes: “ .. . if all 
measuring-rods in the neighborhood of a given point for 
one or another reason were dilated at the same rate inde- 
pendently of the material of which they are made, it would 
be impossible to observe this dilation and in an unambig- 
uous way to correct for it. Therefore, there is no well- 
defined meaning in the statement that such an expansion 
of the standard rods has taken place, and from a physicist’s 
point of view the metric tensor and the geometry obtained 
by measurements with the natural standard rods must 
be the ‘real’ geometry on the surface.’’ One wonders about 
this. If we suppose that, by assuming a suitable expansion 
of the measuring rods at each point, one could simplify the 
geometry, why not assume such an expansion? This, 
however, would lead us to developments lying outside the 
scope of the book. 

NATHAN ROSEN 
University of North Carolina 
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Vector Analysis. Eart C. Rex. Pp. 88+-viii, Figs. 45, 
84X11 in. Lithographed. Wm. C. Brown Company, 
Dubuque, 1952. Price $3.25. 


This text covers vector algebra, dyadics, and vector 
calculus. Fifty-three illustrative problems are worked out, 
and there are numerous exercises, many of a routine 
computational nature. No answers are given. At the end of 
each chapter is a lengthy summary, which, in some cases, 
includes additional discussion and new proofs. 

The third chapter contains two innovations, both ap- 
parently original with the author. The first is the use of 
vector equations to describe all of the points interior to 
and on a surface. Applications of this technique are given 
in both the third and fifth chapters. The other innovation 
is the definition of ‘‘dot’”’ and “‘cross’’ products of complex 
numbers. 

The author appears to have put a great deal of thought 
into many parts of the book, particularly on vector algebra. 
In spite of this, the book has flaws, of which the most 
serious are numerous errors and frequent lack of rigor. 
In his derivation of the scalar and vector products of two 
vectors in component form from the definitions, the author 
includes every step, except that he fails to mention (let 
alone prove) that these products are distributive with 
respect to addition, and that his derivation depends on 
this fact. Most of the proofs are given in great, not to say 
excessive, detail. A notable exception is the divergence 
theorem, which the author considers ‘‘almost axiomatic.” 
This might, perhaps, be justifiable, had the integral 
definition of divergence been used. Divergence is, however, 
defined using the del operator in differential form, and the 
discussion of the physical meaning of divergence is obscure 
and insufficient to render the divergence theorem obvious. 

The following are typical of the errors which occur. 
On page 19, the vector equation of a hyperbola, in the 
standard position, is given as r=aicos#—bjsin@. In 
illustrative example 43, the direction of the magnetic 
field of a long straight wire carrying a current is taken to 
be radially outward from the wire. In illustrative example 
51, the author obtains the value zero for /“.dr/(1+r)? 
because he overlooks the singularity of the integrand at 

= —1. It is unfortunate that naiveté of this sort permeates 
the entire work. 


The author uses unusual terminology in many places, 
and frequently introduces terms without definition. The 
most glaring instance of the latter is the use in Chapters I 
and II of the word ‘‘component,’’ which is not defined 
until the end of Chapter III. The discussion is frequently 
obscure, and the wording ranges, on occasion, from clumsy 
to ungrammatical. There are a number of misprints. 


Mary H. PAYNE 
Michigan State College 
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American Physical Society 
New England Seciton 


The fall meeting of the New England Section of the 
American Physical Society was held at the University of 
New Hampshire on October 25, 1952. The following pro- 
gram was presented: 


Welcoming remarks. President RoBErRT F. CHANDLER, 
Jr., University of New Hampshire. 


Symposium on High Energy Physics 


Design considerations for very high energy accelerators. 
M. S. Livincston, Brookhaven National Laboratory and 
Massachusetts Institute of Technology. 

New particles in cosmic rays. R. W. WILLIAMS, Massa- 
chusetts Institute of Technology. 

Interaction of protons with free and bound neutrons in 
the 100-Mev region. Kart Straucn, Harvard University. 


Address of Retiring Chairman 
Some optical filters. Nora MouLeErR, Smith College. 


Contributed Papers 


1. An inexpensive, easily controlled process for de- 
veloping 400-micron nuclear emulsions. Mary B. SUMMER- 
FIELD, University of California, Berkeley. 

2. Quantum-mechanical methods in classical physics. 
Davin L. FaLkorF, Brandeis University and Massachusetts 
Institute of Technology. 

3. Activation energies of the selenium-tellurium alloys. 
Harry H. Hatt, DANIEL T. HEDDEN, and THomas J. 
TURNER, University of New Hampshire. 

4. Crystal diodes in the electronics laboratory. IRVING 
L. Korsky, Smith College—Some experiments using ger- 
manium diodes, suitable for use in the undergraduate 
electronics laboratory, were described. Among these are 
the familiar tubeless unpowered crystal radio; a simple 
clipper circuit for producing trapezoidal waves, which are 
used in an experiment designed to illustrate the response 
of several circuit elements to pulses; and simple voltage- 
and frequency-multiplier circuits. 

The question of where in the traditional electronics 
course to introduce the newly developed semiconductor 
substitutes for vacuum tubes was discussed. 

W. My, Preston, Secretary 


The Tape Recordings of Important Speeches 


E. U. Connon, K. K. Darrow, E. Fermi, J. C. SLATER 


Free tape recordings are available on loan to members 
of the Association who wish to use them for presentation 
to classes, seminars, or science clubs. The subjects are 
four of the six lectures presented as a Symposium on 
Physics Today at the Twentieth Anniversary meeting of 
the American Institute of Physics, Chicago, Illinois, on 
October 25, 1951, as follows: 


The atom. E. U. Conpon. 

Physics as science and art. K. K. DARRow. 
The nucleus. E. FERMI. 

The solid state. J. C. SLATER. 


Professor Slater’s speech has been added to the list 
since our first announcement. 


Recordings are now available in two forms: 


1. Double track. Each speech is complete on one reel, 
which must be played at 3}inches per sec. No 
rewinding is necessary. Track No. 1 is played to 
the end, then the reels are reversed (and turned 
over) to play track No. 2. This procedure rewinds 
the tape into its original condition. 


2. Single track. Two 7-inch reels are required for 
each address. They must be played at 7} inches 
per sec, and, of course, rewound to restore them 
to their original condition. 


Requests for recordings should be made to THomas H, 
Oscoop, Editor, American Journal of Physics, Michigan 
State College, East Lansing, Michigan. To allow some 
flexibility in scheduling, a request should specify a pre- 
ferred date and two alternative acceptable dates, as well 
as the type of recording (single or double track, or either). 

No charge is made for the loan of recordings. The bor- 
rower is expected, however, to pay the return postage to 
the Editor’s office or to any other address the Editor 
specifies. Transportation is to be by first class mail, special 
delivery, which will cost less than one dollar per reel at 
present rates. Reels will be mailed in convenient reusable 
boxes. 

Prospective borrowers may consult the published papers 
in Physics Today, November, 1951, January, 1952, and 
March, 1952. 
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OUTSTANDING 


NEW SECOND EDITION 


CO L L ie G ie ROBERT L WEBER 


e MARSH W. WHITE 


o ‘onl YS i CS KENNETH V. MANNING 


The Pennsylvania State College 
901 pages, $6.50 


Highly acclaimed and widely adopted in the first edition, this outstanding text has been completely revised. 
Changes in the new edition reflect the experience gained by the authors in their use of the original edition as 


well as the suggestions of scores of teachers who have found College Physics an invaluable foundation for their 
courses. 


This text presents the basic ideas of physics for students of science and engineering at the college level. Simple, 
direct, and concise explanations, and solved numerical examples combine to give him an exact knowledge of 
physical principles. Derivations using calculus notation are given in the Appendix for those to whom they will 
be helpful. Facility in the application of physical principles is encouraged by discussion questions, careful treat- 
ment of units, proper use of significant figures, and a pairing of problems with answers provided for the even- 
numbered. The vitality of physics is emphasized especially in chapters on atomic and nuclear physics and in 
the biographical notes and portraits of all the Nobel prize winners in physics. 


Send for a copy on approval 


McGRAW-HILL BOOK COMPANY, Inc. 


330 West 42nd Street New York 36, New York 
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